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Abstract: We provide a complete classification of asymptotic quasinormal frequencies for static, spher- 
ically symmetric black hole spacetimes in d dimensions. This includes all possible types of gravitational 
perturbations (tensor, vector and scalar type) as described by the Ishibashi-Kodama master equations. 
■^J - \ The frequencies for Schwarzschild are dimension independent, while for Reissner-Nordstrom are dimension 
dependent (the extremal Reissner-Nordstrom case must be considered separately from the non-extremal 
J . case). For Schwarzschild de Sitter, there is a dimension independent formula for the frequencies, except 
in dimension d = 5 where the formula is different. For Reissner-Nordstrom de Sitter there is a dimen- 
sion dependent formula for the frequencies, except in dimension d = 5 where the formula is different. 
_C Schwarzschild and Reissner-Nordstrom Anti-de Sitter black hole spacetimes are simpler: the formulae 
\ for the frequencies will depend upon a parameter related to the tortoise coordinate at spatial infinity, 
and scalar type perturbations in dimension d = 5 lead to a continuous spectrum for the quasinormal fre- 
quencies. We also address non-black hole spacetimes, such as pure de Sitter spacetime — where there are 
quasinormal modes only in odd dimensions — and pure Anti-de Sitter spacetime — where again scalar type 
perturbations in dimension d = 5 lead to a continuous spectrum for the normal frequencies. Our results 
match previous numerical calculations with great accuracy. Asymptotic quasinormal frequencies have also 
been applied in the framework of quantum gravity for black holes. Our results show that it is only in 
the simple Schwarzschild case which is possible to obtain sensible results concerning area quantization or 
loop quantum gravity. In an effort to keep this paper self-contained we also review earlier results in the 
literature. 
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1. Introduction 



Black holes are an undeniable landmark in the road that connects classical to quantum gravity. Having 
been first discovered as static solutions of classical general relativity, they were later shown to actually 
radiate and evaporate once quantum effects were properly taken into account ||, gj. Thus, it seems that in 
a true quantum theory of gravity, non-extremal black holes will actually be unstable. In this setting, an 
important question that immediately comes to mind is whether the black hole solution under consideration 
was really a stable solution of the classical theory, to start off with. 

Researchers first focused on analyzing the linear stability of four dimensional black hole solutions of 
general relativity in |j[ ||, ||. The linear perturbation theory for the Schwarzschild black hole was set up 
in ||, where the classical stability of the solution was also proven. The equation derived in that paper to 
describe the linear perturbations, and their frequencies, is a Schrodinger-like equation and is now known 
as the Regge-Wheeler (RW) equation. The procedure to derive the RW equation is based on a study 
of the linearized Einstein equations in the given background and proceeds with a decomposition of the 
perturbation in tensor spherical harmonics (for spherically symmetric backgrounds) in order to obtain 
a radial equation describing the propagation of linear perturbations. This procedure, as applied to the 
analysis of perturbations to the Schwarzschild metric, was brought to firmer grounds in Q with an explicit 
construction of four dimensional tensor spherical harmonics. This work was further completed in || with 
an extension to the Reissner-Nordstrom (RN) black hole solution. There, it was shown that the complete 
Einstein-Maxwell system for each type of multipole (electric or magnetic) could be reduced to two second 
order Schrodinger-like equations, generalizing the RW equation of the Schwarzschild case. 

Part of the physical picture that emerged from this study of linear perturbations to black holes is 
the following. After the onset of a perturbation, the return to equilibrium of a black hole spacetime is 
dominated by damped, single frequency oscillations, which are known as the quasinormal modes (we refer 
the reader to || [7j for recent reviews on quasinormal modes, and a more complete list of references). 
These modes are quite special: they depend only on the parameters of the given black hole spacetime, 
being independent of the details concerning the initial perturbation we started off with. Moreover, modes 
which damp infinitely fast do not radiate at all, and can thus be interpreted as some sort of fundamental 
oscillations for the black hole spacetime. We shall return to this point in a moment. 

It was not until recently that the black hole stability problem was addressed within a d-dimensional 
setting || These papers tried to be as exhaustive as possible, studying in detail the perturbation 

theory of static, spherically symmetric black holes in any spacetime dimension d > 3 and allowing for the 
possibilities of both electromagnetic charge and a background cosmological constant. The set of equations 



describing linear perturbations in d-dimensions was derived in |8|, 1C]. These equations generalize the RW 
equation and will be denoted as the Ishibashi-Kodama (IK) master equations. The perturbations come in 
three types: tensor type perturbations, vector type perturbations and scalar type perturbations. It should 
be noted that this nomenclature refers to the tensorial behavior on the sphere, S d_2 , of each Einstein- 
Maxwell gauge invariant type of perturbation, and is not related to perturbations associated to external 
particles. For instance, one should not confuse vector type perturbations with perturbations associated to 
the propagation of a spin-1 vector particle in the background spacetime, or scalar type perturbations with 
perturbations associated to the propagation of a spin-0 scalar particle. The IK master equations were used 
in [nj] to study the stability of cf-dimensional black holes, and although many known solutions were 
shown to be stable, stability of some other solutions is still an open problem. In this work we shall make 
use of the IK master equations to analytically compute asymptotic quasinormal frequencies, and thus focus 
on the d-dimensional Einstein-Maxwell classical theory of gravity. 
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Having thus acquired a list of stable black hole solutions, the next question to address within this 
general problem are the quasinormal modes — the damped oscillations which describe the return to the 
initial configuration. As we have said, modes which damp infinitely fast do not radiate, and they are 
known as asymptotic quasinormal modes. Besides their natural role in the perturbation theory of general 
relativity, asymptotic quasinormal modes have recently been focus of much attention following suggestions 
that they could have a role to play in the quest for a theory of quantum gravity [JO], 12]. It was suggested 
in that an application of Bohr's correspondence principle to the asymptotic quasinormal frequencies 
could yield new information about quantum gravity, in particular on the quantization of area at a black 
hole event horizon. It was further suggested in [12] that asymptotic quasinormal frequencies could help 
fix certain parameters in loop quantum gravity. Both these suggestions lie deeply on the fact that the 
real part of the asymptotic quasinormal frequencies is given by the logarithm of an integer number, a 
fact that was analytically shown to be true, for Schwarzschild black holes in d-dimensional spacetime, in 
p"3| , PI ]. A question of particular relevance that immediately follows is whether the suggestions in 
are universal or are only applicable to the Schwarzschild solution. Given the mentioned analysis of 
one has at hand all the required information to address this problem and compute asymptotic quasinormal 
frequencies for d~dimensional black holes 1 . A preliminary clue is already present in 14], where the analysis 
of the four dimensional RN solution yielded a negative answer: the asymptotic quasinormal frequencies 
obeyed a complicated relation which did not seem to have the required form. Another clue was presented 
in |i~6f| , where the analysis of the four dimensional Schwarzschild de Sitter (dS) and Schwarzschild Anti-de 
Sitter (AdS) black holes also yielded a negative answer; again the asymptotic quasinormal frequencies did 
not seem to have the required form. It is the goal of this paper to carry out an extension of the techniques 
in [14, [16|] to static, spherically symmetric black hole spacetimes in any dimension d > 3, including both 
electromagnetic charge and a background cosmological constant. Besides the intrinsic general relativistic 
interest of classifying these asymptotic quasinormal frequencies, it is also hoped that our results can yield 
conclusive implications for the proposals of [|ll], 12], dealing with the application of quasinormal modes to 
quantum gravity. We shall later see that it is only in the simple Schwarzschild case which is possible to 
obtain sensible results concerning area quantization or loop quantum gravity. 

It is important to stress that even if the ideas in [11, 12] turn out not to be universal, it is still the case 
that quasinormal frequencies will most likely have some role to play in the quest for a theory of quantum 
gravity. Indeed, quasinormal frequencies can also be regarded as the poles in the black hole greybody 
factors which play a pivotal role in the study of Hawking radiation. Furthermore, the monodromy technique 



introduced in |T4| to analytically compute asymptotic quasinormal frequencies was later extended, in [17], 
so that it can also be used in the computation of asymptotic greybody factors. It was first suggested in 
[17] that the results obtained for these asymptotic greybody factors could be of help in identifying the dual 
conformal field theory (CFT) which microscopically describes the black hole, and these ideas have been 
taken one step forward with the work of fief . It remains to be seen how much asymptotic quasinormal 
modes and greybody factors can help in understanding quantum gravity. 

The organization of this paper is as follows. In section 2 we provide a summary of our results for the 
reader who wishes to skip the technical details on a first reading. Section 3 represents the main body of 
the paper, where we use both the analytic continuation of the IK master equations to the complex plane 
and a method of monodromy matching at the several singularities in the plane, in order to analytically 
compute asymptotic quasinormal frequencies for static, spherically symmetric black hole spacetimes in 
dimension d > 3. This includes a brief review of quasinormal modes, as well as the Schwarzschild case 
of p4| , for both completeness and pedagogical purposes. In section 4 we address some exact solutions 

lr The least damped quasinormal modes associated to the IK master equations were addressed numerically in Jl^ ]. 
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for quasinormal frequencies, dealing with non-black hole spacetimes. We address Rindler, dS and AdS 
spacetimes, providing full solutions in all cases 2 . It should be noted that there has been some confusion in 
the literature concerning the computation of quasinormal modes using the monodromy technique, as well 
as on the computation of quasinormal modes in pure dS spacetime. We hope that this paper will serve 
to lay these confusions to rest. Section 5 reviews [|ll], 12] and the applications of asymptotic quasinormal 



frequencies to area quantization and loop quantum gravity. We show that these applications only seem 
to work in the simple Schwarzschild setting. In section 6 we conclude, listing some future directions of 
research. We also include three appendices. In appendix A we present a thorough list of conventions for 
black holes with mass M, charge Q and background cosmological constant A, in d-dimensional spacetime. 
Appendix B includes all required formulae for the IK master equations, and appendix C makes a complete 
study of the tortoise coordinate in the spacetimes in consideration, alongside with analysis of the IK master 
equation potentials at several singularities in the complex plane. Throughout we will show that our results 
match earlier numerical calculations with great accuracy and will review some of the earlier literature on 
each case considered, as an effort to make this paper self-contained. 

2. Summary of Results 

For the reader who wishes to skip the main calculation on a first reading, we present in the following a 
summary of our results for the asymptotic quasinormal frequencies, in any spacetime dimension d > 3. In 
the appendices we review the black hole spacetime solutions we wish to consider, as well as the perturbation 
theory for these spacetimes which leads to the quasinormal mode analysis. The set of equations describing 
quasinormal modes in d-dimensions was derived in H [Tc[] , and the perturbations come in three types: 
tensor type perturbations, vector type perturbations and scalar type perturbations. For each spacetime 
in consideration we compute asymptotic quasinormal frequencies, given each type of <i-dimensional per- 
turbation. In the case of black hole spacetimes the computation involves a detailed monodromy analysis, 
alongside with some simple differential equations. In fact, the IK master equation describing quasinormal 
modes is of Schrodinger type, where the potential is associated to the perturbation under study. When 
using the monodromy technique one needs to solve the Schrodinger-like master equation at points where 
the potential is usually either zero (yielding simple plane wave solutions) or of Bessel type 

-^(x) + J -^Hx)=u; 2 Hx), (2.1) 
ax z 4ar 

where j is determined for each different case (for each type of perturbation and for each background 
considered). This equation can be solved in terms of Bessel functions, J u (x), with the result 

$(x) = A + V 27nox J L (lux) + A_ V 27nox J_ L (ujx) , (2.2) 

2 2 

where A + and A_ are constants. Let us list the values of j which one has to deal with, when considering the 
black hole singularity region, at the origin of the coordinate frame. For all uncharged black hole solutions 
one finds jj = for tensor type perturbations, jy = 2 for vector type perturbations and js = for scalar 
type perturbations. For all charged black hole solutions one finds jj = for tensor type perturbations, 
jv+ = jy- = f° r vector type perturbations and j$+ = jg- = ^3 f° r scalar type perturbations. It is 
furthermore simple to observe that for all black hole spacetimes j' v + = j v - — iv> is+ = is- — Js, Jt = is 7 
and jj + jv = 2. This will ultimately imply that the asymptotic quasinormal frequencies will depend only 



2 Other exact calculations of quasinormal frequencies can be found in, e.g., J^, |5c| ]. 
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on the background, being the same for the three types of perturbations. Let us also recall that one can 
prove (for spacetimes which are not asymptotically AdS), without explicit computation of the quasinormal 



frequencies, that j = and j = 2 perturbations must have identical quasinormal spectra [21]. 

For spacetimes without a black hole (but with cosmological constant) the calculation of quasinormal 
frequencies follows in an analytic way, without any asymptotic restrictions. Here one need not use the 
monodromy method (even though it can be simply generalized to those cases as well) and we have proceeded 
by solving the wave equation directly. This is accomplished by first finding an appropriate change of 
variables that brings the quasinormal master equation to a hypergeometric form. Let us finally list our 
results on quasinormal frequencies: 

The Schwarzschild Solution: This case was first studied in |l4|] and we address it in this paper for the 
sake of completeness. For all types of perturbations, tensor, vector and scalar type perturbations, the 
algebraic equation for the asymptotic quasinormal frequencies is the same and is 

e T7F + 3 = 0, 

where Th is the Hawking temperature in the Schwarzschild spacetime. As is well known, this case is 
particularly simple and one can moreover solve for the asymptotic quasinormal frequency as 

u = T H log3 + 2mT H (n+ ~J (n G N, n > 1) . 
The result is independent of spacetime dimension 3 . 

The RN Solution: This case was studied, in the particular d = 4 case, in |l4j . Here we extend those results to 
arbitrary dimension. For all types of perturbations, tensor, vector and scalar perturbations, the algebraic 
equation for the asymptotic quasinormal frequencies is the same and is 



e T H +|^l + 2 cos (vrj) J + [2 + 2 cos (vrj) je t h = 0, 
where are the Hawking temperatures at outer and inner horizons (notice that Tj^ < 0), and where 



d- 3 

J 



2d -5 

There is no known algebraic solution in uj for the above equation. 

The Extremal RN Solution: It is important to realize that, in general, quasinormal frequencies of extremal 
solutions cannot be obtained from the corresponding expression for the non-extremal solution. In fact, the 
monodromy technique deployed in this paper is very sensitive to both the location of the complex horizons 
and the structure of the tortoise at the origin. Thus, as one changes the background solution there will 
be a change of topology in the complex plane and the solution to the quasinormal mode problem will be 
different. We present the extremal RN solution as an example, but one should keep this in mind if also 
interested in extremal solutions with a cosmological constant (which we do not address in this paper, but 
list the possibilities in appendix A). For all types of perturbations, tensor, vector and scalar perturbations, 
the algebraic equation for the asymptotic quasinormal frequencies is the same and is 



3 By "independent" we mean that there is no explicit dependence on the spacetime dimension, d, in the above formula. Of 
course if one wishes to compute the Hawking temperature in this Schwarzschild background, then one finds that it relates to 
the mass M via an expression which also involves the dimension d. 
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sin(f) e*-sm(^)=0, 

where j is as in the previous RN case and T is not a temperature — in fact in the extremal case there is no 
Hawking emission. Rather, it is given by 

„ d-3 ( d-3 \ 



d-3 



where n is related to the black hole mass (see appendix A). It is simple to solve for the asymptotic 
quasinormal frequency as 




(n € N, n > 1) 



Observe that in dimension d = 5 there is no solution for uj. This is in fact the only dimension where there 
is no solution for the asymptotic quasinormal frequencies of the extremal RN geometry. 



The Schwarzschild dS Solution: This case was studied, in the particular d = 4 case, in [16|. Here we extend 
those results to arbitrary dimension. For all types of perturbations, tensor, vector and scalar perturbations, 
the algebraic equation for the asymptotic quasinormal frequencies is the same and is 

f UJ UJ \ f UJ UJ \ 

cosh — — + 3 cosh — — + — — = 0, 

\2T H 2T C ) \2T H 2T C ) 

where Tjj is the Hawking temperature at the black hole event horizon and Tc is the [negative] Hawking 
temperature at the cosmological horizon. There is no known algebraic solution in uj for the above equation. 
The result is independent of spacetime dimension, except in dimension d = 5 where the formula above 
must be replaced by: 

f UJ uj \ ( UJ UJ \ 

sinh — ■ — — — 3sinh — ■ 1 — — — = 0. 

\2T H 2T C ) \2T H 2T C ) 

Observe that for this solution one can actually take the Schwarzschild limit A — > without provoking 
any topology change in the complex plane where the monodromy analysis is performed. Thus, the result 
obtained for Schwarzschild dS includes the pure Schwarzschild solution once one sets the cosmological 
constant to vanish (both expressions have the same, correct, Schwarzschild limit). 

The RN dS Solution: For all types of perturbations, tensor, vector and scalar perturbations, the algebraic 
equation for the asymptotic quasinormal frequencies is the same and is 



cosh (2% - 4) +(1 + 2 cos (7rj)) cosh + 4) +(2 + 2 cos cosh (i + 2% + 4) = °' 

where j is as in the previous RN case, Tjj are the Hawking temperatures at outer and inner black hole event 
horizons and Tq is the Hawking temperature at the cosmological horizon. There is no known algebraic 
solution in uj for the above equation. While this result explicitly depends on the spacetime dimension 
(because j depends on d), the formula is not valid in dimension d = 5, where it must be replaced by: 
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. ( uj w \ 1 + y/E . , / u uj \ 3 + a/5 . , / u lj u) \ 

smh — j- H smh — - + — — H smh 1 + — — = 0. 

\2T C 2T+7 2 V2T+ 2T c ; 2 2T+ 2T c y 

Similarly to the Schwarzschild dS case, one can take the pure RN limit A — > without provoking any 
topology change in the complex plane where the monodromy analysis is performed. Thus, the result 
obtained for RN dS includes the pure RN solution once one sets the cosmological constant to vanish (both 
expressions have the same, correct, RN limit). 

The Schwarzschild AdS Solution: This case was studied, in the particular case of d = 4 and large black hole, 



in [p.q| . It was further studied, in the particular case of d = 5 and large black hole, in [|22j, |23|, [2j, |25] 
using a variety of different analytical methods. Here we extend those results for arbitrary dimension and 
away from the large black hole approximation. For all types of perturbations, tensor, vector and scalar 
perturbations, the algebraic equation for the asymptotic quasinormal frequencies is the same and is 



tan ^— (d + 1) — ujxq^ 



i 

3' 

where xq is a parameter related to the tortoise coordinate at spatial infinity and is given by 



where the R n are the (d — 1) complex horizons and the k n the surface gravities at each of these complex 
horizons. There is no general analytic solution for xq. However, in some cases, it can be computed exactly. 
For instance, for large black holes we compute it to be 

I / 7r \ ( in 
4Th sin I — I exp 



Xq \d — 1 J \d—l 

where T# is the Hawking temperature in the Schwarzschild AdS spacetime. In spite of not having a general 
analytic solution for xq one can still solve for the asymptotic quasinormal frequency as 

7T f %\ 

ujxq = — (d + 1) — arctan ( - j + nir (n G N, n S> 1) . 

If one concentrates on large Schwarzschild AdS black holes, of particular relevance to describe thermal 
gauge theories within the AdS/CFT framework, then the formulae above lead to the following analytical 
result for the leading term, as n — > +oo, in the asymptotic quasinormal frequencies: 

For the most popular AdS/CFT dimensions, d = 4, d = 5 and d = 7, one obtains 

uJd=i ~ 2v SirnTiie 3 , w<f=5 ~ 2v 2%nTne 4 and u>d=7 ~ InnTue s . 

Notice that the above formulae are not valid for scalar type perturbations in dimensions four and five. 
Instead, one finds for these perturbations in dimension d = 4 
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3"7T / i \ . 

ojxq = — — arctan I - I + nir (n £ N, n ~S> 1) , 

and more surprisingly, a continuous spectrum in five dimensions, 

The RN AdS Solution: For all types of perturbations, tensor, vector and scalar perturbations, the algebraic 
equation for the asymptotic quasinormal frequencies is the same and is 



sm 



where xq is as before (only one should recall that this time around there are (2d — 4) complex horizons as 
we are in a charged situation) and j is as in the previous RN cases. Again, there is no general analytic 
solution for xq. In spite of this, one can still solve for the asymptotic quasinormal frequency as 

ojxq = y (d + 1) - ^log ( 2 cos ( ) ) + nvr (n G N, n > 1) . 



Notice that the above formula is not valid for scalar type perturbations in dimensions four and five. 
Instead, one finds for these perturbations in dimension d = 4 

37T 1 j ( 71 J \ \ 

lox = — — - log ( 2 cos ( — j j + n-K (n € N, n > 1) , 
and more surprisingly, a continuous spectrum in five dimensions, 

uj G C. 

Besides the previous results, concerning black hole spacetimes, we have also addressed the case of 
spacetimes without a black hole but with a cosmological constant, i.e., the cases of AdS (where there are 
only normal modes) and of dS (where one finds quasinormal modes only in odd spacetime dimensions). 
Let us list those results as well: 



The AdS Solution: This case was studied before in |2(| for a massless scalar field, and in [^] for d = 4 tensor 
and vector type perturbations. Here we extend those results to arbitrary dimension and perturbation. The 
first thing to notice is that AdS spacetime acts as an infinite potential well and the Schrodinger-like 
equation yields real frequencies only. In other words, there are no quasinormal modes in pure AdS, only 
normal modes. The normal frequencies one finds are (no asymptotic restrictions here, this is an exact 
result) 



u n = (2n + d + l- j) , n £ N, 

where I is the angular momentum quantum number (eigenvalue of the spherical laplacian) and j is j = 1 
for tensor type perturbations, j = 2 for vector type perturbations and j = 3 for scalar type perturbations. 
Thus, for fixed £, tensor and scalar type perturbations yield the same normal frequencies, although different 
from the vector type frequencies. Notice that the above formula is not valid for scalar type perturbations 
in dimensions four and five. Instead, one finds for these perturbations in dimension d = 4 
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= ^J\\\{2n + l + 2), 



and more surprisingly, a continuous spectrum in five dimensions, 



The dS Solution: This case was studied in several different works, discussed in the main text, all with con- 
tradictory results. Here, we solve this apparent confusion in the literature with the following results. There 
are no dS quasinormal modes when the spacetime dimension is even. However, when the spacetime di- 
mension is odd, there are quasinormal modes, with the quasinormal frequencies (no asymptotic restrictions 
here, this is an exact result) 



where £ is the angular momentum quantum number (eigenvalue of the spherical laplacian) and j is j = 1 
for tensor type perturbations, j = 3 for vector type perturbations and j = 5 for scalar type perturbations. 
Thus, for fixed I, tensor and scalar type perturbations yield the same quasinormal frequencies, although 
different from the vector type frequencies. The above formula is independent of spacetime (odd) dimension. 

Finally, we have studied the Rindler solution where we found — without great surprise — that there are 
no quasinormal modes. A detailed discussion concerning applications of these results to quantum gravity 
is included in the text, and we shall make no attempt of summarizing it in here. 

3. Asymptotic Quasinormal Frequencies 

In this section we first wish to review and set notation on both the perturbation theory for spherically 
symmetric, static ci-dimensional black holes 4 , with mass M, charge Q and background cosmological con- 
stant A, and the computation of quasinormal modes and quasinormal frequencies. We refer the reader to 
appendix A for a full list of conventions on the black hole spacetimes we shall consider. In the general 
(M, Q, A) case, there are two different types of fields which can be excited: these are the electromagnetic 
vector field and the gravitational metric tensor field g^ v , and we shall study perturbations to both 
these fields. Because there is no scalar field present in the Einstein-Maxwell system, there are no scalar 
field perturbations to consider. Nevertheless, we shall start by studying the scalar wave equation in our 
black hole backgrounds, in order to set notation on quasinormal modes. 

Consider a massless, uncharged, scalar field, <f>, in a background spacetime described by a metric g^v 
Its wave equation is well known 



where g is the determinant of g^. The question we wish to address is what is the form of the wave equation 
for a background spacetime metric of the type 





(3.1) 



g = -f(r) dt®dt + f(r) 1 dr ® dr + r 2 dft 2 d 



d-2i 



4 Throughout this work we consider only dimension d > 3. 
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i.e., in a d-dimensional spherically symmetric background. This issue was addressed in [28], with the fol- 
lowing result. First perform a harmonic decomposition of the scalar field as <fi = Y^e m r_5_ ^( r ; t) ^im 
where the 6i are the (d — 2) angles and the Y^ m {9i) are the ci-dimensional spherical harmonics. Then, if 



^(r,t) = * w (r)e iw * 

is the [time] Fourier decomposition of the scalar field, the wave equation can be recast in a Schrodinger-like 
form as 

+ V ( x )*»( x ) = ^M*). ( 3 - 2 ) 

where x is the so-called tortoise coordinate and V{x) is the potential, both determined from the function 
f(r) in the background metric. The tortoise coordinate is defined so that || 0] 

f(r)— ( f(r)±) = *- 
dr \ dr J dx 2 

and is thus given by 

dx = ^- (3.3) 

This new coordinate x keeps infinity (or the cosmological horizon, Rc, in the dS case) at x = +oo and 
sends the black hole event horizon, Rh, to x = — oo (in the charged cases this refers to the outer horizon). 
The region of positivity of f(r) thus becomes the real line when in tortoise coordinates: f(x) > for 
The potential in the Schrodinger-like equation can be determined as one moves from the general form of 



the wave equation, (3.1), to its Schrodinger-like form, (|3.2j) , and is given by [28] 



I ( , , /, , , H£ + d-3) + (d-2) (d-4)/(r) + {d - 2) f{r) 



It 

Here, t (£ + d — 3) (with £ € N) are the eigenvalues of the Laplacian on the E> d ~ 2 sphere, and the potential 
still needs to be re-written in terms of the tortoise coordinate in order to be used in the Schrodinger-like 
equation. Once all this is done, the question still remains on what are the allowed values for u, i.e., what 
is the spectrum of the Schrodinger-like operator above. It turns out that the spectrum will contain both a 
continuous and a discrete part, this last one being found when imposing "out going" boundary conditions: 
nothing arrives neither from infinity nor from within the black hole horizon. The [spherical] waves are out 
going at both extrema in x: 



3>u(%) ~ e luJX as x —> -oo, 
$ u (x) ~ e~ %u)X as n +oo. 

Solutions to the Schrodinger-like equation with the previous "out going" boundary conditions lead to a 
discrete set of allowed frequencies, the quasinormal frequencies, with corresponding solutions, the quasi- 
normal modes (6|, f?]]. These quasinormal frequencies, uj, are complex numbers, the real part representing 
the actual frequency of oscillation, the imaginary part representing the damping. It can be further shown 
for many cases that frequencies of quasinormal modes with negative imaginary part do not exist, meaning 
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that these solutions are actually stable. Our interest here is to study the asymptotic behavior of quasi- 
normal frequencies, i.e., the case where the imaginary part of co grows to infinity. It turns out that in 
some cases the real part of the frequency approaches a finite limit, while the imaginary parts grow linearly 
without bound. These generics on quasinormal modes will obviously hold true for the vector field pertur- 
bations and metric tensor field perturbations, the main difference being the change in the potential of the 
Schrodinger-like equation. A complete description of the required potentials to describe the most general 



situation, as derived in ||, 10], is presented in appendix B, to which we refer the reader for further details. 



For a four dimensional Schwarzschild black hole, one has the asymptotic quasinormal frequencies 

1 



lim uj n ~ [offset] + in [gap] + O 

where the real part of the offset is the frequency of the emitted radiation, and the gap are the quantized 
increments in the inverse relaxation time. Here, the gap is given by the surface gravity. One can try 
to extend this analysis to more general situations and also include spacetimes with two horizons, but 



then generic results become much harder to obtain [^, 30, 31, [3^] (curiously, for spacetimes with multiple 
horizons, there is a unique definition of temperature only when the ratio of surface gravities is in Q |32| , j33|). 
Another property of quasinormal modes is a reflection symmetry uj <-> — U) which changes the sign of Re(cc>). 
Indeed, is a quasinormal mode corresponding to —uj. 

We have chosen the time dependence for the perturbation to be e lujt , so that Im(u>) > for stable solu- 
tions (implying that the perturbation vanishes exponentially in time). This implies that Re(±zux) — > +oo 
as x — ► =Foo and thus, while solutions can oscillate, they must exponentially increase with In other 
words, & w (x) is not a normalizable wave function and this is, ultimately, the reason why we speak of quasi- 
normal modes rather than normal modes (which would form a complete set of stationary eigenfunctions 
for the IK master differential operator). 

Before embarking in the actual calculation, let us address the general strategy of the method introduced 
in [jH], which begins with the question of how to properly define and impose quasinormal boundary 
conditions. This is known to be somewhat complicated, at least at an operational level, because as long 
as we restrict x S R the quasinormal boundary conditions above amount to distinguishing between an 



exponentially vanishing term and an exponentially growing term. The idea of [14] is to do an analytic 
continuation to the complex plane, taking both r € C and x £ C. We will see that in the complex plane 
one can impose quasinormal boundary conditions in a completely new way. Indeed, if one picks the complex 
contour Im(ux) = in C, then the asymptotic behavior of e ±luJX is always oscillatory on this line and there 
will never be any problems with exponentially growing versus exponentially vanishing terms. One should 
thus restrict to studying the boundary conditions on the so-called Stokes line, Im(ojx) = 0. 

For Schwarzschild, RN, Schwarzschild dS and RN dS black hole spacetimes, numerical tests generically 
indicate that the asymptotic quasinormal frequencies are such that Im(a;) 3> Re(u;). In other words, one 
should have u = ujr + ma;/, with ujr, ujj £ R, ujj > and n —* +oo, in the asymptotic case. Thus, in the 
n ~ +oo limit, the contour Im(ujx) = can be approximated by the curve Re(x) = which is immediate to 
plot in C. This replacement of contours is selecting asymptotic conditions for the quasinormal modes. For 
Schwarzschild AdS and RN AdS things are different, as numerical tests generically indicate that asymptotic 
quasinormal frequencies behave as Im(a;) ~ Re(w), and thus one should have uj = n(ujR + iujj) + u>o, with 
ur,uji G R, u>i > and n — > +oo, in the asymptotic case. We shall later see with greater detail how to 
plot the Stokes lines for each black hole spacetime under consideration. 

To fully understand the advantage of the analytic continuation in the exact computation of asymptotic 
quasinormal frequencies, let us first address regions of the complex r-plane around an event horizon (be 
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it a black hole horizon or a cosmological horizon). The horizons themselves are defined by the zeroes of 
/(r), i.e., /(Rh) = 0, and besides the physical real horizons Rh £ R there can be other, non-physical 
complex horizons Rh £ C. We shall denote horizons which are not real as fictitious horizons. Power series 
expansion near any of these horizons simply yields f(r) ~ (r — Rh) /'(Rh) + • • • > and it follows for the 
tortoise coordinate 

f dr f dr 1 1 1 

X = JW)~J (r-RH)f'(RH) = 7(RH) l0g(r - Rh) = 2kn- bg(r - RH) = 4^ l0g(r " W> (3 ' 4) 

locally near the chosen horizon (provided that the horizon is nondegenerate, i.e., that Rh is a simple zero 
of f(r)). Here kn is the surface gravity and Th is the Hawking temperature. 

One learns that around any nondegenerate horizon the tortoise coordinate will thus be multivalued 
and it makes sense to ask for the monodromy of tortoise plane waves in clockwise contours around any 
given horizon f(Rn) = 0. Let us consider an horizon Rh £ C, and a clockwise contour 7 C C centered 
at Rh and not including any other horizon. As we take the tortoise coordinate around 7 it increases by 
— y~. This immediately implies the monodromy of the plane waves along the selected contour: 

^i,R„ [e ±iulX ] =e ±f ^. 

This result is quite interesting as it now allows one to recast quasinormal boundary conditions as mon- 
odromy conditions [ |14| ; if one wants, as quasinormal mode monodromy conditions. At the black hole event 
horizon the quasinormal boundary condition is 

$u(x) ~ e luJX as x — > -00, 

where x — * — 00 as r — ► Rh, with Rh the black hole horizon. One immediately re-writes this boundary 
condition as a monodromy condition for the solution of the master equation at the black hole event horizon: 

WI^ Rh [* w (x)] = e ^, 
with 7 a clockwise contour. The other quasinormal boundary condition, 

$ u (x) ~ e~ luJX as x — ► +00, 

lives at r ~ +00, as x — > +00 when r — > +00 (for all spacetimes except asymptotically dS spacetimes). 
Thus, this boundary condition cannot be recast as a quasinormal monodromy condition. Instead, if one 
considers asymptotically dS spacetimes, the above boundary condition is located at the cosmological event 
horizon, as x — > +00 when r — > Re- Then, it is immediate to recast this boundary condition as a 
monodromy condition for the solution of the master equation at the cosmological event horizon: 

9rt 7l n c [* w (a:)] =e k c, 

with 7 a clockwise contour. We shall see in the following how these simple ideas about boundary conditions 
will allow for analytic calculations of asymptotic quasinormal frequencies in all static, spherically symmetric 
black hole spacetimes. 
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3.1 Asymptotically Flat Spacetimes 

pl discusses the stability of black holes in asymptotically flat spacetimes to tensor, vector and scalar 



perturbations. For black holes without charge, all types of perturbations are stable in any dimension. For 
charged black holes, tensor and vector perturbations are stable in any dimension. Scalar perturbations are 
stable in four and five dimensions but there is no proof of stability in dimension d > 6. As we work in 
generic dimension d we are thus not guaranteed to always have a stable solution. Our results will apply if 
and only if the spacetime in consideration is stable. 

3.1.1 The Schwarzschild Solution 

For completeness, we first present a computation of the asymptotic quasinormal frequencies for the 
Schwarzschild <i-dimensional black hole, using the monodromy method. This calculation was first done in 



[14]. The interesting result shown in [14] is that, if one is only interested in the asymptotic quasinormal 
modes, there is no need to solve the IK master equation exactly. Rather, there is a method which explores 
the analytic continuation of the master equation to the complex plane and demands only for approximate 
solutions near infinity, near the origin, and near the black hole event horizon. Knowledge of the solutions 
in these regions, together with monodromy matching along a specially chosen contour, then yields the 
quasinormal frequencies. Let us carefully explain this method in the simplest Schwarzschild example, as 
we shall employ it several times in the following. 

We consider solutions of the Schrodinger-like master equation ( fT^ ) 

d 2 $> 



(x 



+ V[t(x)]$(x) = w 2 <S>(x) 



dx 2 

in the complex r-plane. Let us begin at infinity. Since the potential V{r) vanishes for r ~ +oo, we will 
have 

$(x) ~ A + e iujx + A^e~ iuJX 
in this region. The boundary condition for quasinormal modes at infinity is then 

A+ = 0. (3.5) 
Next we study the behavior of $(x) near the singularity r = 0. In this region, the tortoise coordinate is 

r d-2 

X ~ 



2(d- 2)/i' 

and the potential for tensor and scalar type perturbations is 



V[r(x)] 



1 j 2 — 1 



Ax 2 Ax 2 ' 

with j = (see appendix^]). The Schrodinger-like master equation approximates to 

d 2 $ , . j 2 - 1 , 2 , , , 

whose solution is (this is the solution for j ^ — more on this in a moment) 

TTUJX J j (LUX) , 
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where J u (x) represents a Bessel function of the first kind and B± are (complex) integration constants. One 
would next like to link this solution at the origin with the solution at infinity. 

For the Schwarzschild asymptotic quasinormal modes one has Im(u;) 3> Re(w), with Im(o;) — > +00, 
and hence u is very large and approximately purely imaginary. Consequently, one has ux G R for x G «R; 
in a neighborhood of the origin, the above relation between x and r tells us that this happens for 



rp — p g 2(d—2) d—2 ^ 

with p > and n = 0, 1, . . . , 2d — 5. These are half-lines starting at the origin, equally spaced by an angle 
of ^2 . Notice that the sign of ux on these half-lines is (— l) n ; in other words, starting with the half-line 
corresponding to n = 0, the sign of ux is alternately positive and negative as one goes anti-clockwise 
around the origin. 

Precisely because we are interested in these asymptotic modes, we may consider the following asymp- 
totic expansion of the Bessel functions 

/ 2 / w f \ 
•w) ~ V ~ cos [ z ~ ~h~ ~ 4 J ' z ^ ' ^ > 

from where we learn that 



<£(x) ~ 2B + cos (wx — a + ) + 2i?_ cos (ux — a_) 

= (B + e~ ia + + B^e~ ia -) e iu)X + (B + e ia+ + B_e ia ~) e~ iulx (3.7) 

in any one of the lines corresponding to positive ux, and where we have defined 

This asymptotic expression for $ near the origin is ideal to make the matching with its asymptotic expres- 
sion at infinity. This matching must however be done along the so-called Stokes line, defined by ux G R 
(or Im(ux) = 0), so that neither of the exponentials e ±tuJX dominates the other. In this Schwarzschild case 
the Stokes line definition corresponds to x G iM. or Re(x) = 0. 

To trace out the Stokes line Re(x) = let us first observe that we already know its behavior near the 
origin. Furthermore, this is the only singular point of this curve: indeed, since x is a holomorphic function 
of r, the critical points of the function Re(x) are the zeros of 



dx 
dr 



1 



r.d-3 



f(r) r d ~ 3 - 2/i 



(i.e., r = only). We have an additional problem that x is a multivalued function: each of the "horizons" 

2m 



Rn 



exp 



d-3 



n 



n 



0,l,...,d-4, 



is a branch point. In fact, near such points one has (see (3.4)) 

1 



x ~ 



f'(Rn) 



log(r - R r , 
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contour 



branch cut 




Figure 1: Stokes line for the Schwarzschild black hole, along with the chosen contour for monodromy matching, in 
the case of dimension d = 6. 







Figure 2: Numerical calculation of the Stokes lines for the Schwarzschild black hole in dimensions d = 4, d = 5, 
d = 6 and d = 7. The different shadings also illustrate the various horizon singularities and branch cuts (note that 
these branch cuts are not necessarily equal to the ones used for the calculation in the main text). 



Thus we see that although the function Re(x) is well defined around R n with n = and n = (if d is 
odd), as f'(R n ) is real in these cases, it will be multivalued around all the other fictitious horizons. 

For r ~ oo one has x ~ r (see appendix |C|). Consequently, x is holomorphic at infinity and we can 
choose the branch cuts to cancel out among themselves. Therefore Re(x) is well defined in a neighborhood 
of infinity, and moreover Re(a;) = will be approximately parallel to Re(r) = in this neighborhood. 
Two of the 2(d — 2) branches of the Stokes line starting out at the origin must therefore be unbounded. 
The remaining can either connect to another branch or end up in a branch cut. On the other hand, it is 
easy to see that the Stokes line must intersect the positive real axis exactly in one point, greater than Rq. 
Using this information plus elementary considerations of symmetry and the sign of Re(x), one can deduce 
that the Stokes line must be of the form indicated in Figure H. These results are moreover verified by the 
numerical computation of the same Stokes lines, as indicated in Figure ||. 

Let us now consider the contour obtained by closing the unlimited portions of the Stokes line near 
r ~ oo, as shown in Figure |l|. At point A we have ujx 3> 0, and therefore the expansion (3/7) holds at this 
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point. Imposing condition (|3.5|) one obtains 

B +e - ia + +B„e~ ia - = 0. (3.8) 

For z ~ one has the expansion 

J v {z) = z u w(z), (3.9) 

where w is an even holomorphic function. Consequently, as one rotates from the branch containing point 
A to the branch containing point B, through an angle of x ~ ~ 2{d-2)n r °tates through an angle of 
3ir, and since 

\/2ne 37ri cox J,, (e 3ni cox) = e^ (1±i) V2ircox J, l (cox) ~ 2e 6ia± cos(wx - o±), 

=t 2 y =t 2 

one has (notice that e 3m cox = —cox) 



$(x) ~ 2 J B+e 6m + cos ( — cox - a+) + 2B_e &ia ~ cos (-m - a_) 
= (B + e 7ia + + B_e 7ia -) e iu3X + (B + e 5ia + + B_e 5ia ~) e~ iuJX 

at point B. As one closes the contour near r ~ oo, one has x ~ r and hence Re(x) > 0; since Im(u>) 3> 0, 
it follows that e lu3x is exponentially small in this part of the contour, and therefore only the coefficient of 
e -iuix ghQuid b e trusted. As one completes the contour, this coefficient gets multiplied by 

B+e 5ia+ + B_e 5ia ~ 
B + e ia + + B_e ia - 

On the other hand, the monodromy of e~ tuJX going clockwise around this contour is e~~ , where 

k = \fm 

is the surface gravity at the horizon. Since 

x ~ 7<k) log{r ~ Ro) 

for r ~ Rq, x will increase by — jr^^ = — T as one & oes clockwise around Rq, and hence e~ tuJX will get 
multiplied by 

-iUlf- \ 7VUJ 

e V T(Ro)J = e — ~. 
Thus the clockwise monodromy of $ around the contour depicted in Figure || is 

B+e 5ia+ + B_e 5ia - _ m 
B + e ia + + B_e ia - 6 

The important point to realize now is that one can deform this chosen contour — without crossing any 
singularities — so that it becomes a small clockwise circle around the black hole event horizon r = Rq. Near 
the horizon r = Rq, we have again V(r) ~ 0, and hence again 
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$(x) ~ C+e luJX + C-e~ lu)X . 
The condition for quasinormal modes at the horizon is therefore C_ = 0. Again using the fact that 

1 



log(r - Rq) 



f'(Ro) 

for r ~ Rq, we can restate this boundary condition as the condition that the monodromy of $ going 
clockwise around the contour should be 



Because the monodromy is invariant under this deformation of the contour, the condition for quasinormal 
modes at the horizon follows as 

B 5ia ++B _ e 5ia- 

— ; : e k = e k . (3.10) 

B + e ia + +B_e ia - v ' 

The final condition for equations Q3.8| ) and ( |3.10| ), which reflect quasinormal mode boundary conditions 
both at infinity and at the black hole horizon, to have nontrivial solutions (B+,B-) is simply 



e 5ia +e -^ 



e k 



& 



sin (2a+ + ?f) 



-Aid- 



sin (2a_ + !f) 



0. 



This equation is automatically satisfied for j = 0. This is to be expected, as for j 
J 



(3.11) 

the Bessel functions 



=C 2 



coincide and do not form a basis for the space of solutions of the Schrodinger-like master equation 
near the origin. Following [14|, we consider this equation for j nonzero and take the limit as j — ► 0. 
This amounts to writing the equation as a power series in j and equating to zero the first non-vanishing 
coefficient, which in this case is the coefficient of the linear part. Thus we just have to require that the 
derivative of the determinant above with respect to j be zero for j = 0: 



-nre 



nre 



44> 



sin(f + T) shr(f + T 

7TLJ\ . . /"ITLO 



+ 



2 cosh 



- ) — sinh 







fcos(f + T) -§«*(f + ^) 

2ttlo 



27TL, 

44> e k 



4^- 



log3 + i(2n + l)vr, n G N. 



For vector type perturbations, the potential for the Schrodinger-like master equation near the origin 
is of the form 



V[r(x)] 



1 



4rc 2 



with j = 2 (see appendix |C|). Repeating the same argument, one ends with the same equation as ( 3.11| ), 
except that now 2a± = \ ± ir, 4a± = tt ± 2ir. This equation is exactly the same as in the j = case, for 
which 2a± = ^, 4a± = tt, and consequently we end up with the precise same quasinormal frequencies. 
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In the calculation above we have obtained the asymptotic quasinormal frequencies uj with Me(u;) > 0. 
However, one knows that if uj is a quasinormal frequency then so is — uj. Consequently ^ = - log 3 + 
i(2n + l)ir (for n G N) must also be a solution for the asymptotic quasinormal frequencies. To understand 
how these comes about, and why we did not obtain them above, we recall that we imposed condition ( |3.5[ ) 
at point A, where ix — ► +00. Notice that we could instead have imposed it at point B, where ix — ► —00. 
It is in choosing one of these points that we automatically choose the sign of Me(u;). Indeed, if Me(c<j) > 
then e luJX 3> e~ l ^ x at point A and e tuJX <C e~ luJX at point B. Consequently, only at point A is it meaningful 
to impose condition (^). On the other hand, if Me(w) < then e luJX <C e~ tuJX at point A and e tuJX 3> e~ luJX 
at point B, and only at point B is it meaningful to impose condition (|3.5| ). Had we done this, we would 
have obtained the second set of asymptotic quasinormal frequencies. 

The results above have been thoroughly checked in the literature. Actually, an analytic calculation of 
the asymptotic quasinormal frequencies was first done in |l3| , in four dimensions, and then extended to d- 
dimensions in |14}| , at least for tensor type perturbations (and already using the monodromy method). The 
result was shown to be dimension independent and equal to Tjj log 3. It was already expected that these d- 



dimensional frequencies would scale linearly with the Hawking temperature, from the earlier results in [34]. 



Later, in [35], it was shown that the same result of Tjj log 3 holds for both vector and scalar perturbations. 



All these results were later subject to a detailed numerical check in j36|, 37], with fully positive results. 



3.1.2 The Reissner Nordstrom Solution 

Here we compute the asymptotic quasinormal modes of the RN (i-dimensional black hole using the mon- 
odromy method. This calculation was done for d = 4 in [O]. We consider solutions of the Schrodinger-like 
equation fl3,2|) in the complex r-plane. Since the potential V vanishes for r ~ +00, we will again have 



$(x) ~ A + e iuJX + A„e- iulx 
in this region, the boundary condition for quasinormal modes at infinity being 

A + = 0. (3.12) 
Next we study the behavior of $(x) near the singularity r = 0. In this region, the tortoise coordinate is 

r 2d-5 

X ~ 



(2d-5)0 2 ' 

and the potential for tensor type and scalar type perturbations is 



V[r(x)] ~ 



Ax 2 



with j = 2d-b ( see a PP en dix y). The solution of the Schrodinger-like equation in this region is therefore 
well approximated by 

B + V2mJx J i (ujx) +B-\/2 TTUJX J_ L (UJX) , 

2 2 

where J u (x) represents a Bessel function of the first kind and B± are (complex) integration constants. 

For the asymptotic quasinormal modes one has Im(u;) 3> Me(o'), and hence uj is approximately purely 
imaginary. Consequently, one has ujx E R for x £ iW; in a neighborhood of the origin, this happens for 
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j. — p g 2(2d— 5) 2d-5 

with p > and n = 0, 1, . . . , Ad — 11. These are half-lines starting at the origin, equally spaced by an 
angle of 2 d-b ' Notice that the sign of ujx on these lines is (— l) n , in other words, starting with the 
line corresponding to n = 0, the sign of ujx is alternately negative and positive as one goes anti-clockwise 
around the origin. 

From the asymptotic expansion ( |3. 6| ) we see that 



$(x) ~ 2B+ cos (ujx — a+) + 2i?_ cos (lux — q_) 

= [B+e~ ia + + B^e~ ia ') e iuJX + (B+e ia + + B_e ia ') e~ iwx (3.13) 

in any one of the lines corresponding to positive ujx, where again we define 

a± = j(l±j). 

This asymptotic expression for <I> near the origin is to be matched with its asymptotic expression at infinity. 
This matching must again be done along the Stokes line ujx G R (that is, x E iM), so that neither of the 
exponentials e ±lu3X dominates the other. 

To trace out the Stokes line Re(x) = we observe that we already know its behavior near the origin. 
Furthermore, this is the only singular point of this curve: indeed, since x is a holomorphic function of r, 
the critical points of the function Re(x) are the zeros of 

dx _ 1 r M " 6 
dr ~ 7(0 ~ r M " 6 - 2nr d ~ 3 + 2 
(i.e., r = only). We have the additional problem that £ is a multivalued function: each of the "horizons" 

r = R„, n = 0, 1, . . . ,d - 4 

(see appendix |C|) is a branch point. From ( |3.4| ), we see that although the function Re(x) is still well defined 
around with n = and n = (if d is odd), as f'(R^) is real in these cases, it will be multivalued 
around all the other fictitious horizons. 

For r ~ oo one has f(r) ~ 1 and hence x ~ r. Consequently x is holomorphic at infinity, and we can 
choose the branch cuts to cancel out among themselves. Therefore Re(x) is well defined in a neighborhood 
of infinity, and moreover Re(x) = will be approximately parallel to Re(r) = in this neighborhood. Two 
of the 2(2d — 5) branches of the Stokes line starting out at the origin must therefore be unbounded. The 
remaining can either connect to another branch or end up in a branch cut. On the other hand, it is easy 
to see that the Stokes line must intersect the positive real axis exactly in two points, one in each of the 
intervals (Rq ,Rq) and (Rq, +oo). Using this information plus elementary considerations of symmetry and 
the sign of Re(x), one can deduce that the Stokes line must be of the form indicated in Figure ||. These 
results are moreover verified by the numerical computation of the same Stokes lines, as indicated in Figure 
I 

We consider the contour obtained by closing the unlimited portions of the Stokes line near r ~ oo 
and avoiding the enclosure of the inner horizon, as shown in Figure [3|. At point A we have ujx > 0, and 
therefore the expansion fl3.13|) holds at this point. Imposing condition ( |3.12| ) one obtains 
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Figure 3: Stokes line for the Reissner-Nordstrom black hole, along with the chosen contour for monodromy matching, 
in the case d = 6. 




Figure 4: Numerical calculation of the Stokes lines for the Reissner-Nordstrom black hole in dimensions d = 4, 
d = 5, d = 6 and d = 7. The different shadings also illustrate the various horizon singularities and branch cuts. 



B + e~ ia+ + B-e~ ia - = 0. (3.14) 

For z ~ one has the expansion ([TI]). Consequently, as one rotates from the branch containing point A 
to the next branch in the contour (through an angle of 23^5)1 00 ~ (2d-b)o' 2 T °^ a ^ es through an angle of 2ir, 
and since 

V2ne 2m ujx J, L (e 2m ujx) = e^r {1±j) V2itljx J, l (us) ~ 2e 4ia± cos(ujx - a±) 

=t 2 ± 2 

one has 
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$(x) ~ 2B + e 4ia + cos (cjx - a+) + 25_e 4ia ~ cos {cox - a_) 

= ( J B + e 3iQ + + B„e 3ia -) e iujx + (B+e 5 ^ + B_e 5ia ~) e~ iuJX (3.15) 



on that branch. As one follows the contour around the inner horizon r = R , it is easily seen from ( p.4| 
that x approaches the point 

2iri iri 



f'(Ro) k 

where 

k- = \f'(R ) 

is the surface gravity at the inner horizon (notice that k~ < 0). Consequently, after going around the inner 
horizon, the solution will be of the form 



$(x) ~ C+J2iruj{x - 5) J L (u(x - 5)) + C-J2iru{x - 5) J_ L (co(x - 5)) , 

2 2 

as one approaches the origin. Since u(x — 6) is negative on this branch, from the asymptotic expansion 

Jv{z) = \—cos[z + — + -), z<-l, (3.16) 



TTZ \ 2 4 

we have 



<&(x) ~ 2C+ cos (co(x - 5) + a+) + 2C_ cos (u(x - 5) + a_ 



+ C-e ia -e- iujS ) e iuJX + ( C+e- ia +e iujS + CLe - **-^ 



e 



on this branch. This must be matched to the expression ( 3.15 ) for $>(x), and hence we must have 



B+e 3ia + + B.e Ma - = C+e ia +e~ iui5 + C-e ia -e~ lu)S (3.17) 

Finally, we must rotate to the branch containing point B. Again x — 5 rotates through an angle of 2ir, and 
since 



2m 2 ^uo{x - 5) J ±1 (e 2ni co(x - 5)) = e^ {1±j) ^2ttoo(x - 5) J ±L (u(x - 5)) ~ 2e 4ia± cos(w(x - 5) + a±) 
(as co(x — 5) < on the branch containing point B), one has 

$(x) ~ 2C + e 4ia + cos (u(x -5)+ a+) + 2C-e 4ia ~ cos (v(x - 5) + a_) 
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on that branch. As one closes the contour near r ~ oo, one has x ~ r and hence Me(x) > 0; since 
Im(uj) ^> 0, it follows that e luJX is exponentially small in this part of the contour, and therefore only the 
coefficient of e~ luJX should be trusted. As one follows the contour, this coefficient is multiplied by 



C 4 



+ C-e 



B + e ia + + B_e l 



TYLO 



On the other hand, the monodromy of e~ tu)X going clockwise around this contour is e~ 

We can now deform the contour to a small clockwise circle around the black hole event horizon. Near 
the outer horizon r = Rq , we have again V ~ 0, and hence again 



<f>(x) ~ D + e iuJX + D_e~ iujx . 

The condition for quasinormal modes at the horizon is therefore D- = 0. Using (|3.4| ), we restate this 
condition as the condition that the monodromy of going clockwise around the contour should be 



where 



/'(«+)/ _ eI + 



is the surface gravity at the outer horizon. Therefore, one finally concludes that the condition for quasi- 
normal modes at the horizon is 



C+e 3i a+e iu>5 + c _ e 3ia. e iu,5 

B + e ia + +B_e ia - ' 



(3.19) 



The final condition for equations (3.14), (3.17), (3.18) and (3.19), which reflect quasinormal mode 
boundary conditions both at infinity and at the black hole event horizon, to have nontrivial solutions 

(B + ,B_,C+,C_)is 



e -«o+ 
e 5ia + 



_^—ia.— giujS 



o <^ 



e 2ia+ 
e 4ia + 


e- 2ia - 

e 2ia. 




2 TV 

e k- 

e -2ia + 




e k- 

e -2ia 


e k+ 




e 2ia+ 


e 2ia. 


2nuj 

e k- 


= 







2-rru 

e k+ 



-(1 + 2cos(vrj)) - (2 + 2cos(7rj))e" 



2ttu 
k~ 



(3.20) 



As in the Schwarzschild case, the fact that we imposed condition (|3.12| ) at point point A means that 
equation ( |3.20| ) yields the asymptotic quasinormal modes to with Re(i<j) > 0. The remaining modes are of 
the form — uJ, and could be obtained by imposing condition ( p. 12 ) at point point B. 
For vector type perturbations, the potential near the origin is of the form 



V[r{x)] 



J' 



Ax 2 ' 
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with j = Repeating the same argument, one ends up with the same equation ( |3,20| ) for this new 

value of j. Since 

3d -7 _ d - 3 



2d - 5 2d - 5 

we see that this equation is exactly the same as in the j = case, and consequently we end up with 
the same quasinormal modes. 

One may wonder if it is possible to recover the Schwarzschild quasinormal frequencies from our result 
above. This would arise from the # — > limit of (|3.20| ). In this limit we have Rq — > and hence 

(d-3)/x 

k ~ -7-= — > -oo, 

while 



K + ) d " 2 

remains fixed. In this case it is simple to check that the RN quasinormal equation reduces to 



ek+ = - (3 + 4 cos (nj)) , 

which is not the quasinormal equation for the Schwarzschild spacetime. For example, in d = 4, this formula 
yields e >=+ = —5. The reason why this is so traces back to the Stokes line. Comparing the Stokes lines of 
the Schwarzschild and RN spacetimes one immediately realizes that there is a change in the topology of this 
curve at the origin, thus invalidating the limit above. Another limit one could try to obtain concerns the 
quasinormal frequencies of the extremal RN black hole, when — > /z. This time around the computation 
is slightly trickier and in the end one will obtain the equation in the extremal limit as 



2ttuj 



e k = - (2 + 2 cos (7rj)) , 
where we have defined 



2{d - 2)R 



i 



with Rq = /x 53 ^ . As we shall see in the following this equation yields the wrong quasinormal frequencies for 
the extremal RN spacetime (although it actually predicts the correct gap). In this case there is no change in 
the topology of the Stokes line, but it turns out that one cannot deform the monodromy contour of the RN 
black hole into the monodromy contour of the extremal RN black hole without crossing horizon singularities 
in the complex plane, thus invalidating the limit above. Interestingly enough, the same parameter k as 
above will also appear when performing the extremal RN calculation from scratch. The lesson one should 
retain is that in many cases one must be careful when taking limits of the final equation for the asymptotic 
quasinormal frequencies. 

Let us briefly review the literature concerning asymptotic quasinormal frequencies in the RN space- 
time, as we would like to compare our results to what has been previously accomplished on this subject. 
Quasinormal frequencies for near extremal RN solutions were first studied numerically in p8[ , focusing on 
the four dimensional case. The highly damped modes were found to have a peculiar behavior, as black 
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holes with different charge apparently shared a specific mode frequency, and some of the modes seemed 
to be multivalued. The first analytic solution, in four dimensions and for non-extremal RN solutions, was 
given in [jl4| |l7fl . These authors found a formula (which we have generalized above) for the asymptotic 
quasinormal frequencies, but they have not managed to find its general solution. Also, simple numeric 
analysis seemed to point towards the peculiarities previously found. This is in full agreement with what 
we have found, and just seems to be pointing in the direction that there is no simple solution for the 
frequencies in this case, unlike the previous Schwarzschild solution. Numerical checks to the formula in 



|14|| were performed in |3Sj, [39| again found that very highly damped quasinormal modes of RN black 
holes have an oscillatory behavior as a function of the charge. [^] found that RN asymptotic quasinormal 
frequencies are typically not periodic in lm(uj) (in contrast with the Schwarzschild black hole case). They 
also understood why the 9^0 limit of RN quasinormal frequencies does not yield Schwarzschild quasi- 
normal frequencies: the limit is a singular one as it involves topology change at the level of the contours 
in the complex plane (the same thing happening in the attempt to go extremal, 2 —* fi 2 , as we shall 
see in detail in the following). In summary, our results completely match what was previously known in 
the literature. However, it would be very interesting to produce further numerical data concerning higher 
dimensional RN quasinormal modes to match against all our analytical results. 

3.1.3 The Extremal Reissner Nordstrom Solution 

We now compute the asymptotic quasinormal modes of the d-dimensional extremal RN black hole using 
the monodromy method. This example is illustrative of the fact that the monodromy calculation is very 
sensitive to the location of the singularities in the complex plane, as well as to the topology of the Stokes 
line at the origin. Thus, when taking limits of the parameters, one should always have in mind whether 
one is crossing singularities or changing the topology of the contour. In an affirmative case, the limit on 
the parameters will not be valid and one will have to address the calculation from scratch. 



Again we consider solutions of the Schrodinger-like equation (3J2) in the complex r-plane. Since the 
potential V vanishes for r ~ +oo, we will have 

in this region. The boundary condition for quasinormal modes at infinity is then 

A + = 0. (3.21) 

Next we study the behavior of $(x) near the singularity r = 0. In this region, the tortoise coordinate is 
still 



r 2d-b 



X ~ 



(2d-5)6> 2 ' 

and the potential for scalar type and tensor type perturbations is still 

with j = ^3 ( see appendix |C|). Therefore the solution of the Schrodinger-like equation in this region 
approximates 

J l (ra) +B_V2 1TLOX J j (LUX) , 
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where J v {x) represents a Bessel function of the first kind and B± are (complex) integration constants. 

For the asymptotic quasinormal modes one has Im(uj) 3> Me (a;), and hence u is approximately purely 
imaginary. Consequently, one has ojx G IR for x G iM,; in a neighborhood of the origin, this happens for 



r — p g 2(2d-5) ~ 2d-5 

with p > and n = 0, 1, . . . , Ad — 11. These are half-lines starting at the origin, equally spaced by an 
angle of 2^5- Notice that the sign of ojx on these lines is (— l) n+ ; in other words, starting with the 
line corresponding to n = 0, the sign of ojx is alternately negative and positive as one goes anti-clockwise 
around the origin. 

From the asymptotic expansion ( |3.6[ ) we see that 



<&(x) ~ 2B + cos (ojx — a + ) + 2B_ cos (ojx — q_) 

= (B + e- ia + + B„e~ ia -) e iuJX + (B + e ia + + B_e ta -) e^ x , (3.22) 

in any one of the lines corresponding to positive ojx, where again 

a ± = j(l±j). 

This asymptotic expression for <3? near the origin is to be matched with its asymptotic expression at infinity. 
This matching must again be done along the Stokes line ojx G M. (that is, x G iR), so that neither of the 
exponentials e illu)X dominates the other. 

To trace out the Stokes line Re(x) = we observe that we already know its behavior near the origin. 
Furthermore, this is the only singular point of this curve: indeed, since x is a holomorphic function of r, 
the critical points of the function Re (2;) are the zeros of 

dx 1 r M ~ 6 
dr ~ W) ~ ( r d-3 _ 0)2 

(i.e., r = only). We have the additional problem that x is a multivalued function: each of the "horizons" 

r = R n , n = 0, 1, . . . , d — 4 

(see appendix [C]) is a branch point. From ( |3.4| ) we see that although the function Me(x) is still well defined 
around Rk with n = and n = (if d is odd), it will be multivalued around all the other fictitious 
horizons. In analogy with the non-extremal case, we are led to define 

Jd-3)^ 
2(d - 2)R Q 

For r ~ 00 one has f(r) ~ 1 and hence x ~ r. Consequently x is holomorphic at infinity, and we can 
choose the branch cuts to cancel out among themselves. Therefore Me(a;) is well defined in a neighborhood 
of infinity, and moreover Me(x) = will be approximately parallel to Me(r) = in this neighborhood. Two 
of the 2(2d — 5) branches of the Stokes line starting out at the origin must therefore be unbounded. The 
remaining can either connect to another branch or end up in a branch cut. On the other hand, it is easy to 
see that the Stokes line must intersect the positive real axis exactly in two points, the horizon r = Rq and 
a point on the interval (Rq, +00). Using this information plus elementary considerations on symmetry and 
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contour 




Figure 5: Stokes line for the extremal Reissner-Nordstrom black hole, along with the chosen contour for monodromy 
matching, in the case d = 6. 







Figure 6: Numerical calculation of the Stokes lines for the extremal Reissner-Nordstrom black hole in dimensions 
d = A, d = 5, d = 6 and d = 7. Different shadings also illustrate the various horizon singularities and branch cuts. 

the sign of Me (a:), one can deduce that the Stokes line must be of the form indicated in Figure [|. These 
results are moreover verified by the numerical computation of the same Stokes lines, as indicates in Figure 



We consider the contour obtained by closing the unlimited portions of the Stokes line near r ~ oo, as 
shown in Figure ||. At point A we have ujx > 0, and therefore the expansion ( 3.22j ) holds at this point. 



Imposing condition (3.21) one obtains 



B +e - ia + + B^e~ ia - = 0. (3.23) 

For z ~ one has the expansion ( p.9| ). Consequently, as on rotates from the branch containing point A to 
the branch containing point B (through an angle of 2 ^ 5 ), x ~ r 2d _^ rotates through an angle of 5ir, and 
since 
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one has (notice that e 5m ujx 



± i (e 5m to.r 



LUX 



2e 



10ia± 



cos(lox — a±) 



-UJX) 



~ 2B + e l0ia + cos {-ux - a+) + 2B^e 10ia ' cos {-ux - a_) 
= (B + e 1Ua + + B_e 1Ua ~) e iuJX + (B + e 9ia + + B_e 9ia -) eT wx 

at point B. As one closes the contour near r ~ oo, one has x ~ r and hence Me(x) > 0; since Im(ct>) 3> 0, 
it follows that e tuJX is exponentially small in this part of the contour, and therefore only the coefficient of 
e -iujx ghouls De trusted. As one follows the contour, this coefficient is multiplied by 

B+e 9ia+ + B_ e 9ia - 
B + e ia + +B_e ia - 

On the other hand, the monodromy of e~ lu)X going clockwise around this contour is e~~ . We can now 
deform the contour to a small clockwise circle around the black hole event horizon. Near the horizon 
r = Rq, we have again V ~ 0, and hence again 



The condition for quasinormal modes at the horizon is therefore C_ = 0. We restate this condition as 
the condition that the monodromy of going clockwise around the contour should be e~ . Therefore the 
condition for quasinormal modes at the horizon is 



B+e 9ia+ +B e 

B + e ia + + B_e 



9ia_ 



e fc . 



(3.24) 



The final condition for equations ( |3.23| ), ( |3.24 ), reflecting quasinormal boundary conditions, to have non- 
trivial solutions (B + , S_ ) is 



2tvuj 

e +e k 



27TUJ 

e ia -e— 



4^ 



2^ Sil1 

e k = — 



8111 



(3.25) 



As before, the fact that we imposed condition (|3.2l[) at point point A means that equation ( 3.25| ) yields 
the asymptotic quasinormal modes uj with Me(w) > 0. The remaining modes are of the form —u>, and 
could have been obtained by imposing condition (|3.2l| ) at point B. This result also illustrates the fact, 
which we have advertised earlier, that the extremal limit of the RN quasinormal equation does not yield 
the quasinormal equation for the extremal RN spacetime. 

For vector type perturbations, the potential for the Schrodinger-like equation near the origin is of the 
form 



V[r(x)] 



J' 



Ax 2 



with j 



3d- 7 



2rf _ 5 . Repeating the same argument, one ends up with the same equation ( |3,25| ) for this new 
value of j. Since 
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3d -7 _ 2 d-3 
2d -5 ~~ ~ 2d - 5' 

we see that this equation is exactly the same as in the j = ^3 case, and consequently we end up with 
the same the quasinormal modes. 

As one reviews the literature searching for earlier numerical checks on our results, one observes that 
there has not been much work on this matter. Quasinormal frequencies for extremal RN solutions were 
studied numerically in [fjjj ], focusing on the four dimensional case. However, the interest of that paper 
lied in the least damped modes (opposite to what we want). The authors found a coincidence in quasi- 
normal frequencies for gravitational perturbations with multipole index I and those for electromagnetic 
perturbations with index t — 1 (at extremality) . This fact was later explained in [^] as arising from the 
[hidden] supersymmetry of the extremal solution. So, in this extremal limit, the black hole responds to 
the perturbation of each field in the same manner. This is certainly in agreement with our results, but 
as we have seen for tensor, vector and scalar type perturbations this is actually a very generic feature. 
An attempt at computing asymptotic quasinormal frequencies for the extremal RN black hole, using the 



monodromy method of [14], was recently performed in [43|. Unfortunately, these authors have misidentified 
the topology of the tortoise coordinate at the origin and have thus wrongly computed the monodromies for 
their black holes, obtaining an incorrect result for the quasinormal frequencies. An attempt at computing 
numerically quasinormal frequencies of the d = 4 extremal RN black hole was recently done in , but un- 
fortunately the numerical method was not stable enough to find the asymptotic values for the frequencies. 
It would thus be very interesting to produce further numerical data concerning extremal RN quasinormal 
modes in order to match against our analytical results. 

3.2 Asymptotically de Sitter Spacetimes 



|10| discusses the stability of black holes in asymptotically dS spacetimes to tensor, vector and scalar 
perturbations. For black holes without charge, tensor and vector perturbations are stable in any dimension. 
Scalar perturbations are stable up to dimension six but there is no proof of stability in dimension d > 7. 
For charged black holes, tensor and vector perturbations are stable in any dimension. Scalar perturbations 
are stable in four and five dimensions but there is no proof of stability in dimension d > 6. As we work in 
generic dimension d we are thus not guaranteed to always have a stable solution. Our results will apply if 
and only if the spacetime in consideration is stable. 



Quantization of a scalar field in dS space was first addressed in |45 ] . While these authors found that 
the cosmological event horizon is stable, they also found that there is an isotropic background of thermal 
radiation. The emitted particles are, however, observer dependent, as is the "cosmological sphere" of dS 
(we refer the reader to the recent review [46|, for more details and a list of current open problems). Analysis 
of the wave equation in dS spacetime also led to the natural boundary conditions on quasinormal modes: 
incoming waves at the black hole horizon and outgoing waves at the cosmological horizon. 

3.2.1 The Schwarzschild de Sitter Solution 

We now compute the quasinormal modes of the Schwarzschild dS (f-dimensional black hole. This calculation 
was first done in for the particular case of d = 4. Again we consider solutions of the Schrodinger-like 
equation ( p.2| ) in the complex r-plane. We start by studying the behavior of $(x) near the singularity 
r = 0. In this region, the tortoise coordinate is 
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r d-2 

X ~~2(d-2)^ 

and the potential for tensor and scalar type perturbations is 

with j = (see appendix |C|). The solution of the Schrodinger-like equation in this region (with j ^ 0) is 
therefore 

ttujx J i (ujx) + B- V 2ixuix J_i (ujx) , 

2 2 

where J v (x) represents a Bessel function of the first kind and B± are (complex) integration constants. 

For the asymptotic quasinormal modes one has Im(u;) 3> Re(a>), with Im(cj) — > +oo, and hence cj is 
approximately purely imaginary. Consequently, ujx £ R for x G zR; in a neighborhood of the origin, the 
above relation between x and r tells us that this happens for 



y — p g 2(d— 2) T d-2 

with p > and n = 0, 1, . . . , 2d — 5. These are half-lines starting at the origin, equally spaced by an angle 
of ^2 . Notice that the sign of lux on these half-lines is (— l) n ; in other words, starting with the half-line 
corresponding to n = 0, the sign of ujx is alternately positive and negative as one goes anti-clockwise 
around the origin. 

From the asymptotic expansion ( |3.6[) we see that 



~ 2B + cos (ujx — a+) + 2B- cos (ujx — q_) 

= (B+e~ ia + + B-e~ ia -) e iuJX + (£+e ia + + £_e ia ") e"^ (3.26) 
in any one of the lines corresponding to positive ujx, where again 

a± = ^(l±j). 

We shall use this asymptotic expression to make the monodromy matching. This matching must be done 
along the Stokes line ujx £ R (or x G iM), so that neither of the exponentials e illu)X dominates the other. 

To trace out the Stokes line Re(x) = let us first observe that we already know its behavior near the 
origin. Furthermore, this is the only singular point of this curve: indeed, since x is a holomorphic function 
of r, the critical points of the function Re(x) are the zeros of 



dx 1 r 



d-3 



dr f(r) -Ar^- 1 + r d - 3 - 2/i 
(i.e., r = only). We have an additional problem that x is a multivalued function: each of the "horizons" 



R n (n = 1, . . . ,d — 1) is a branch point. From (3.4) we see that although the function Re(x) is still well 
defined around the real horizons, it will be multivalued around all the other fictitious horizons. 

For r ~ oo one has x ~ xq + ^ (see appendix |C]). Thus, x is holomorphic at infinity, and we can 
choose the branch cuts to cancel out among themselves. Therefore Re(x) is well defined in a neighborhood 
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Figure 7: Stokes line for the Schwarzschild de Sitter black hole, along with the chosen contour for monodromy 
matching, in the case of dimension d = 6. 




Figure 8: Numerical calculation of the Stokes lines for the Schwarzschild de Sitter black hole in dimensions d = 4, 
d — 5, d = 6 and d = 7. The different shadings also illustrate the various horizon singularities and branch cuts (note 
that these branch cuts are not necessarily equal to the ones used for the calculation in the main text). 

of infinity, and moreover Re(ar) = cannot extend out to infinity as, generically, xq is not real. Thus the 
Stokes line branches starting at the origin must either connect to another branch or end up in a branch cut 
(and one would expect them to connect in the region where Re(x) is well defined). On the other hand, it 
is easy to see that the Stokes line must intersect the positive real axis exactly in two points: one between 
the black hole horizon Rh and the cosmological horizon Rc, the other between the cosmological horizon 
and infinity. Using this information plus elementary considerations on symmetry and the sign of Re(x), 
one can deduce that the Stokes line must be of the form indicated in Figure [?]. These results are moreover 
verified by the numerical computation of the same Stokes lines, as indicated in Figure |8[ 

Let us now consider the contour defined by the closed part of the Stokes line (contained in the region 
where Re(x) is well defined), as shown in Figure]^. At point A we have ujx 3> 1, and therefore the expansion 
( 3.26| ) holds at this point. For z ~ one has the expansion (^fi|). Consequently, as one rotates from the 
branch containing point A to the branch containing point B, through an angle of , x ~ ~ 2(d-2)(i ro ^ a ^ es 
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through an angle of 3ir, and since 

^2ne 3ni ujx J ±l (e 3nl u;x) = e^ {1±j] V2ttwx J ±l (cox) ~ 2e 6ia± cos(wx - a±) 



one has at point B (notice that e ujx 



-ojx) 



$(x) ~ 2B + e eia+ cos {-ujx - a+) + 2B„e 6ia - cos (-ws - a_) 

= (B + e 7ia + + B-e 7ia ~) e iujx + (S+e 5ia + + £_e 5fo ") e"^. (3.27) 

If $ corresponds to a quasinormal mode, its clockwise monodromy around r = Rjj must be the same 
as the clockwise monodromy of e luJX , that is, e k H , where 

*H = \f{RH) 

is the surface gravity at the black hole horizon. Similarly, its clockwise monodromy around r = Rc must 
be the same as the clockwise monodromy of e~ UJX , that is, e k c , where 

kc = \f{Rc) 

is the surface gravity at the cosmological horizon (notice that we are taking kc < 0). Since these are the 
only singularities of $ inside the contour, the monodromy of <I> around the contour must be 

7TOJ 7TU 

e k H k c . 

7VUJ j 7TUJ 

On the other hand, the clockwise monodromy of e luJX around the contour is e k H k c . Moreover, we have 
just learned from ( p.26| ) and ( p. 27 ) that as one goes around the contour the coefficient of e tuJX in the 
asymptotic expansion of $ gets multiplied by 

B + e 7ia + + B_e 7ia - 
B +e - ia + + B-e~ ia - ' 

For this term to have the required monodromy we must impose 



B+e 1%a + +B_e 



TVUJ 

e k H^ k c =e k a K c 



B+e 



7ictj 



+ B_e 



Hot- 



B + e~ ia + + B_e~ ia 
Similarly, for the term in e~ wx we get the condition 

B+e Ua + +B_e 5ia 



TTUJ _7Ti 

k H k, 



g H K C = e H K C 44> 



TYUJ TTUJ 

k H k r 



B + e- ia + + B_e~ 



B e 5ia+ +B _ e 5ia. 



B + e ia + + B_e ia - B + e ia + + B_e ta - 

The condition for these equations to have nontrivial solutions {B + , i?_ ) is then 



2ttuj 

e k H . 



K c e 



e K H e 



e k c e 



e K n e 



O 



sin ( 4a 



l TYUJ 

kc 



sin ( 2a + + f± 



sin ( 4a_ - ^ 



sin(2a_ + ^ 



0. 
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As in the Schwarzschild case, this equation is automatically satisfied for j = 0. This is to be expected, 
as for j = the Bessel functions J, i coincide and do not form a basis for the space of solutions of the 

=t 2 

Schrodinger-like equation near the origin. As in [14|, we consider this equation for j nonzero and take 
the limit as j — > 0. This amounts to writing the equation as a power series in j and equating to zero the 
first nonvanishing coefficient, which in this case is the coefficient of the linear part. Thus, we just have to 
require that the derivative of the determinant above with respect to j be zero for j = 0. This amounts to 



7T COS TT 



kc 



-TT COS TT 



kc 



sin(f + ^ sin f + ^ 



+ 



sin TT 



kc 



sin TT 



kc 



0. 



from where we obtain our final result as 



/ TTU) TTU>\ ( ITU! TTU)\ 

cosh — +3 cosh h - — 

\k H kc ) \k H k c ) 



0. 



(3.28) 



Notice that if a; is a solution of this equation then so is —u), as must be the case with quasinormal modes. 

From the expression above it is possible to recover the Schwarzschild quasinormal frequencies. We first 
write the equation as 



1YUJ I 7T(jJ 



7TCJ 7TLtJ 



e k H k c + e k H k c -\-'5e k H k c + 3e k H 



0. 



and take the limit as A — > + . In this limit, we have Rc ~ A 2, and hence Rc — > +oo; therefore 



k c 



- XRc = ~Rc- 



0~ 



If we assume that Re(o>) > 0, we see that the two middle terms are exponentially small, and hence the 
quasinormal equation reduces to 



e k n + 3e k H = 44> e k H 



which is exactly the equation obtained in [13, 14]. Therefore, the Schwarzschild black hole is not a singular 
limit of the Schwarzschild dS black hole as far as quasinormal modes are concerned, unlike what happens 
with the RN black hole solution. The reason for this is clear from the monodromy calculation: whereas 
the structure of the tortoise near the singularity r = in the RN solution depends crucially on whether 
the charge is zero or not, in the Schwarzschild dS case it does not depend on A. Thus, as R — > +oo, the 
cosmological horizon approaches the point at infinity and the contour approaches the contour first used in 
[14]. One could also try to derive the /x — ► limit, trying to obtain quasinormal frequencies for pure dS 
spacetime. In this case Rh — > and 



k H 



d-3 



Rh -» 0+ 



while kc — > — vA. The quasinormal equation for the Schwarzschild dS spacetime becomes 



e k c 
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from where one obtains 



V 2 2vri J 



Unfortunately, as we shall later see in section 4, this is the incorrect value for the dS quasinormal frequencies. 
While there was no expectation that the offset would come out right, one could have hoped that the gap 
could have been correct. We shall later learn that the correct gap in pure dS spacetime is 2^f\ and not 
\/A as obtained above. Thus, while the Schwarzschild spacetime is a good limit of the Schwarzschild dS 
black hole, pure dS spacetime is not, at least as far as asymptotic quasinormal frequencies are concerned. 

For vector type perturbations, the potential for the Schrodinger-like master equation near the origin 
is of the form 

with j = 2 (see appendix |C]). Repeating the same argument, one ends up with the same equation, except 
that now 2a± = ^ ± tt, 4a± = ir ± 2tt. This equation is therefore exactly the same as in the j = case, 
for which 2a± = ^, 4a± = tt, and consequently we end up with the precise same quasinormal frequencies. 

Some remarks are due for d = 4 and d = 5, as in these cases the Stokes lines are not quite of the 
general form we assumed in the previous calculation. We start by studying the behavior of solutions to 
the Schrodinger-like equation ([T^) near infinity. For r ~ oo one has x ~ xq for a (complex) constant xq, 
and the potential is 

? 2 - 1 



4(x — xo) 2 

with = d — 1, d — 3, d — 5 for tensor type, vector type and scalar type perturbations (see appendix |C]). 
Hence, 



<I>(x) ~ A + y / 27TUj(x - X ) Jjoo (U)(x - X )) + A_^2lT0j(x - Xq) J_ioo (u(x - Xq)) 



2 2 



in this region, where J u (x) represents a Bessel function of the first kind and A± are (complex) integration 
constants. Given expansion Q3.9| ), we see that $ has monodromy (— l) d around infinity. The d = 4 case 
was actually done in |^| . In this case, one can take advantage of the fact that $ has zero monodromy at 
infinity to conclude that as one follows the Stokes line the (anti-clockwise) monodromy of $ around the 
fictitious horizon r = R equals the (clockwise) monodromy around the real horizons r = Rh and r = Rq', 
since the same is true for the plane wave functions e ±lu)X , the whole calculation reduces to what was done 
above, and we get the same formula. Another way to see this is that since x, $ are regular at infinity one 
can deform the d = 4 contour into, say, the d = 6 contour in the Riemann sphere C U {oo}. 

For d = 5 the situation is more subtle. As shown in Figure ||, the Stokes line approaches infinity, and 
hence one must close the contour near infinity. In the branch containing point B, the solution is then of 
the form 

<£(x) ~ A + \/2mu) (x - xq) Jjoo (u>(x - x )) + A_ y / 27rw(x - x ) J_joo (uj(x - x )) 



e -iuj(x-x ) 
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where [5± = ^(1 ± joo), with joo =4,2,0 for tensor type, vector type and scalar type perturbations. As r 
rotates clockwise towards the branch containing point A, by ^, x — xq ~ rotates anti-clockwise by the 
same amount, and hence 

$(a;) ~ (A + e 3i(3 + + A_e 3i ^ e ^-*o) + (^A + e if3 + + A_e i/3 ~) e -M*-*o) 

at point B. Whereas the coefficient of e~ luJX is unchanged in this process, it is easy to see that the 
coefficient of e lulx reverses sign for all three values of joo. Therefore, there appears an extra minus sign 
in the corresponding monodromy matching equation, and the condition for nontrivial solutions (B+,B-) 
changes slightly. Carrying out the computation, the condition for quasinormal modes in d = 5 turns out 
to be 



. , f-KUJ 7TuA n . , / TTU 7TuA 

sinn — — 6 sinn h - — = U. 

\k H k c ) \k H k c ) 

It is easily seen that this equation still has the correct Schwarzschild limit. 

Let us now review the literature concerning asymptotic quasinormal frequencies in the Schwarzschild 
dS spacetime, so that we can compare our results to what has been previously accomplished on this 



subject. For Schwarzschild dS, early results on quasinormal modes were studied in [47], without great 
emphasis on the asymptotic case. The first analytical results in d = 4 were derived in the near-extremal 
situation, where event and cosmological horizons are nearly coincident [48], but the approximation used 



therein is not expected to hold in the asymptotic limit, at least on what concerns the real part of the 
asymptotic frequencies. These results were later extended to ^-dimensions (and to RN dS black holes) 
in p9| , with similar results. Further approximations were studied in J50[ |, in a limit where the black hole 
mass is much smaller than the spacetime radius of curvature, but focusing explicitly on the time dependent 
transient situation. An attempt at an analytic solution for the asymptotic quasinormal frequencies, using 



the monodromy technique of Ej], was done in 51]. However, an erroneous identification of the relevant 
contours led these authors to a wrong result — in fact, the authors assumed the monodromy contour to be 
the same as in the Schwarzschild case, while as we have seen the contour changes dramatically in non- 
asymptotically flat spacetimes (see also p2l , where other arguments were given trying to explain the failure 



of [51] to reproduce available numerical data). Focusing on d = 4, there are some interesting analytical 
results in ]3^| . These authors find that because there are two different surface gravities, there are also 
two sets of solutions for Im (uj) when the horizons are widely spaced, namely Im (cj) equally spaced with 
spacing equal to ku or Im (u>) equally spaced with spacing equal to kc- It was further claimed in |}^] that 
this lack of consensus on quasinormal frequencies was due to the fact that there is no global definition 
of temperature in this spacetime. Our results appear to confirm this expectation: in the limit where the 
cosmological radius goes to infinity, and one recovers the Schwarzschild modes, we found the spacing to be 
equal to kn- In the limit where the black hole radius is very small, our final formula yields modes with 
spacing equal to kc- 



The four dimensional case was first solved analytically in [16]. As discussed in that paper, the analytical 



results are in complete agreement, and to a large degree of accuracy, with the numerical results provided 



(for this case of d = 4 Schwarzschild dS) in 52, 54 1. Alongside with the fact that our expressions yield 



the correct Schwarzschild limit in any dimension, we believe that there is strong evidence supporting our 
general results. However, it would be very interesting to produce further numerical data concerning higher 
dimensional Schwarzschild dS quasinormal modes to match against all our analytical results. 
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3.2.2 The Reissner Nordstrom de Sitter Solution 



We now compute the quasinormal modes of the RN dS d-dimensional black hole. Again we consider 
solutions of the Schrodinger-like equation ( |3.2| ) in the complex r-plane. We start by studying the behavior 
of <E>(x) near the singularity r = 0. In this region, the tortoise coordinate is 



r 2d-5 



(2d-5)0 2 ' 

and the potential for tensor and scalar type perturbations is 



V[r(x)] 



Ax 2 



with j = T^pl (see appendix |C|). The solution of the Schrodinger-like equation in this region (with j / 0) 
is therefore 



$(x) ~ B+V2 TTUJX TTUJX J 

2 2 

where J v {x) represents a Bessel function of the first kind and B± are (complex) integration constants. 

For the asymptotic quasinormal modes one has Im(u;) 3> Re(u;), with Im(cj) — > +oo, and hence u is 
approximately purely imaginary. Consequently, lux G R for x G iM; in a neighborhood of the origin, the 
above relation between x and r tells us that this happens for 



r — p g 2(2d-5) ~ r 2ii-5 

with p > and n = 0, 1, . . . , Ad — 11. These are half-lines starting at the origin, equally spaced by an 
angle of 2 J_ 5 . Notice that the sign of ujx on these lines is (— l) n+1 ; in other words, starting with the 
line corresponding to n = 0, the sign of ujx is alternately negative and positive as one goes anti-clockwise 
around the origin. 

From the asymptotic expansion ( p.6[ ) we see that 



~ 2B + cos {ujx — a + ) + 2B_ cos (ujx — a_) 

= (B + e~ ia + + B_e~ ia -) e iuJX + (B + e ia + + B_e ia ~) e~ iu>x (3.29) 
in any one of the lines corresponding to positive ujx, where again 

We shall use this asymptotic expression to make the monodromy matching. This matching must be done 
along the Stokes line ujx £ R (or x E iR), so that neither of the exponentials e ±lu3X dominates the other. 

To trace out the Stokes line Re(x) = let us first observe that we already know its behavior near the 
origin. Furthermore, this is the only singular point of this curve: indeed, since x is a holomorphic function 
of r, the critical points of the function Re(x) are the zeros of 

dx _ 1 r M - 6 

dr ~ f(r) ~ -Ar M " 4 + r 2d ~ 6 - 2/ir d ~ 3 + 2 
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Figure 9: Stokes line for the Reissner-Nordstrom de Sitter black hole, along with the chosen contour for monodromy 
matching, in the case of dimension d = 6. 




Figure 10: Numerical calculation of the Stokes lines for the Reissner-Nordstrom de Sitter black hole in dimensions 
d = 4, d = 5, d ~ 6 and d = 7. Different shadings also illustrate the various horizon singularities and branch cuts 
(note that these branch cuts are not necessarily equal to the ones used for the calculation in the main text). 



(i.e., r = only). We have an additional problem that x is a multivalued function: each of the "horizons" 
Rn {n = 1, . . . , 2d — 4) is a branch point. From ( |3.4| ) we see that although the function Me(:z;) is still well 
defined around the real horizons, it will be multivalued around all the other fictitious horizons. 

For r ~ oo one has x ~ xq + ^- (see appendix |C|). Consequently, x is holomorphic at infinity, and 
we can choose the branch cuts to cancel out among themselves. Therefore Me(x) is well defined in a 
neighborhood of infinity, and moreover Me(x) = cannot extend out to infinity, as generically xq is not 
real. Thus the Stokes line branches starting out at the origin must either connect to another branch or 
end up in a branch cut (and one would expect them to connect in the region where Me(x) is well defined). 
On the other hand, it is easy to see that the Stokes line must intersect the positive real axis exactly in 
three points, one in each of the intervals (Rjj, Rjj), (R~hi R c) and (Rc, +oo). Using this information plus 
elementary considerations on symmetry and the sign of Me(x), one can deduce that the Stokes line must 
be of the form indicated in Figure |^. These results are moreover verified by the numerical computation of 



the same Stokes lines, as indicated in Figures nil and 11 
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Figure 11: Close-up on the numerical calculation of the Stokes lines for the Reissner-Nordstrom de Sitter black 
hole in d = 4, d = 6 and d — 7. Locally near the origin these lines resemble the corresponding ones for the 
Reissner-Nordstrom black hole in an asymptotically flat spacetime (see Figure 0). 



Let us now consider the contour defined by the closed portions of the Stokes line in such a way that 
one avoids enclosing the inner horizon, as shown in Figure |9| (this contour is contained in the region where 
Me(x) is well defined). At point A we have ujx > 0, and therefore the expansion ( 3.2S| ) holds at this point. 
For z ~ one has the expansion ([^). Consequently, as one rotates from the branch containing point A 
to the next branch in the contour (through an angle of 2 rf-s )i x ~ ^xi?^ rotates through an angle of 2ir, 
and since 



(2d-5)6' 



one has 



v / 27re 27ri tjx J, 3 



e ¥( 1± JV2vra;x J, j 



UJX 



2g4ia± 



COS (UJX 



a±) 



$(x) ~ 2B + e 4ia + cos (ujx - a+) + 2B_e 4ia ~ cos {ujx - a_) 
= (S+e 3 ^ + B^e 3ia ~) e iwx + (S+e 5 ^ + B_e 5ia ~) e 



(3.30) 



on that branch. As one follows the contour around the inner horizon r = R H , it is easily seen from ( |3.4| 
that x approaches the point 



2vri 



f'(R- H ) 



where 



k~ = \f'(R H ) 

is the surface gravity at the inner horizon (notice that k~ < 0). Consequently, after going around the inner 
horizon, the solution will be of the form 



$(x) ~ C + J2ituj(x - 8) J L (u(x - 5)) + C-J2iru}(x-5) J_ L (u>(x - 5)) , 

2 2 

as one approaches the origin. Since uj(x — 5) is negative on this branch, from the asymptotic expansion 
(|3.16|) we have 
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$(x) ~ 2C+ cos (u(x - 5) + a+) + 2C_ cos - 5) + a. 



_e ~e e 



on this branch. This must be matched to the expression ( 3.30 ) for $>(x), and hence we must have 



B+e 3ia+ + B^e 3ia ~ = C + e ia +e~ iul& + C-e ia 'e' iu)6 (3.31) 
B+e 5ia+ + B_e 5ia ~ = C + e~ ia +e iul& + C-e^-e^ (3.32) 

Finally, we must rotate to the branch containing point B. Again x — 5 rotates through an angle of 2n, and 
since 



2m 2 ^uj{x - 5) J ± i (e 27Ti uj(x - 5)) = e^ {l±j ^2kuj(x - 5) J ±L (oj(x - 5)) ~ 2e 4ia± cos(w(x - 5) + a±) 
(as oj(x — 5) < on the branch containing point B), one has 

$(x) ~ 2C+e 4ta+ cos (u(x -5)+ a+) + 2C-e 4ia ~ cos (lv(x - 6) + a_) 



on that branch. 

If <I> corresponds to a quasinormal mode, its clockwise monodromy around r = must be the same 
as the clockwise monodromy of e iwx , that is, e*^, where 

k + = \f'(R + H ) 

is the surface gravity at the black hole horizon. Similarly, its clockwise monodromy around r = Rc must 
be the same as the clockwise monodromy of e~ tuJX , that is, e k c , where 

k c = \f{Rc) 

is the surface gravity at the cosmological horizon (notice that we are taking kc < 0). Since these are the 
only singularities of <£ inside the contour, the monodromy of <I> around the contour must be 



7TUJ 7TU. 



7YUJ _i_ 1TUJ 



On the other hand, the clockwise monodromy of e lu)X around the contour is e*^ k c . Moreover, we have 
just learned that as one goes around the contour the coefficient of e luJX in the asymptotic expansion of 
gets multiplied by 



B + e~ ia + +B_e~ ia - 
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For this term to have the required monodromy we must impose 



£1 gSiCK-)- r — /~l'-< 



iujd 



g feT k c 



B +e - ia + + B_e- %a - 
Similarly, for the term in e~ lulx we get the condition 



44> 



icuS _j_ £j gSia.— icuS 



B,e- ia + + B_e~ 



e k c . (3.33) 



B+e ia + + B_e ia - 



7TUJ 7TLU 



■KUS TVOJ 



Q +e 3ia +e iu5 + C _ e 3ia- e iu8 

B+e ia + +B_e ia - 



2-Tra 



(3.34) 



The condition for equations ( 3.31] ), fl3.32|) , ( |3.33| ), (|3.34| ) , reflecting quasinormal boundary conditions both 
at cosmological and outer black hole horizons, to have nontrivial solutions (B+, B-, C+, C_) is then 



2TTU1 

e ia +e~k~+ 



s 6ia + 
I k c 



4ia_ 



e fc" 1 



2ia_ 



e fc 4 



_g2Q-|_ g — 20^(5 



-e 



27ra 



_g«a_ g — itod 



_^3ia— ^icuS 



1 



e 2*a_g — 



e"""+e k 



gw— g fc - 



from where we obtain our final result as 



, / 1TU> 1TLU \ , , ., , / 7TO; 7TO; \ , ... , / 2lTU> TTLU TTOJ \ 

COSh \k+ " +(l+2cos(^))cosh (_ + —J +(2+2cos(vr J ))cosh ^_ + — + — J = 0. (3.35) 

Notice that if lo is a solution of this equation then so is —Co, as must be the case with quasinormal modes. 
To recover the RN quasinormal frequencies we first write out the equation as 

1 + k c + (1 + 2cos(7rj)) ( e"^~ + e fe c J + (2 + 2 cos(7rj)) ( e~ + ~^ + e k ~ k c J = 

and next take the limit as A — ► + . In this limit, we have Rc ~ A~2 , and hence Rc — > +oo; therefore 

fee = J/'CRc) ~ -Ai?c = ~ -» 0". 

2 /tc 
If we assume that Re(u;) > 0, we see that the equation reduces to 

+ (1 + 2cos(vrj')) + (2 + 2 cos(7rj))e~^ = 0, 

which is exactly equation ( |3.20| ). Therefore, the RN black hole is not a singular limit of the RN dS black 
hole, as far as the quasinormal modes are concerned. The reason for this is basically the same as for the 
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Schwarzschild dS case: the topology of the tortoise at the origin does not change as we take the cosmological 
constant to vanish. 

For vector type perturbations, the potential for the Schrodinger-like master equation near the origin 
is of the form 

with j = 2d-1 - Repeating the same argument, one ends up with the same equation ( |3.35| ) for this new 
value of j. Since 

3d -7 d-3 



2d -5 2d -5 

we see that this equation is exactly the same as in the j = ^3 case, and consequently we end up with 
the same quasinormal modes. 

Some remarks are due for d = 4 and d = 5, as in these cases the Stokes lines are not quite of the 
general form we assumed in the previous calculation. We start by studying the behavior of the solutions 
of the Schrodinger-like equation (3.2) near infinity. For r ~ 00 one has x ~ xq for a (complex) constant 
xq, and the potential is 

7 2 - 1 



4(x-x ) 2 

with = d — 1, d — 3, d — 5 for tensor type, vector type and scalar type perturbations (see appendix |C]). 
Hence, 



<3?(x) ~ A + y / 27TUj(x - Xq) Jjoo (u){x - Xq)) + A_ y/2iru(x - Xq) J_joo (u){x - Xq)) 

2 " 2 

in this region, where J v (x) represents a Bessel function of the first kind and A± are (complex) integration 
constants. Given expansion ( |3.9| ), we see that $ has monodromy (— l) d around infinity. In the d = 4 
case, one can take advantage of the fact that $ has zero monodromy at infinity to conclude that as one 
follows the Stokes line the (anti-clockwise) monodromy of $ around the fictitious horizon r = R equals 
the (clockwise) monodromy around the real horizons r = R^, r = R~^ and r = Rq; since the same is true 
for the plane wave functions e ±luJX , the whole calculation reduces to what was done above, and we get the 
same formula. Another way to see this is that since x, $ are regular at infinity one can deform the d = 4 
contour into, say, the d = 6 contour in the Riemann sphere C U {00}. 

For d = 5 the situation is more subtle. As shown in Figure [TI], the Stokes line approaches infinity, and 
hence one must close the contour near infinity. In the branch containing point B, the solution is then of 
the form 



3>(x) ~ A + y / 2iruj(x - xq) Jjoo (ui{x - x )) + A_ y / 27rw(x - xq) J_joo {u)(x - xq)) 

2 2 

where /3± = j(l ± joo)> with joo = 4, 2, for tensor type, vector type and scalar type perturbations. As r 
rotates clockwise towards the branch containing point A, by ^, x — xq ~ y- rotates anti-clockwise by the 
same amount, and hence 
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at point B. Whereas the coefficient of e~ lulx is unchanged in this process, it is easy to see that the 
coefficient of e luJX reverses sign for all three values of joo. Therefore, there appears an extra minus sign in the 
corresponding monodromy matching equation, and the condition for nontrivial solutions (£?+, B-, C+, C_ ) 
changes slightly. Carrying out the computation, the condition for quasinormal modes in d = 5 turns out 
to be 



, / TTUJ TTUJ \ , , 2tT . . , / TTUJ TTUJ \ , , 27F . . , / 2lTUJ TTUJ TTUJ \ 

sinh " F ) + (1 + 2 cos( T )} smh (jF + % J + (2 + 2 cos( T )} sinh (iF + F + J = a 

It is easily seen that this equation still has the correct RN limit. 

As one searches the literature in the quest for earlier numerical results which could serve as checks on 
our own results, we find that there has not been much work on this matter. In the particular near-extremal 
situation, where the outer black hole event horizon and the cosmological horizon are nearly coincident, 
an approximation scheme for the d-dimensional RN dS solution was developed in [49|. However, the 



approximation used therein is not expected to hold in the asymptotic limit. Quasinormal modes of spinor 
perturbations in RN dS spacetime were studied in j55|, but this is a type of perturbation we have not 
addressed at all. Because our expressions yield the correct RN limit in any dimension, we believe this is 
strong evidence favoring our general results. However, it would be very interesting to produce numerical 
data concerning higher dimensional RN dS quasinormal modes to match against our analytical results. 

3.3 Asymptotically Anti de Sitter Spacetimes 

]l0| discusses the stability of black holes in asymptotically AdS spacetimes to tensor, vector and scalar 
perturbations. For black holes without charge, tensor and vector perturbations are stable in any dimension. 
Scalar perturbations are stable in dimension four but there is no proof of stability in dimension d > 5. For 
charged black holes, tensor and vector perturbations are stable in any dimension. Scalar perturbations 
are stable in four dimensions but there is no proof of stability in dimension d > 5. As we work in generic 
dimension d we are thus not guaranteed to always have a stable solution. Our results will apply if and 
only if the spacetime in consideration is stable. 

Quantization of a scalar field in AdS was first addressed in where considerable attention was 
given to the question of what are the AdS boundary conditions. In fact, in AdS light rays can reach 
spatial infinity and return to the origin in finite time, as measured by the observer at the origin (crossing 
AdS within half the natural period). This would apparently indicate "reflecting" boundary conditions. 
However, the "walls" of AdS are at [timelike] spatial infinity, and so the concept of "reflecting" boundary 
conditions is actually somewhat obscure. In spite of this, embeddings of AdS spacetime in the Einstein 
Static Universe may lead to Dirichlet or Neumann boundary conditions (on the counterparts of the AdS 
fields) which, however, on the AdS fields themselves, will look rather complicated (these are not "reflecting" 
boundary conditions from a pure AdS point of view, see fl5q| ). As it turns out, the only sensible boundary 
condition to impose on quasinormal modes is the usual incoming waves at the black hole event horizon 
and the new requirement of vanishing of the wave function at infinity. 

3.3.1 The Schwarzschild Anti de Sitter Solution 

We now compute the quasinormal modes of the Schwarzschild AdS ci-dimensional black hole. This cal- 
culation was first done, for the case of dimension d = 4 and large black holes, in flq|. Again we consider 
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solutions of the Schrodinger-like equation ( |3.2| ) in the complex r-plane. We start by studying the behavior 
of <£(x) near the singularity r = 0. In this region, the tortoise coordinate is 



„d-2 



x ~ 



2{d-2)^ 

and the potential for tensor and scalar type perturbations is 



V[r(x)] 



Ax 2 



with j = (see appendix |C|). The solution of the Schrodinger-like equation in this region (with j ^ 0) is 
therefore 



<J>(x) ~ B + V2ttujx J 1 (ujx) + £>_ V 2ttujx J_i (ujx) , 

2 2 

where J y {x) represents a Bessel function of the first kind and B± are (complex) integration constants. 

We now wish to examine the Stokes line, i.e., the curve Im(wx) = 0. For the the Schwarzschild AdS 
solution, however, it is no longer true that Im(u;) 3> Ke(cj) for the asymptotic quasinormal modes: instead, 
one finds that ojxq is asymptotically real, where x ~ a^o for r ~ oo (see appendix |C|; this quantity is well- 
defined, as x has zero monodromy around infinity, and hence one can choose the branch cuts arising from 
the branch points at the horizons to cancel among themselves). Interestingly, this is exactly the condition 
that the Stokes line Im(wx) = should extend out to infinity. From the general expression for the tortoise 
(see appendix |C|), 



d-l 



1 



71=1 



one finds 



d-l 



x 



n=l 



1 

Rn 



(with an appropriate choice of branch cuts). For instance, for large black holes, where one has the approx- 
imate expression 



R r . 



Ia7 



2ni 

ex P (— (n-1) 



one finds 



x 



AT H sin 3=j 



cxp 



ITT 



d-i r 



where Th is the Hawking temperature for large Schwarzschild AdS black holes. In general, however, there 
is no analytic solution for xo- Because ujxq is asymptotically real we will have arg(u) = — arg(xo) = — #o 
and therefore the function Im(o;x) will in general be multivalued around all horizons. To bypass this 
problem, we choose a particular branch and simply trace out the Stokes curve Im(a;x) shifting the branch 
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cuts so that it never hits them. Note that the behavior of e ±lu)X will still be oscillatory along the curve. In 
a neighborhood of the origin, the above relation between x and r tells us that Im(cjx) = for 

tPQ | inn 

r = p e d - 2 d - 2 , 

with p > and n = 0, 1, . . . , 2d — 5. These are half-lines starting at the origin, equally spaced by an angle 
of -^2- Notice that the sign of tux on these half-lines is (— 1)" , in other words, starting with the half-line 
corresponding to n = 0, the sign of ujx is alternately negative and positive as one goes anti-clockwise 
around the origin. 

From the asymptotic expansion ( |3.6[ ) we see that 



<I>(x) ~ 2B + cos (ujx — a + ) + 2i?_ cos (ujx — a_) 

= (B+e~ ia + + B^e~ ia -) e iuJX + (B + e ia + + B_e ia ~) e~ iujx (3.36) 

in any one of the lines corresponding to positive ujx, where again we have defined 

*± = \(l±j). 

We shall use this asymptotic expression for matching. This matching must be done along the Stokes line 
ujx S M, so that neither of the exponentials e ±luJX dominates the other. 

To trace out the Stokes line Im(wi) = let us first observe that we already know its behavior near the 
origin. Furthermore, this is the only singular point of this curve: indeed, since x is a holomorphic function 
of r, the critical points of the function Me(x) are the zeros of 

dx 1 r d ~ 3 

dr ~ f(r) ~ -\r d ~ l + r d ~ 3 - 2^ 

(i.e., r = only). To understand the behavior of the Stokes line near the singularities, we notice that if 
one follows our procedure of shifting the branch cuts, the curve 



Im (a \og(z)) = 

is the curve 



a\og(z) = p O z = e ip e ir}P (pel) 

(with a = This is a spiral that approaches the singularity z = 0, except in the case where £ = 0, 

i.e., when a is purely imaginary. Therefore, generically one expects the Stokes line to hit all d — 1 horizons 
in this spiraling fashion. Recall however that at least one branch must extend out to infinity. In our choice 
of branch cuts, we have arranged things so that this branch corresponds to n = 1; therefore the preceding 
branch (corresponding to n = 0) will hit the black hole horizon r = Rh- Notice also that the formula 

UJ 

UJX ~ UJXq 



X\r 

for r ~ oo implies that the argument of r along the branch which extends to infinity must be asymptotically 
equal to the argument of uj, i.e., —9q. Therefore, one would guess that the Stokes line is as depicted in 



Figure 12. We have verified this guess with a numerical computation of the same Stokes line, and this 



- 43 - 





Im , 


1 . 7! 




72 
• 




contour 


Stokes line 
s. Me 








"71 *~ 




"? 


A 




• 

72 




7i 





Figure 12: Stokes line for the Schwarzschild Anti-de Sitter black hole black hole, along with the chosen contour for 
monodromy matching, in the case of dimension d = 6. 




Figure 13: Numerical calculation of the Stokes lines for the Schwarzschild Anti-de Sitter black hole in dimensions 
<i = 4, d = 5, d = 6 and d = 7. Different shadings also illustrate the various horizon singularities and branch cuts 
(note that these branch cuts are not necessarily equal to the ones used for the calculation in the main text). 

is indicated in Figures [l^ and [14]. It should be noted that due to the branch cuts (which can be readily 
identified in the figure) the spirals at the singularities are not so clearly depicted in the numerical result. 
Now, for r ~ oo we have 



4(x-x ) 2 

where = d — 1, d — 3, d — 5 for tensor type, vector type and scalar type perturbations (see appendix |C]). 
Consequently, 

<I>(x) ~ C + J2-kuj(x - x ) Jjvo (uj(x - x )) + C7 Y / 2vra;(x - x ) J_j°° (^( x ~ x o)) 

2 2 

for r ~ oo. The boundary condition $ = at r = oo (discussed at the beginning of this section) requires 
that C- = 0. Hence, 
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Figure 14: Numerical calculation of the Stokes lines for large Schwarzschild Anti-de Sitter black holes in d — 4, 
d = 5, d = 6 and d = 7. The lines are as before, but horizon singularities are now located at (d — 1) roots of unity. 



~ C+J2itijj{x - xq) Jjoo {u{x - Xq)) 

2 

for r ~ co. This is to be matched with expansion (|3.36| ) along the branch of the Stokes line extending out 
to infinity. To do so, we notice that for r ~ oo one has 

LO{X - X ) ~ -TTT-' 

| A|r 

and that this quantity is negative on this branch. Consequently, expansion ( |3.16| ) yields 
in the same limit, where 

We conclude that B + , Z?_ must satisfy 

(B + e- ia + + B_e~ ia ~) e ^o e -i/3+ _ [B + e ia + + B_e ia -) e~ iuJX ° e i(3 + . (3.37) 

Expansion ( |3.9| ) implies that as one rotates from the branch containing point B to the branch containing 
point A (see Figure |l2|) (through an angle of i n r > corresponding to an angle —it in x), one has from 

V27re- 7ri ojx J, L (e'^iox) = e^r( 1±j) V2irux J, L (ujx) ~ 2e~ 2ia± cos{ujx - a±) 
that $ changes to 

*(a?) ~ 2B + e~ 2ia + cos (-ujx - a+) + 25_e _2ia - cos (-ujx - a_) = 
= (B +e - ia + + B_e~ ia -) e iujx + (£+e~ 3ia + + B_e~ 3ia -) e~ iuJx . 

This form of the solution can be propagated along the corresponding branch of the Stokes line which 
approaches the event horizon, and where we know that $(x) must behave as e luJx . Consequently we obtain 
the second condition on B + , B- , as 

B +e - 3ia + +B_e~ 3ia - = 0. 
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The two conditions on these coefficients can only have nontrivial solutions if and only if 



-3iau 



-3ia- 



44> 



-3iaq 



-3i«- 



sin (a + + /?_)_— wxo) sin (a_ + j3 + — loxq) 



0. 



Again, this equation is automatically satisfied for j = 0. We must therefore take j nonzero and then 
compute the limit as j — ► 0. This amounts to writing the equation as a power series in j and equating to 
zero the first nonvanishing coefficient, which in this case is the coefficient of the linear part. Thus we just 
have to require that the derivative of the determinant above with respect to j be zero for j = 0. This is 



3i7T , 



3mt , 



4 " 4 ° 

sin (I + /3 + - ivx ) sin (f + fi + - ujx ) 
from where we obtain our final result as 



+ 



3i7T 

e * 



3iir 

e 4 



f cos (f + P+ - ojx ) -fcos(f +(3+-u;x ) 



TT 



tan ( - + /9 4 



OJXq 



3 ^ 



wxq = ^ + /3+ — arctan f - ] + n-zr (n £ N). (3.38) 



Notice that in this asymptotically AdS case there was no need to do a monodromy calculation, all we had 
to do was follow the Stokes line from infinity to the black hole event horizon. 

For vector type perturbations, the potential for the Schrodinger-like master equation near the origin 
is of the form 



V[r(x)] 



3 



1 



Ax 2 



with j = 2 (see appendix)^). Repeating the same argument, one ends up with the same equation, except 
that now a± = ? ± 5. Since a+ — a_ = it, this equation is still automatically satisfied for j = 2; taking 
the limit as j — > 2 one obtains 



3vr 



UJXQ 



[ arctan 1-1 + nvr (n G N). (3.39) 

The difference of ^ between this equation and equation ( |3.38| ) exactly compensates the difference between 
the values of (3 + for the different types of perturbations, so that all types of perturbation yield exactly the 
same asymptotic quasinormal frequencies, as given by 



TT 



LUXq 



(d+1) 



arctan 



+ nix (n G N). 



Some remarks are due for scalar type perturbations in dimensions d = 4 and d = 5. Indeed, in the 
previous calculation we have assumed that j'oo is positive. This is not true for scalar type perturbations in 
dimensions four and five. In dimension d = 4 the calculation still goes through as long as one sets joo = +1 
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(rather than —1) as all one needs is = 1. Consequently the formula for asymptotic quasinormal 
frequencies of scalar type perturbations in dimension d = 4 gets changed to 

coxa = — arctan ( - ) + nir (n £ N). 



In dimension d = 5 the situation is more subtle. Here joo = and strictly speaking the two Bessel functions 
in the solution at r ~ oo are not linearly independent. As usual, we choose joo slightly positive and carry 
through with the calculation, taking the limit joo — ► 0. In this case the solution at r ~ oo automatically 
vanishes as i^s= i s 

positive. Thus, there is no constraint arising from infinity on the solution to the 
master equation, and all one is left with is the constraint at the black hole horizon. This condition alone is 
not enough to quantize the asymptotic frequencies. Thus, in five dimensions, the asymptotic spectrum of 
Schwarzschild AdS scalar type perturbations is continuous, i.e., u £ C. This is a rather unexpected result. 
We believe that the consequences of this result should deserve further studies. 

We now need to review the literature concerning the calculation of asymptotic quasinormal frequencies 
in the Schwarzschild AdS spacetime, in order to compare our results to earlier work done on this subject. 
Quasinormal modes of Schwarzschild AdS black holes were first addressed in p% j5S| , |59| , |60| although the 
emphasis lied in the lowest frequency modes (with |j(J] focusing on spinor perturbations). It was shown in 
|58| that quasinormal frequencies of Schwarzschild AdS black holes have a direct interpretation in terms of 
the dual CFT: large AdS black holes correspond to CFT thermal states, so that quasinormal frequencies 
correspond to poles in the retarded Green's function of the dual CFT. Because in this context the decay 
of the perturbations is describing the return to thermal equilibrium, the inverse of the gap is naturally 
associated with the time scale for approaching this thermal equilibrium and the inverse of the offset is 
naturally associated with the oscillation time scale. It was moreover shown in |58|] that these time scales 
for the approach of thermal equilibrium are universal, in the sense that all scalar fields with the same 
angular dependence will decay at this rate. Due to the AdS/CFT correspondence, this work sparkled a 
series of investigations on AdS asymptotic quasinormal frequencies, which naturally concentrated in the 



five dimensional case (see, e.g., |2^, |23|, 61, p2|). 

For the case of dimension d = 4, the first numerical results for the asymptotic quasinormal frequencies 
were published in [63]. These authors found that scalar perturbations are isospectral with both odd and 



even parity gravitational perturbations, and they also found the existence of modes with purely imaginary 
frequency. Later, an extensive study of asymptotic quasinormal frequencies for Schwarzschild AdS black 



holes in d = 4 was done in [27], and numerically produced some numbers which exactly match our analytical 
prediction. While the authors of [^7]] found that the real part of the frequency mode increases with the 
overtone number, n, in what seems to be a characteristic particular to AdS space, they also found the 
following numerical data (having set |A| = 1). For large black holes, with Rh 3> 1, the result for scalar 
field perturbations (corresponding to tensor type perturbations) was found to be (we have translated the 



results in 27] to our own conventions on quasinormal boundary conditions) 



uJ n ~ Rh [(1.299 + 2.25i) n + (0.557 + 0.423i)] , 
while for gravitational odd perturbations (corresponding to vector type perturbations) they found 

w n ~ R H [(1.299 + 2.25a) n + (0.58 + 0.42*)] , 
and for gravitational even perturbations (corresponding to scalar type perturbations) 
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w n ~ Rh [(1-299 + 2.25?) n + (0.581 + 0.4H)] . 

For intermediate black holes, with Rh ~ 1, the result for scalar field perturbations (corresponding to tensor 
type perturbations) was found to be 

oo n ~ (1.97 + 2.35i) n + (0.79 + 0.35i) , 
while for gravitational odd perturbations (corresponding to vector type perturbations) they found 

uj n ~ (1.97 + 2.35i) n + (0.93 + 0.32i) , 
and for gravitational even perturbations (corresponding to scalar type perturbations) 

uj n ~ (1.96 + 2.35*) n + (2.01 + IM) . 

For small black holes Rh <C 1 and in this Rh — ► limit the frequencies were shown to approach the pure 
AdS frequencies (which we address later in the present paper), a result which had been first shown in 
|p4| . The numerical results for scalar field perturbations (corresponding to tensor type perturbations) were 
found to be 

oo n ~ (1.69 + 0.57«) n + (2.29 + 0.46i) , 
while for gravitational odd perturbations (corresponding to vector type perturbations) they found 

uj n ~ (1.69 + 0.59i) n + (2.49 + 0.06z) , 
and for gravitational even perturbations (corresponding to scalar type perturbations) 

Lo n ~ (1.61 + 0.6») n + (2.7 + 0.37i) . 
Let us also recall that the case of large Schwarzschild AdS black holes in dimension d = 4 was first solved 



analytically in [p.q] , fully agreeing with the numerical results in [27], as listed above. 

For the case of dimension d = 5, different methods were applied in [22, 23, to study scalar field 
perturbations (corresponding to tensor type perturbations). The analysis of [22, 23 1 relied on the AdS/CFT 
correspondence, their goal being to find poles in the retarded gauge theory thermal correlators, which in 
turn is equivalent to computing quasinormal frequencies in the dual gravitational background. Their results 
were obtained numerically, via use of the Heun equation. The asymptotic quasinormal frequencies found 
for large AdS black holes were (we have translated the result to our own conventions on quasinormal 
boundary conditions) 

uo = 2imTn (1 + i) + <^0i 

where = ttTh (1.2139 + 0.7775i) and Th is the Hawking temperature in this spacetime. The analysis of 
p5j was similar, in the sense that they also computed asymptotic quasinormal frequencies for large AdS 
black holes in dimensions d = 3 and d = 5 via an approximation of the wave equation by the hypergeometric 
equation, and in the high frequency regime. However, they differed from |22|, ^] as they used a method of 
monodromy matching rather than a numerical scheme. In dimension d = 5 they find a result which is in 
complete agreement with the results of |22|, E3[. 
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Let us now compare with our own results. As we have seen, our Schwarzschild AdS quasinormal 
equation depends on a parameter, xq, for which there is no general analytic solution. However, one can 
find an analytical formula for xq in the case of large black holes, and we shall begin there. For these black 
holes the asymptotic quasinormal frequencies are given by (notice that one needs to be careful with the 
choice of logarithmic branch cuts when computing xq) 

I 7T \ it f 7T \ iT ( 7T ( % 

u> = AttiiTh sin ( - — - I + ATjj sin I - — - j e^- 1 I — (d + 1) — arctan ( - 



so that their leading term is located along the direction of a (d — 1) root of unity. In dimension d = 4, and 

47T I 



using the relation Tjj = j-\\\Rh, it is immediate to write down 



uJ 



W n H + + — (0) 

= (1.299 + 2.250?) n + (0.573 + 0.419i) 
for tensor (corresponding to scalar field) and vector type perturbations, and 

a 'l + t rf\ n+ *£(l +t rf\(*E.- arctan '' 



\\\R H 4V ) 4vrV J \ 4 \3 

= (1.299 + 2.250i) n + (1.222 + 1.544i) 

for scalar type perturbations. Similarly, in dimension d = 5 it is immediate to write down for tensor and 
vector type perturbations: 

uj ,„ ., ,„ ., / 2 ( i 

= 27m (1 + i) + -/r (1 + i) I 1 arctan I — 

T H \ 7T V 3 

= 2vm (1 + i) + 7T (1.22064 + 0.77936i) . 

In either case, for large Schwarzschild AdS black holes, our analytical results are in full and complete 
agreement with earlier numerical calculations to a great degree of accuracy, as long as one restricts to 
tensor (scalar field) or vector type perturbations. For scalar type perturbations in dimension d = 4 our 
results yield complete agreement with the gap, but not the offset. We believe further numerical studies 
should be performed in order to obtain better estimates on the offsets. 



Let us next turn to the case of intermediate black holes as discussed in ]27j, where Rh ~ 1. Here, we 
first proceed numerically in order to determine xq (paying special attention to the choice of logarithmic 
branch cuts) and then plug the numerical result for xq into our general expression for the asymptotic 
quasinormal frequencies. Using the values of [27] for cosmological constant, black hole mass and horizon 
radius it is immediate to obtain, from the Schwarzschild AdS quasinormal equation, 

to = (1.969 + 2.350f) n + (0.752 + 0.370i) 
for tensor (corresponding to scalar field) and vector type perturbations, and 

co = (1.969 + 2.350f) n + (1.737 + 1.545i) 
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for scalar type perturbations. This result is again in full and complete agreement with the numerical calcu- 



lations of intermediate black holes in [27], to a great degree of accuracy and for all types of perturbations. 

Finally, we consider the case of small black holes, where Rh — * 0. If we take this limit directly in 
the Schwarzschild AdS quasinormal equation, it is not too hard to see that only the complex horizons will 
contribute in the calculation of xq which becomes 

TT 

Xq 



2y/\\\ 

with the asymptotic frequencies immediately following as 

rr—7 ( d+1 1 (l 

uj = 2\/\\\ n H arctan — 

V ' ' V 4 TT V 3 

As we shall later show, in section 4, this yields the correct value for the gap of pure AdS normal frequencies 
(while the offset comes out incorrect, but there is no reason why it should be preserved in this limit). Once 
again this confirms the numerical expectations in [27, 64] where it was observed that the Schwarzschild 



AdS asymptotic frequencies approached the pure AdS frequencies in the Rh — * limit. We can moreover 
check the numerical prediction of our analytical formula. Using the values of [27] for cosmological constant, 



black hole mass and horizon radius, it is immediate to obtain for small black holes 

co = (1.696 + 0.571i) n + (2.182 + 0.529i) 
for tensor (corresponding to scalar field) and vector type perturbations, and 

u) = (1.696 + 0.571t) n + (3.030 + 0.815t) 
for scalar type perturbations. This result is again in full and complete agreement with the numerical 



calculations for the gap of small black holes in [27|, to a great degree of accuracy and all types of perturba- 
tions. The agreement with the numerical calculations for the offset is reasonable, where we believe further 
numerical studies should be performed in order to obtain better estimates on the offsets. 

In summary, our results completely match what was previously known in the literature. However, it 
would be very interesting to produce further numerical data concerning higher dimensional Schwarzschild 
AdS quasinormal modes to match against all our analytical results. 

3.3.2 The Reissner Nordstrom Anti de Sitter Solution 

We now compute the quasinormal modes of the RN AdS <i-dimensional black hole. Again we consider 
solutions of the Schrddinger-like equation Q3.2j ) in the complex r-plane. We start by studying the behavior 
of <E>(x) near the singularity r = 0. In this region, the tortoise coordinate is 



r 2d-5 



x ~ 



(2d - 5) 2 ' 

and the potential for tensor and scalar type perturbations is 

with j = 2d-<~, ( see appendix |C|). The solution of the Schrodinger-like equation in this region (with j ^ 0) 
is therefore 
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~ B + V2ttujx J i (tux) + B- V 2ttujx J _j_ (ujx) , 

2 2 

where J v (x) represents a Bessel function of the first kind and B± are (complex) integration constants. 

We now wish to examine the Stokes line, i.e., the curve Im(wx) = 0. For the the RN AdS solution, 
however, it is no longer true that Im(cj) 3> Re(a>) for the asymptotic quasinormal modes: instead, one finds 
that ojxq is asymptotically real, where x ~ xq for r ~ oo (see appendix |C| this quantity is well-defined, as 
x has zero monodromy around infinity, and hence one can choose the branch cuts arising from the branch 
points at the horizons to cancel among themselves). Interestingly, this is exactly the condition that the 
Stokes line Im.(ux) = should extend out to infinity. From the general expression for the tortoise (see 
appendix [C]) , 



one finds 




(with an appropriate choice of branch cuts). Thus, because ujxq is asymptotically real, we will have 
arg(cj) = — arg(xo) = — 0$ and therefore the function lm.{ux) will in general be multivalued around all 
horizons. To bypass this problem, we choose a particular branch and simply trace out the Stokes curve 
Im(wx) shifting the branch cuts so that it never hits them. Note that the behavior of e^ luJX will still be 
oscillatory along the curve. In a neighborhood of the origin, the above relation between x and r tells us 
that Im(ws) = for 

■ in-K 

if — n g2d— 5 2d— 5 ; 

with p > and n = 0, 1, . . . , Ad — 11. These are half-lines starting at the origin, equally spaced by an 
angle of 2 d-5 • Notice that the sign of ujx on these lines is (— l) ra ; in other words, starting with the line 
corresponding to n = 0, the sign of cox is alternately positive and negative as one goes anti-clockwise 
around the origin. 

From the asymptotic expansion fl3.6|) we see that 

<I>(x) ~ 2B + cos {ojx — a + ) + 2i?_ cos (lux — q_) 

= (B + e~ ia + + B_e~ ia -) e iuJX + (B + e ia + + B_e ia -) e~ iu3X (3.40) 
in any one of the lines corresponding to positive ujx, where again 

We shall use this asymptotic expression for matching. This matching must be done along the Stokes line 
ujx G M, so that neither of the exponentials e ±tuJX dominates the other. 

To trace out the Stokes line Ihi(ujx) = let us first observe that we already know its behavior near the 
origin. Furthermore, this is the only singular point of this curve: indeed, since x is a holomorphic function 
of r, the critical points of the function Me(x) are the zeros of 
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dx 1 r 



2d-6 



dr f(r) —\ r 2d-4 + r 2d-6 _ 2nr d -3 + 2 

(i.e., r = only). Generically one expects the Stokes line to hit all d — 1 horizons in a spiraling fashion. 
Recall however that at least one branch must extend out to infinity. In our choice of branch cuts, we have 
arranged things so that this branch corresponds to n = 2; therefore the preceding branch (corresponding 
to n = 1) will hit the black hole horizon r = R^. Notice also that the formula 

u 

UJX ~ UJXq 



X\r 

for r ~ oo implies that the argument of r along the branch which extends to infinity must be asymptotically 
equal to the argument of u, i.e., —9q. Therefore, one would guess that the Stokes line is as depicted in 



Figure 15. We have verified this guess with a numerical computation of the same Stokes line, and this 
is indicated in Figures |l^ and [l?]. It should be noted that due to the branch cuts (which can be readily 
identified in the figure) the spirals at the singularities are not so clearly depicted in the numerical result. 
Now, for r ~ oo we have 

7 2 - 1 
Joo 1 



V 



4(x-x ) 2 

where joo = d — 1, d— 3, d— 5 for tensor type, vector type and scalar type perturbations (see appendix Q). 
Consequently, 



<3?(x) ~ C + ^2ttu(x - x ) J ho (uj(x - x )) + C_ \J2ttuj{x - x ) J_joo {w{x - x )) 



2 2 




Figure 16: Numerical calculation of the Stokes lines for the Reissner-Nordstrom Anti-de Sitter black hole in 
dimensions d = A, g? = 5, d = 6 and d = 7. Different shadings illustrate various horizon singularities and branch cuts 
(note that these branch cuts are not necessarily equal to the ones used for the calculation in the main text). 




Figure 17: Close-up, near the origin of the complex plane, on the numerical calculation of the Stokes lines for the 
Reissner-Nordstrom Anti-de Sitter black hole in dimensions d = 4, d = 5, d = 6 and d = 7. 

for r ~ oo. The boundary condition = at r = oo (discussed at the beginning of this section) requires 
that C_ = 0. Hence, 

$(x) ~ C+ y / 2mj(x - xo) Jjoo {<jj{x - xq)) 

2 

for r ^ co. This is to be matched with expansion (|3.40| ) along the branch of the Stokes line extending out 
to infinity. To do so, we notice that for r ~ oo one has 

LO{X ~ X ) ~ 7T~j J 

| A|r 

and that this quantity is negative on this branch. Consequently, expansion ( |3.16| ) yields 

$(x) ~ c + e il3 +e luJ ^- x ^ + c+e-^+e-^-^ 

in the same limit, where 

P+ = ~(l+j o0 ). 

We conclude that B + , Z?_ must satisfy 

(B + e~ ia + + B^e~ ia -) e^e' 1 ^ = (B + e ia+ + B„e ia -) e'^e^. (3.41) 
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Expansion implies that as one rotates from the branch containing point B to the branch containing 
point A (see Figure |15|) (through an angle of ~ 2d-b m r ' corresponding to an angle — tt in x), one has from 



\/2vre- 7ri cJ3; J ±3 _ (e~ %i ux) = e^r (1±j) \/2ttujx J ±l (ux) ~ 2e~ 2ia± cos(wx - a±) 
that <!> changes to 

*(a?) ~ 2 J B+e~ 2iQ + cos (-ux - a+) + 2B_e~ 2iOL - cos (-wx - a_) = 
= (B + e~ ia + + B^e~ ia -) e iujx + (B + e~ 3ia + + B_e -3ia -) e~ luJx . 

This form of the solution can be propagated along the corresponding branch of the Stokes line which 
approaches the outer event horizon, and where we know that $(x) must behave as e luJX . Consequently we 
obtain the second condition on B + , B- , as 

B + e- 3ia + + B^e~ 3ia ~ = 0. 
The two conditions on these coefficients can only have nontrivial solutions if and only if 

g— 3io+ g— 3ia_ 

g-ia+g-i/3+giajxo _ g«a+ e i/3+ e -iu>xo e~ ia ~ e~^ + e iuJX ° — e ia ~ e i @ + e~ iuJX ° 



O sin(vrj) e -^ + e"^ sin ( ^ ) e*"*" 



2 J 



from where we obtain our final result as 



ux = log \2 cos (^\\ +j+ 13+ + niT (iiEN). (3.42) 

Notice that in this asymptotically AdS case there was no need to do a monodromy calculation, all we had 
to do was follow the Stokes line from infinity to the black hole event horizon. 
For vector type perturbations, one has 

3d - 7 „ d-S 

J 



2d -5 2d -5 

(see appendix |C|) . Therefore it is easily seen that the values of ujxq — j3 + for the corresponding quasinormal 
modes are shifted by 5 relative to the quasinormal modes for tensor and scalar type perturbations. This 
difference exactly compensates the difference between the values of /?+ for the different types of pertur- 
bations, so that all types of perturbation yield exactly the same asymptotic quasinormal frequencies, as 
given by 

ojxq = J (d+ 1) + ^log ^2 cos fyj J + n7T ( n G N )- 

Some remarks are due for scalar type perturbations in dimensions d = 4 and d = 5. Indeed, in the 
previous calculation we have assumed that joo is positive. This is not true for scalar type perturbations in 
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dimensions four and five. In dimension d = 4 the calculation still goes through as long as one sets = +1 
(rather than —1) as all one needs is = 1. Consequently the formula for asymptotic quasinormal 
frequencies of scalar type perturbations in dimension d = 4 gets changed to 

ujxq = ^ + log ^2 cos (^Y^j 1 + n7T (neN). 

In dimension d = 5 the situation is more subtle. Here joo = and strictly speaking the two Bessel functions 
in the solution at r ~ oo are not linearly independent. As usual, we choose joo slightly positive and carry 
through with the calculation, taking the limit joo — * 0. In this case the solution at r ~ oo automatically 
vanishes as is positive. Thus, there is no constraint arising from infinity on the solution to the master 
equation, and all one is left with is the constraint at the black hole horizon. This condition alone is not 
enough to quantize the asymptotic frequencies. Thus, in five dimensions, the asymptotic spectrum of RN 
AdS scalar type perturbations is continuous, i.e., to E C. This is a rather unexpected result. We believe 
that the consequences of this result should deserve further studies. 

Let us now review the literature on asymptotic quasinormal frequencies for the RN AdS spacetime, 
so that we can compare our results with earlier work done on this subject. The stability of RN AdS black 
holes was addressed in [65, pfifl , in a very different way from the analysis in These authors focused on 



the d = 4 case and found that, while instabilities in the linearized theory are harder to find in this RN AdS 
case, one can still decide on the stability of the solution if one tries to uncover thermodynamic instabilities. 
Such instabilities would correspond to the onset of Bose condensation (or other sort of phase transitions) 
in the dual field theory. Interestingly enough, at the perturbation theory level unstable modes of the RN 
AdS4 black hole are not found in metric fluctuations at linear order, but rather in the gauge fields and 
scalars of N = 8 gauged supergravity. [35, 36 argue that this instability is in the same universality class 



as the one of instabilities where the horizon does fluctuate. As to the metric, it also fluctuates, but only at 
subleading order. The thermodynamic instabilities can only be studied for large black holes and, in this 
limit, they do coincide with the dynamical instabilities. For a discussion of this sort of instabilities, in a 
broader context, we refer the reader to |67[| . All in all, one should have the above discussion of instabilities 
in mind when addressing RN AdS black holes. 

The first numerical results for the asymptotic quasinormal frequencies of RN AdS black holes were 



computed in [63|, in the particular case of dimension d = 4 and focusing on scalar field, electromagnetic 
field and gravitational field perturbations. The most significant finding of these authors was that purely 
damped RN AdS modes behave in a very peculiar way: their damping seems to go to infinity in the 
extremal black hole case, suggesting the possibility that extremally charged RN AdS black holes may 
be marginally unstable to electromagnetic and gravitational perturbations (unless, of course, the mode 
amplitude goes to zero in this extremal limit). This question was later solved in p8] . Indeed, |68|] presents 
an extensive numerical study of quasinormal frequencies for massless scalar fields in d = 4 RN AdS 
spacetime, corresponding to tensor type perturbations in our results. They find that as one increases the 
charge, at first the scalar field decays exponentially and oscillates but for > # cr iti C ai = O.38950 max the 
decay becomes purely exponential. So, for 6 > ^critical the asymptotic scalar perturbations are dominated 
by non-oscillatory modes and |Imu;| keeps decreasing monotonically and smoothly up to = mSuX . As 
the black hole becomes extremal, |Ima;| — > 0. However, in the extremal limit the decay of the scalar field 
perturbation changes from exponential to power law, and the extremal RN AdS black hole is thus stable 
to scalar field perturbations (not marginally unstable as worried in |63| and mentioned above). For the 
higher modes both Mew and Imw increase with overtone number, n, although as one increases charge 0, 
Mew increases slower and Imw increases faster, with n. In the large black hole regime, with |A| = 1 and 
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R^j = 100, asymptotic frequencies become evenly spaced and as n — ► oo (we have translated the results 
adapting the conventions of [6q] to our own) 



u) n+ i - Lo n ~ 129.9 + 225i, = 
w n+ i - cj„ ~ 120.6 + 233i, = 0.150 
u n+1 -oj n ~ 111.1 + 256i, = 0.30 



max • 



max • 



As one increases charge, the real part of the spacing decreases and the imaginary part of the spacing 
increases. We can now match these predictions to our RN AdS quasinormal equation. This is done in a 
similar fashion to what we did in the case of Schwarzschild AdS black holes: one computes xq numerically 
according to its general expression (as there is no analytical solution for xq, and being careful with the choice 
of branch cuts for the several logarithms) and then plugs into our analytical formula for the frequencies 
in order to obtain (using the values of |38| for cosmological constant, black hole mass, charge and horizon 
radius) 

lo = (122.362 + 231.978i) n + (50.871 + 47.297i) , 
when = O.150 max , and also 



uJ 



(110.973 + 255.82i) n + (50.108 + 54.253i) 



when = O.30 max . Our results are in full and complete agreement with the numerical calculations of [68], 
to a great degree of accuracy. However, it would be very interesting to produce further numerical data 
concerning higher dimensional RN AdS quasinormal modes to match against all our analytical results. In 
particular, for the results above, we have produced both gap and offset, where there is only numerical 
data for the gap. A point to mention concerns the critical charge cr itical mentioned above. Unfortunately 



asymptotic quasinormal modes cannot see this critical charge: as described in 68], the critical charge is 
found by following the frequencies of a specific mode, as one varies the charge. For the asymptotic fre- 
quencies n — > +oo and the choice of focusing on a specific mode becomes somewhat obscure. Nevertheless, 
we believe this matter should be further explored within our general setting. 

Let us make one final remark concerning the case above where = and where the result exactly 
matches the Schwarzschild AdS result. One can prove that this must be the case by analyzing the limit 
of the RN AdS quasinormal frequencies when 0^0. Here, it is not too hard to check that the (2d — 4) 
RN AdS complex horizons have the following limit: (d — 1) of these horizons go to the (d — 1) complex 
horizons of Schwarzschild AdS, while the remaining (d — 3) horizons go to zero as 

'0 2 \^ 

when 0^0. This limit is such that these remaining (d — 3) horizons will have zero contribution to xq. 
Thus, in the 0^0 limit the RN AdS xq reduces to the Schwarzschild AdS xq. In this case, the gap of 
the asymptotic quasinormal frequencies of RN AdS spacetime will reduce to the gap of the asymptotic 
quasinormal frequencies of Schwarzschild AdS spacetime (although the offset comes out incorrect, as ex- 



pected in this sort of limits). This shows, in particular, why the numerical result of [68] must match the 
Schwarzschild AdS result when = 0. 
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4. Exact Solutions for (Quasi) normal Frequencies 



While most black hole solutions do not allow for an exact calculation of quasinormal frequencies, there 
are some spacetimes where one can actually perform such an analytic calculation (and without restricting 
oneself to the asymptotic case). We shall address some of these examples in this section; in particular, we 
shall deal with the case where the parameters fi and 9 are both zero but A is unconstrained (thus being 
in a situation where there is no black hole), and while still remaining in the <i-dimensional realm. Let us 
stress that we present these results for the sake of completeness, as they do not describe full fledged black 
holes. Still, they are valid solutions to Einstein's equations and are definitely worth studying within the 
classification problem we address. The asymptotically flat case we shall consider is Rindler space, and it 
is without surprise that we shall find a trivial solution for the quasinormal frequencies (after all, Rindler is 
but flat space). The spacetimes with non-vanishing cosmological constant we will consider are, obviously, 
AdS and dS. Using the analysis in ||, we write down the equations for the several types of perturbations 
and show how they can be exactly solved with quasinormal mode boundary conditions. For AdS, we find 
an unexpected continuous spectrum for scalar perturbations in five dimensions. For dS, we find that there 
are only quasinormal frequencies in odd spacetime dimensions. 

4.1 The Rindler Solution 

Let us begin with the Rindler spacetime solution. It actually happens that this solution does not belong 
to the class of backgrounds analyzed in ||] , and one thus cannot use their perturbation theory formulas. 
Instead of developing a full perturbation theory for this spacetime, we shall simply calculate the quasinormal 
frequencies for a massless scalar field (meaning a field obeying the scalar wave equation in this spacetime). 
As we shall see, the solution to this problem is rather simple and does not justify by itself an extension of 
the analysis in ||. The ^-dimensional Rindler spacetime is described by the metric 

g = —p 2 dt ® dt + dp ® dp + dz\ g) dz\ + • • • + dzd-2 <8> dzd-2, 

where the coordinate p must take positive values. In this geometry, it is immediate to find the massless 
Klein-Gordon equation for a scalar field <j). We Fourier decompose the field in both time and transverse 
space (meaning the transverse coordinates z = (z%,..., Zd-2)) as 

0(p,z,i) = $(p) e -V^, 

with k the wave vector, and use this decomposition in the wave function. It then simply follows 

od 2 ® d$ f 2 2 , 2 n , . . 

P lp T + P lp~ + ( 9 ) ' ( } 

with k 2 = k ■ k. The general solution to ( |4.l| ) is given in terms of the modified Bessel functions of the first 
and second kind, I u {z) and K u (z), as {C\,C2 S C are constants) 



HP) = Ci lUkp) + C 2 K iw (kp). 

Next, in order to define quasinormal mode boundary conditions, one first finds the tortoise coordinate. 
This is a particularly simple case, where x = log p. Starting with the behavior near x = —00, one observes 
that the effective potential V(x) = k 2 e 2x vanishes and the Schrddinger-like wave equation is simply solved 
by plane wave solutions of the type 
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$(x) ~ e 



Even though there is no black hole, we shall nevertheless demand the standard quasinormal mode boundary 
conditions, i.e., an in-going wave near the Rindler horizon x = — oo, i.e., <J?(x) ~ e iajx . Let us see what 
this implies in our general solution above. Near the origin, p = 0, the modified Bessel functions behave as 

Iiu(kp) ~ p iuJ , 

Ktoikp) ~ C+p^ + C-p-™, 

where C + ,C_ G C are specific constants. Because the boundary condition is <£(/?) ~ j0 l£J as p — > 0, one 
discards the modified Bessel function of the second kind Ki U (kp), keeping only h^ikp) as an allowed 
solution. Of course one must still analyze whether this is a good solution near the other boundary. At 
x = +00 the effective potential blows up and the wave equation is certainly not solved by plane waves. 
The natural requirement is to demand for regularity, which in this case translates to the vanishing of the 
wave function. However, the modified Bessel function of the first kind we have kept satisfies 

Iiu(kp) ~ — e kp -> +00 (p -» +00). 

One thus comes to the conclusion that this function is not a good solution, as it does not satisfy the 
boundary conditions. The one modified Bessel function satisfying the necessary requirements at infinity 
describes both in-going and out-going waves at the origin, while the one describing the in-going wave at 
the origin blows up at infinity. In other words, there are no solutions with quasinormal mode boundary 
conditions in Rindler space. This is not an entirely unexpected result as Rindler space is but flat space, 
where we know that there are no quasinormal frequencies (although these are different boundary conditions 
to the usual ones for Minkowski spacetime). 

4.2 The Anti de Sitter Solution 

We now turn our attention to the pure AdS spacetime solution. As a Schrodinger-type problem, this 
background effectively behaves as a bounding box, a sort of an infinite well potential with finite extent. 
This characteristic raises the question of what boundary conditions one should use; these were also earlier 
discussed in section 3. This question was first addressed in p6[ within the context of quantization of 
scalar fields in a four dimensional AdS background. Even though an embedding of AdS spacetime in the 
Einstein Static Universe can lead to both Dirichlet or Neumann boundary conditions (on counterparts of 
AdS scalar fields), the only sensible conditions found in j56| were vanishing of the fields at infinity. This 
issue was further addressed in []69| , |26|| , within the study of the ci-dimensional wave equation for a massive 
scalar field in AdS. There, it was first realized that in d-dimensions the AdS wave equation reduces to a 
hypergeometric equation and is thus exactly solvable. The authors also realized that only normal modes 
exist. Due to the finite extent of spacetime in the tortoise coordinate, it is actually not surprising to 
find that any field evolving in an AdS geometry may be decomposed in a normal mode harmonic analysis 
(the basic difference between normal and quasinormal modes being that the former are a complete set of 
stationary eigenfunctions, while the latter are not). Observe that the normal mode problem is different 
from the quasinormal mode one: for normal modes one requires the vanishing of the wave function at both 
boundaries, while for quasinormal modes one requires plane wave solutions at these same boundaries. In 
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Figure 19: Potentials for AdS vector type perturbations in d = 4 and d = 5. 



what follows, we shall describe how to find both normal modes and normal frequencies for cf-dimensional 
AdS spacetime. Let us note that this problem has been previously addressed for a massless scalar field in 
1 26], and for tensor and vector type perturbations of d = 4 AdS spacetime in [^7]]. Our results naturally 
extend these, to arbitrary perturbations and arbitrary dimension. 

Let us begin with a "tour" of the different Schrodinger-type effective po- 
tentials one can find in AdS, in order to set the stage for the results we shall 
later find. The AdS solution has parameter A < (with /i = = 9), so that 
f(r) = 1 — Ar 2 = 1 + |A|r 2 . In any dimension d one can thus simply find the 
tortoise coordinate as 



0.2 0.4 0.6 0.8 1 1.2 1.4 

Figure 18: Potential for 
AdS tensor type perturba- 
tions in d = 5. 



dp 

W) 



1 



tan 



A 

A| x 



arctan 



A| r , 



when described via the tortoise, as lim^o^M = and lim r ^ +00 x[r] = ^ T^ . The potential for tensor 



The first thing one observes is that [infinite] space becomes a finite interval, 
coise, as lim^o^M = and 

type perturbations is thus written as (see the appendix) 

| A| (21 (£ + d-3) + {d-2) 2 + 2^e(e + d-3)-(d-2)^j cos (2^/\\\ cc)) 



V T (x) 



In its range, x £ (0, 



A| x I cos 



A| x 



WW 



4 sin 2 



), this potential looks like an infinite well potential with rounded edges. In 



Figure [18| we present an illustrative plot, where we take | A| = 1, d = 5 and £ = 1 (as one should recall that 
for tensor type perturbations I > 1). The shape of the potential does not change much as we change its 
parameters. For vector type perturbations, the potential is (see the appendix) 



I A| [21 {£ + d - 3) + (d - 2) (d - 4) + 21 (I + d - 3) cos ( 2y/\ A| x 

Vy(x) = ^ r ^ 

4 sin ( A| x I cos 2 ( y|A| x 

This time around, the potential looks like the infinite well potential with rounded edges only at dimension 



d > 5. In d = 4 the potential is actually different. In Figure 19 we present two illustrative plots, where 
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Figure 20: Potentials for AdS scalar type perturbations in dimensions d = 4, d = 5, d = 6 and d = 7, respectively. 



we take |A| = 1 and £ = 2 (as one should recall that for vector type perturbations t > 2), and sequentially 
d = 4 and d = 5. Observe that as x — > — in four dimensions, the potential does not go to zero but 



rather to |A| £ {£ + 1). Except for the change at dimension d = 4, again the shape of the potential will not 
change much as we change its parameters. We will see that even though the potential in d = 4 is different, 
there will be no substantial change in the normal modes. The same will not happen as we turn to scalar 
type perturbations, where the potential is (see the appendix) 

| A| (it (£ + d - 3) + (d - 4) 2 + 2 (l U + d - 3) + (d - 4) ) cos (2J\X\ x 
V s (x) = ^ ^ 



4 sin (^\/|A| xj cos 2 ^a/|A| x 
In this case, the potential looks like the infinite well potential with rounded edges only at d > 7. In 



d = 4, 5, 6 the potential is different and changing. For the illustrative plots of Figure 20 we shall take 



|A| = 1 and I = 2 (as one should recall that for scalar type perturbations £ > 2), and sequentially d = 4, 
d = 5, d = 6 and d = 7. Observe that as x — > — 7= both in four dimensions (the first figure on the left) 

and in six dimensions (the third figure) the potential does not go to zero. Rather, we have 

lim V s (x)\ d ^ = \\\£(e+1), 



lim V s (x)\ d=e = \X\ (£(£ + 3) + 2). 



WW 

The potentials in d = 4 and d = 6, will not produce substantial changes in the results for the normal 
modes. But the d = 5 potential will, and in fact will lead to a continuous normal mode spectrum, which 
turns out to be an unexpected result. Let us next proceed with the actual computation. 

It turns out that the AdS wave equation can be reduced to a hypergeometric equation, in which case 
it pays off not to start with the potentials as written above, but in the original r coordinate. In this case, 
for the pure AdS spacetime the potential governing generic d-dimensional gravitational perturbations may 
be written as 

^M = (l + lV)(a|A| + ^), 

where, in the expression above, 
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d + 2t - 4 



and a depends on the type of perturbation under study. Explicitly: 



( \d{d — 2), Tensor Type Perturbations, 



j (d — 2) (d — 4) , Vector Type Perturbations, 
k j (d — 4) (d — 6) , Scalar Type Perturbations. 



Changing to the new variable z = sin ^ylAlxJ , with x the tortoise coordinate, the wave equation gets 
transformed into a hyper geometric type of differential equation, 



d 2 § . , 

4z(l-*) TT + 2 l-2z — + 
dz z dz 



d§ 4& 2 z(l- «)-4oz-46(6+l)(l- z) 



$ = 0, 



(4.2) 



where we have introduced the notation cD = 



4z(l - z) 

and where the new coordinate, z, ranges from to 1. 



Because one knows the solution to the hypergeometric differential equation, this yields the solution for an 
arbitrary perturbation in the d-dimensional AdS spacetime. We shall now make a separate analysis of this 
equation, for each of the different type of perturbations. 

We begin with the analysis of normal modes for tensor type perturbations. For these perturbations, 
it is not too hard to find that the solution to (|4.2| ) is 



$(z) = Ci x 2 F X 
+C 2 x 2 F X 



+ 6j -l + d + 2£ 



2 ' 2 
3-d-l 



2-d -2+d+2l 
(Z - 1) 4 Z * + 



3-d-£ + u 5-d~2£ 



2-d A-d-2l 
(Z — 1) 4 Z 4 



(4.3) 



where Ci,C 2 G C are constants and 2 Fi [a,/3, 7I z] is the standard hypergeometric function. Two basic 
properties of 2 F\ [a, (3, 7) z], which shall be needed in the following, are 



lim 2 F\ [a, A7l z \ = 1, 

>0 



lim 2 Fi [a,/3,j\ z] 

Z—>1 



r[ 7 ] rjT-a-fl 
r[ 7 - «] r[ 7 - /?] 



(4.4) 



Next, one imposes boundary conditions for normal modes. As we are restricting ourselves to dimension 
d > 4, the solution to (4.3) which is regular at the origin, z = 0, is the first term on the right hand side 
of (f4.3j). One thus sets C 2 = 0. In order to satisfy boundary conditions at the other boundary, z = 1, one 
still needs to demand for a vanishing wave function at this point. This can only be achieved if 



2F1 



+ Cj -l + d + 2i 
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as the poles of the Gamma function will dominate over the power law singularity in (z — l) - * - . It is simple 
to observe, using (|4.4j) , that the condition above is satisfied when the frequency satisfies the quantization 
condition 



Cj = 2n + d + 



1, 



n = 0,1,2, 



Let us note that we have chosen only positive values for the frequency. It thus follows that the normal 
frequencies for tensor type perturbations of d-dimensional AdS spacetime are given by 



A] (2n + d + l-l), n = 0,1,2,.... 



This result agrees with the one in |^], which is a result for scalar fields. This is expected as the Klein- 
Gordon evolution equation is actually identical to the one for tensorial gravitational perturbations. 

Let us proceed with the analysis of normal modes for vector type perturbations. This follows along 
very similar lines to the previous case. In d = 4 the parameter a in ( |4.2j ) actually vanishes, and this 
equation has for solution (this is the general solution when a = 0, we have not yet set d = 4) 



$(z) = Ci x 2 Fi 
+C 2 x 2 F X 



-2 + d+2£-2u -2 + d + 2£ + 2Cj -l + d + 2£ 



4 ' 4 ' 2 

4- d-2£-2u A-d-2£ + 2u 5-d-2£ 



-2+d+2£ 
Z 4 + 



Z 4 



(4.5) 



This four dimensional case has actually been already solved in [27]. A similar line of thought to the one 
above leads to choosing C2 = and to the quantization of frequency as 

d 



n = 0,1,2,.... (4.6) 
Thus, the normal frequencies for vector type perturbations of four dimensional AdS spacetime follow as 



u) = 2n H h 

2 



(jj= y/\\\ (2n + £ + 2), n = 0,1,2,.... 
For any other dimension d > 4, the solution for vector type perturbations is given by 



Hz) = d x 2 Ft 
+C 2 x 2 F 1 



l + l-Qj l + l + Qj -l + d + 2£ 



2 

4-d 



2 2 
u 4-d-e + u; 5-d-2£ 



4-d -2+d+2t 
(Z — 1) * Z 4 



(z-iy 



4 -d 4-d-2£ 
4 Z 4 



and an analogous line of reasoning leads to setting C 2 = and to the normal frequencies for vector type 
perturbations of cf-dimensional AdS spacetime as 



A] (2n + d + £-2), 71 = 0,1, 2,. 



If we set d = 4 in the equation above we obtain the result we had previously found concerning the four 
dimensional case. Thus, even though the normal mode eigenfunctions differ from d = 4 to d > 4, the 
normal frequencies are all given by the above formula. 
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To end the analysis of AdS spacetime, we now discuss the normal modes for scalar type perturbations, 
again via the same line of reasoning. As we have seen in our description of AdS effective potentials, the 
parameter a in fl4.2|) will vanish when either d = 4 or d = 6. The solution to the wave equation in these 
cases will then be given by (^1]), with the prescribed frequencies ([O]) (the four dimensional case was dealt 
with in |27j]). So, in dimension d = 4 the normal frequencies of scalar type perturbation in AdS spacetime 
are 

uj = (2n + t + 2) , n = 0, 1, 2, . . . , 

while in d = 6 they are given by 

uj = yf\\\ (2n + 1 + 3) , n = 0, 1, 2, . . . . 

Let us proceed with the study of the wave equation when d > 6 (we shall address the special d = 5 case at 
the end). In this case the solution to the wave equation is 



Hz) = Ci x 2 F 1 
+C 2 x 2 Fi 



2 + £-oj 2 + 1 + uj -l + d + 2£ 

2 " 

\5-d- 



2 2 
Qj h-d-t + G) 5-d-2£ 



6-d -2+d+2t 
Z — 1) 4 Z 4 + 



6-d A~d-2l 
[Z — 1) 4 Z 4 



and the standard analysis leads to C 2 = and to the normal frequencies for scalar type perturbations of 
d-dimensional AdS spacetime as (d > 6) 



to = y/\\\ (2n + d + t-3), n = 0, 1, 2, . . . , 

where as usual we have only selected the positive values for the frequency. We also observe that the d = 6 
case is included in the previous expression, but the d = 4 case is not (and thus one needs to apply the 
expression we have previously obtained). In order to conclude our analysis, one must still address the case 
where d = 5. As we have seen in the description of AdS effective potentials, this dimension seems rather 
special as the effective potential is negative in certain regions of spacetime. This will lead to an unexpected 
result for scalar type perturbations in five dimensions. If we use the solution above (where, as we have 
already seen, C 2 = 0) in the particular case of d = 5, we find 



2 + 



UJ 



2 + £ + uj 4 + 21 



3+21 
Z 4 . 



One immediately realizes that given the well behaved solution at the origin, it automatically satisfies 
the other boundary condition of $(z = 1) = 0. Thus, in five dimensions, the spectrum of AdS scalar 
type perturbations is continuous, i.e., uj £ K. This is a rather unexpected result. We believe that the 



consequences of this result should deserve further studies. Let us also point out that it was shown in |p7 |, 
in the particular d = 4 case, that as one progressively decreases the size of a Schwarzschild AdS black hole, 
its quasinormal frequencies approach the pure AdS frequencies. We have made use of this property, as well 
as of our results above, in section 3, when discussing Schwarzschild AdS quasinormal frequencies. 
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4.3 The de Sitter Solution 

The last case of exact solutions we wish to address is that of the pure dS spacetime solution. The analysis 
of wave propagation in d-dimensional dS spacetime closely parallels that for AdS, but with one crucial 
difference: there is now a cosmological horizon and this will have defining consequences in the boundary 
conditions one should pick. Indeed, one now chooses quasinormal mode boundary conditions at the cosmo- 
logical horizon, i.e., an out-going plane wave solution. Given this, while in pure AdS we only found normal 
modes, in pure dS we shall only find quasinormal modes. The question of pure dS quasinormal modes has 
actually been recently subject to some discussion, and also to some confusion, in the literature. For the case 



of dimension d = 4, the analysis we shall follow was presented in detail in [47, 32] . We believe this is a clear 
and rigorous derivation, and the authors found that there are no quasinormal modes in four dimensional 
dS spacetime. The open issue is what happens as we extend these results to d-dimensions. Some authors 
[70] claim there will always be dS quasinormal modes. This is in contradiction with the results in |47], |32|] 
and a closer inspection at [70] reveals that the authors do not actually impose quasinormal mode boundary 
conditions at the cosmological horizon (they demand for vanishing of the wave function), thus falsifying 
their conclusions. Some other authors |7l| claim there will never be dS quasinormal modes. While this 
result seems not to contradict any other earlier authors, truth is that upon closer inspection the authors of 
]7"H are actually looking for real frequencies, which is clearly not the case for quasinormal frequencies, thus 
falsifying their conclusions. Further authors have recently addressed quasinormal frequencies of massive 
scalar fields in pure de Sitter spacetime, in d dimensions [72]. We hope our results will help to settle the 



dust, as we shall show that, in general, even dimensional pure dS spacetime never has quasinormal modes, 
but odd dimensional pure dS always has quasinormal modes. 

As in the previous section, we begin with a quick "tour" of the different Schrodinger-type effective 
potentials one finds in dS spacetime. The dS solution has parameter A > (with fi = = 0), so that 
f(r) = 1 — Ar 2 . In any dimension d the tortoise is 

dp -iarctanh^Ar 



f(p) 



r = —= tanh (\/X x] . 
y/X V ' 

This time around r < Rq = — ;= (see the appendix) so that the tortoise is in R + as lim^o^M = and 
lim r ^R c x[r] = +oo. The potential for tensor type perturbations is thus written as (see appendix B) 

A ( 2i (1 + d - 3) + Id - 2) 2 + 2 ( I (£ + d - 3) - (d - 2) ) cosh ( 2yf\ 

Vt(x) = -± -^-^ ? ' 

4 sinh 2 ( y/\ x ) cosh 2 f \f~X 



In its range x 6 M + , this potential looks like the one plotted in Figure |2l| (where we set A = 1, d = 5 and 
I = 1, as one should recall that for tensor type perturbations I > 1). The shape of the potential does not 
change much as we change its parameters, and at infinity the potential decays to zero. For vector type 
perturbations the potential is 



A ( 2£ (£ + d - 3) + (d - 2) (d - 4) + It {£ + d - 3) cosh ( 2y/X x 

WW = — ^ r - ( 

4 sinh 2 ( VA x ) cosh 2 ( 
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In its range x € M. + , this potential looks very much like the previous one. 
Again, the shape of the potential is basically unaffected as we change its pa- 
rameters and at infinity the potential decays to zero. Finally, the potential for 
scalar type perturbations is (see appendix B) 




Vs(s) 



A (21 (£ + d - 3) + (d - 4) 2 + 2(l (£ + d - 3) + (d - 4) ) cosh x)) 
4sinh 2 x^j cosh 2 (VA x^j 



, In the range x € M + this potential again looks like the previous ones, at infinity 

tensor type perturbations m , . TTT , , ,., , , 

^ _ g decaying to zero. We learn that all potentials are very much alike and they do 

not change substantially as we change dimension (unlike the AdS counterpart). 

Let us then proceed with the quasinormal mode computation. 

As in the AdS case, for pure dS spacetime the wave equation can be reduced to a hyper geometric 

equation, and again it pays off to start with the potentials as written in the original r coordinate. In pure 

dS spacetime, the potential governing d-dimensional gravitational perturbations may be written as 



V'M = (l-^)(-«A + M), 



where a and 6 were defined in the previous section. The dS problem is slightly simpler than the AdS 

one, as the factorization <3?(r) = r b+1 (l — Ar 2 )* 2 ^(r) transforms the dS wave equation for <£(r) in a 
hypergeometric equation for ^(r), namely 



(1 - Ar 



d 2 ^ 



dV 



+ 2 (1 + (1 - Ar 2 ) b - (2 + iCj) Ar 2 ) — - Ar (2 - a + b 2 + b (3 + 2iu) + 3iu - Cj 2 ) ^ = 0. 



dr 2 y ' y ' dr 

Given the hypergeometric solution for ^(r), it immediately follows the general solution for <E>(r) as 



$(r) = d x 2 F 1 
+C 2 x 2 F 1 



3 - Vl + 40 + 2b + 2iQj 3 + VI + 4a + 2b + 2iCj 3 + 2b 



4 ' 4 '2 

1 - VI + 4a - 2b + 2iu 1 + VI + 4a - 2b + 2iu 1 - 26 



Ar 
Ar 2 



r 6+1 (l-Ar 2 )^ + 
r- b (1 - Ar 2 )^ . 



Interestingly enough, with the allowed values for a (which depend on the type of gravitational perturbation), 
one will always get rid of the square root. In particular, one obtains the general expression 



$(r) = d x 2 F 1 
+C 2 x 2 Fi 



3-d + j + 2b + 2iu> 3 + d - j + 2b + 2iuj 3 + 26 



4 ' 4 1 2 

1 - d + j - 26 + 2i(b 1 + d - j - 26 + 2iuj I -2b 



Ar 2 



„f>+i 



Ar z 



(1- Ar 2 ) 2 + 
'(1-Ar 2 )^. 



where j = 1, 3, 5, for tensorial, vectorial and scalar type gravitational perturbations, respectively. The fact 
that j is always an odd number will have significant consequences in the final result, as we shall see shortly. 
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Next, one needs to impose quasinormal mode boundary conditions. At the origin, there is no horizon 
whatsoever, and one demands for a regular solution. This immediately sets C2 = 0. We are thus left 
with the first term, and the quasinormal mode boundary condition of out-going plane waves near the 
cosmological horizon. At the cosmological horizon, r ~ Rq, one has (l — Ar 2 ) — > and Ar 2 — > 1, while the 
hypergeometric function is particularly simple when its argument goes to zero, i.e., lim z ^o 2-F1 [a, /?, 7I z] = 
1. Therefore, in order to conveniently analyze our solution near r ~ Rc, it is helpful to make use of the 
following transformation law for hypergeometric functions 



(1 - Ar 2 )^ 2 F l 



3- d + j + 2b + 2iQj 3 + d - j + 26 + 2iw 3 + 26 
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(1- Ar 2 )"^r[i^] r[ 
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4 



r 



3-d+j+2b-2iw 
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4 ' 4^ : 

3 - d + j + 26 + 2ub 3 + d - j + 26 + 2iu 



1 — i&i 



1 + iuj 



1- Ar z 



1 - Ar" 



+ 



Near the cosmological horizon (l — Ar 2 ) — > and lim^o 2F1 [oj/SjTI z \ = 1> so that making use of this 
identity the hypergeometric piece is no longer relevant and one can concentrate only on the remaining 
factors. At the cosmological horizon the tortoise is x — * +00 and one can write 



(l-Ar 2 )^ 



(1-Ar 2 ) + ^ 

Out-going waves at infinity (waves that behave as e~ luJX ), satisfying dS quasinormal mode boundary 

conditions, must thus have the form $(r) ~ (l — Ar 2 )* 2 . In other words, one can ensure quasinormal 
mode boundary conditions by requiring that the in-going wave term is zero and that the out-going term 
is non-vanishing. One is thus lead to the conditions 




iuix 



— iujx 
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~ 3-d+j+2b-2iu~ 




4 




4 



+ 0. (4.7) 



The first of the conditions in (4.7) implies that one of the Gamma functions in the denominator must have 
a pole. This is actually not enough: if one of the Gamma functions in the denominator has pole, then so 
will r [iu] have a pole (of the same degree), and so they will cancel each other out. One must thus require 
both Gamma functions in the denominator to have a pole. This will guarantee the vanishing equality. 
Recalling that 26 = d + 2i — 4, it is simple to see that one will have the required poles whenever 
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2icu + 21 + j - 1 = -An, 

2iuj + 2d + 2£ - j - 1 = -4m, 

with both n, m G N. On the one hand both these equations imply that icD will be a negative integer (recall 
that j is an odd number) so that the Gamma function in the numerator will have a pole, as we have stated 
previously. On the other hand, both these equations imply that none of the Gamma functions on the 
second condition in (|4.7|) can have poles in such a way that that expression could vanish. The requirement 
of out-going boundary conditions is thus satisfied. We still have to deal with the fact that there are two 
conditions for the one unknown, the frequency. These equations are compatible, and yield a solution for 
the quasinormal frequency, if and only if 



d - j = 2 (n - m) e 2Z. 

Because j = 1,3,5, d — j can only be even if d is odd. One is led to the conclusion that the pure dS 
spacetime will have quasinormal modes only in odd dimensions. In this case, the quasinormal frequencies 
for tensor (j = 1), vector (j = 3) and scalar (j = 5) type perturbations are given by 



Vx(^2n + e+^^j, n = 0,1,2,.... 



to = i 

One might worry about this result for odd spacetime dimension, as in this case the solution to the wave 
equation is not given by the general expression above, involving two hypergeometric functions. In fact, the 
two corresponding solutions became dependent 5 , and the second hypergeometric function must be replaced 
by a Meijer G-function in order to generate the space of solutions. However, the Meijer G-function has 
the same singular behavior at the origin as the hypergeometric function it is replacing, and hence it must 
be discarded. The remaining hypergeometric function is exactly the same as we had above and thus our 
reasoning remains unchanged, leading to the final result above. Observe that this result is certainly in 



agreement with the recent work in [32], where it was shown that four dimensional dS spacetime has no 
quasinormal modes. It would be very interesting to obtain numerical evidence in support of our results, 
both in this dS case and our previous study of AdS. 



5. Applications to Quantum Gravity 

The quest for asymptotic quasinormal frequencies was subject to increased momentum after the works 
of Hod [O], showing how these frequencies could be used to infer on the quantum of area at the black 



hole event horizon, and of Dreyer [12], showing how the Barbero-Immirzi parameter |73| , 74, [7f|] in loop 
quantum gravity could be fixed starting from these frequencies. Having computed quasinormal frequencies 
in a wide variety of situations, we would now like to understand whether the ideas of p| |l2| are only valid 
for Schwarzschild four dimensional black holes, or if they have a broader scope. Let us point out that it is 
the scalar type perturbations of the gravitational mode which seem to be the most relevant for applications 
in quantum gravity, as it is this type of perturbations which corresponds to changes in the geometry of 
the black hole event horizon |1C|| . As we have also shown that asymptotic quasinormal frequencies for 
scalar type perturbations always coincide with the ones for tensor type, this will actually not be an issue. 
We will next swiftly review these different approaches to the study of quantum gravity from asymptotic 
quasinormal frequencies and see what we c an further learn by making use of our own results. 
We thank Bin Wang for pointing this out to us. 
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5.1 The Quantum of Area from Quasinormal Frequencies 

It has been long advocated that black holes should have a discrete area spectrum |f76| , 77] and that one 
could expect, on general grounds, horizon area to be quantized as A n = f(n), with n G N. The authors of 



1 76, 77] actually make the case for a linear function, with the proposed non-extremal black hole discrete 



spectrum for area being of the harmonic oscillator type (in Gd = c = h = 1 units) 

A n = ntt, n £ N, 

where f2 was an undetermined number. Observe that [76] if the energy at level n is degenerate, with 
multiplicity g(n), one can further identify the black hole entropy at this level with log g(n). So, if one 
makes use of the Bekenstein-Hawking (BH) entropy, S = 4, one can immediately obtain 

g(n) = 5(1) exp j (n - 1) , 

where we have appropriately normalized g(n). Because g(n) is an integer, this further constraints = 
4 log k with k = 2, 3, 4, . . ., The choice of k is, however, subject to some debate. 

Along these lines, JLT| ] made a pertinent observation which allowed for a first principles determination 
of the above value for k (and, in some ways, for an independent check of the above arguments concern- 
ing horizon area quantization), at least in the Schwarzschild case. What [11] did was to use asymptotic 



quasinormal frequencies in order to determine the unknown $7. We have seen that these frequencies are 
located in the complex plan and that the real part of the quasinormal frequencies approaches a constant 
value as the mode n is increased. Making use of Bohr's correspondence principle which states that "at 
large quantum numbers, the [quantum] transition frequencies should be equal to the classical oscillation 
frequencies'", one thus expects that the asymptotic behavior of the quasinormal frequencies may contain 
relevant information about the quantum black hole (the problem with this approach to black hole quan- 
tum mechanics is, of course, to correctly identify the physically relevant oscillation frequencies). For the 
Schwarzschild black hole, the asymptotic behavior for the highly damped frequencies is |l3|, [l4|, |35| 

^ = log3 + 2™ (n + l ) +0 



for any type of perturbation and in any dimension d, with Tjj the Hawking temperature at the Schwarzschild 
horizon. Clearly, Re(w) depends only on the black hole parameters. One thus interprets the highly damped 
quasinormal frequencies (specifically, lim n ^ +00 Me(u; n )) as a characteristic of the black hole itself. 



In order to remain fully general, we now proceed with a slight deviation from the arguments in [11]. 
Let us consider the first law of black hole mechanics 0], which states that a static black hole, with mass 
M, charge Q, surface gravity k, and electric surface potential when perturbed so that M — > M + dM 
and Q — > Q + dQ, reacts as 

k 

dM = — dA+ $\ H dQ, 



basically stating "energy = heat + work". If one makes use of Bohr's principle as suggested in [11], 
dM = Re(a; 0O ), and further understanding that dQ = 0, it immediately follows 

dA = -liRe^oo) = 4 log 3. 

J-H 
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All together, via a quasinormal mode calculation the unknown number Q in the horizon quantum mechan- 
ical area spectrum formula for a given black hole has now been pinpointed to be = 4 log 3. Attempts to 
extend these ideas to other black holes were first performed in |78|, |7s!J but, as we have computed in the 
present work, the asymptotic frequencies used in those papers turned out not to be the correct ones to use. 

When moving to the general case (which is the main focus of this paper), some questions concerning 
the previous procedure come to mind. Let us first address the question of the value of dQ, when in the 
Einstein-Maxwell situation. Recall that we are computing perturbations of both the gravitational and 
electromagnetic modes at the electrovacuum, and that these perturbations are themselves uncharged. So, 
any application of Bohr's principle we may make will always have to satisfy dQ = 0. In this case, for an 
asymptotic quasinormal frequency to be of any relevance in the determination of the area spectrum as 
described above, it better be of the form Me(u; 00 ) = Tjj logx, with x G N. While the proportionality to 
the temperature is somewhat generic from dimensional analysis, obtaining the term in log a; is a rather 
special event. Another source of confusion are situations with more than one horizon. Again we recall that 
one is always working with massless, uncharged, perturbations of the electrovacuum, in which case the 
application of the first law as above still holds. Indeed, generically, if one has two horizons — name them 
1 and 2 — then k\dA\ = &2GL42 and it is irrelevant which horizon we choose to concentrate on. With these 
points cleared and out of the way, let us see what happens in the different situations. 

Actually not much happens, besides the simple Schwarzschild situation which we have just reviewed. 
Indeed, for RN, Schwarzschild dS, or RN dS, there is no algebraic solution for the asymptotic quasinormal 
frequencies and thus no "logarithm" to use in the arguments above. For Schwarzschild AdS and RN AdS, 
there is an algebraic solution for the asymptotic frequencies (depending on the parameter xq), but this 
solution does not have the required logarithmic structure to be used in the arguments above. Finally, in 
the extremal RN situation, there is both a solution and the required logarithmic structure. Unfortunately, 
the argument of the logarithm is only an integer in dimensions d = 4 and d = 5 where it takes the values 
1 and 0, respectively, thus again invalidating the arguments above. Thus, it does not seem likely that the 
arguments in |jll| will have application away from the realm of the Schwarzschild spacetime. 



5.2 The Barbero Immirzi Parameter in Loop Quantum Gravity 

The entropy of a large class of four dimensional 6 non-extremal black holes has been computed within the 
formalism of loop quantum gravity 
entropy with area, S cx A 
Barbero-Immirzi parameter, 7 [73, 74 



81]. While this calculation reproduces the expected scaling of 
it fails to reproduce the numerical factor of \ up to an arbitrary factor, the 
]. This parameter appears in the definition of the [real] phase 



space variables and is thus present in the quantization, fixing the spectrum of the area operator. Because 
it is arbitrary (positive real), it can be adjusted so that the entropy calculation is correct, but then this 
result is no longer a prediction of loop quantum gravity. It has been suggested that one should regard 
this parameter as an analogue of the 9 parameter in QCD Q|. It has also been suggested that one could 
use quasinormal mode considerations in order to fix this parameter semi-classically [12]. Let us see how 



this parameter comes about and how it can be fixed for the four dimensional Schwarzschild black hole 
using quasinormal frequencies. In order to remain fully general in what respects the spectrum of the area 
operator ^2} |83l (for the upcoming analysis of our results), we proceed with a slight deviation from 
the arguments in [|12|1 (including at the same time the supersymmetric results of |$4| ) . 

In loop quantum gravity the degrees of freedom associated to a four dimensional quantum black hole 
are described by a Chern-Simons gauge theory living at the event horizon [80, pi}] . One constructs surface 



Let us stress that this subsection is strictly four dimensional. 
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states with the aid of spin networks such that when an edge of the spin network intersects and punctures 
a surface it associates to it a quantum of area. In order to compute the black hole entropy one starts by 
computing the dimension of the boundary Hilbert space which describes the surface states, dimTCs, and 



then simply writes S = log dim TC^. For a large number of punctures, Chern-Simons theory yields [8C] 

N 

dimH E = II(2j p + l), 

P =l 

where there are iV punctures at the event horizon and where j p labels the SU(2) representation associated 



to the spin network edge piercing at puncture p (see [81] for further details). A similar result holds for 



J\f = 1 supersymmetric spin networks, where this time around J85| ] 

N 

dimWs = JJ(4j p + l), 

P =i 

and j p now labels representations of Osp(l\2). For large horizon area, one further expects the counting of 
surface states to be dominated by states where all the spin network edges at the horizon punctures carry 



the same [fundamental] representation j [8C, 81 1. In this case, for bosonic spin networks 



S = iVlog (2j + 1) , 
while for M = 1 supersymmetric spin networks 

S = iVlog (4j + l). 

Clearly, one still needs to estimate the number of horizon punctures, N, and to determine the fundamental 
representation, j. This can be accomplished via use of the area operator. 

In loop quantum gravity the area operator is quantized. However, it happens that in the classical 



geometrical theory there is more than one operator describing area [32|. Upon quantization, this fact leads 
to two different spectra of eigenvalues: there is a standard area operator [30, 31 1 which leads to a quantum 
of area — for a surface intersected by one edge of a spin network with label j — with eigenvalues (again, we 
use units where all physical constants are set to unity, so that tp = 1) 



A(j) = 87r 7 Vj(j + l), (5-1) 

but there is also a different area operator, introduced in ]82] ], which leads to the alternative eigenvalue 
spectrum 

A(j) = 87T7j. (5.2) 

In the expressions above 7, 7 is the Barbero-Immirzi parameter, finally making its appearance, and the 
total area of the black hole event horizon is given by a sum over all the punctures, A = Y2p=i A(j p )- While 
it is operator (|5.1| ) which is used in most applications, there is no a priori reason to rule out operator fl5.2|) . 
Yet another area operator appears when dealing with supersymmetric spin networks, which has eigenvalue 
spectrum given by Q 




A(j) = 8iryJj + - . (5.3) 
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Here 7 is the supersymmetric Barbero-Immirzi parameter and j is now understood to label Osp(l\2) 
representations instead of SU(2) ones. To keep all possibilities together, we shall proceed with the following 
expression for the area spectrum 



A n (j) = 8Tr<y n ^j (j + n), 

where notation should be clear. When n = 0, 5 and 1, one obtains expressions ( |5.2| ), (|5.3|) and ( |5.lD , 
respectively. Going back to our calculation of the entropy, recall that one expects the counting to be 
dominated by states where all spin network edges carry the same fundamental representation j. Thus, for 
a black hole with total area A and quantum of area A(J), the number of punctures N can be computed 
simply as N = ^jy ; where we still do not know j. The entropy for bosonic and supersymmetric spin 
networks then follows as 

_ log (2j + 1) A - = log (4j + 1) A 

' 4 ' 



2vr 7 „ V / iO' + n ) 4 2vT7^/i (j + |) 



The standard lore ^] now states that the physically fundamental representation j should be the 
mathematically fundamental one, and so j = \. In this case 

s _ log2 A g_ log3 A 
7 n vrV2n + 1 4' 77^ 4' 

and one can obtain S = ^ with an adequate choice of the (super) Barbero-Immirzi parameter 7. On the 
other hand, one would much rather go without the final arbitrary choice of 7 (even though one should 
point out that the very same choice of Barbero-Immirzi parameter holds for the large class of black holes 
studied in [80, plJl , including charged black holes and cosmological situations). 

An alternative path, which we wish to study, was presented in [|l^] (and extended to supersymmetric 
spin networks in [84]). The point of departure from the previous reasoning lies in the computation of the 
number of punctures, N. Dreyer [12] uses instead Hod's |ll| result for the quantum of area, 



4 

dA = — Me^), 

and from dA = A(j) (with an unknown fundamental representation j) it follows a direct expression for the 
Barbero-Immirzi parameter 



In 



{j + n) 

At the same time, from = -A, it also follows a direct and independent expression for the entropy, 

_ r g log(2j + l) A d T H log(4j + l) A 
o = =— ; : and b = — -. 

Me(woo) 4 Me(woo) 4 

The interest of these expressions is best appreciated in the Schwarzschild case, where they become 

S=^ 2j + 1 ^± and S= l °^ + 1) ±. 
log 3 4 log 3 4 



The fundamental representation for the bosonic spin network must then be chosen as j = 1 [12] while the 
one for the supersymmetric spin network must be j = | |84]] . This also leads to an independent result for 
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the Barbero-Immirzi parameter. These calculations have initiated a debate on the correct choice of gauge 
group at the spin network edges, but we do not wish to proceed into that. Instead, we just note that if it 
is not the case that Me(u; 00 ) = Tjjlogx, with i£N, then the previous argument will fail and one is back 
at the initial arbitrary choice for 7. 

It so happens that the previous calculations apply to a broad range of horizons, including cosmological 
horizons []80| , 81]. This is an immediate problem, as we have shown that, with the exception of the 
Schwarzschild case, one cannot find solutions with the required logarithmic structure and with x in the 



integers. So, it does not seem likely that the arguments of [12] will have application away from the realm 
of the Schwarzschild spacetime and, generically, one will thus not be able to fix the arbitrary Barbero- 
Immirzi parameter using asymptotic quasinormal frequencies. Yet another problem has recently emerged 
directly within loop quantum gravity [^, 87]. These authors have shown that the initial calculation of 
black hole entropy in loop quantum gravity |8(| ^] was actually incorrect: there was a miscounting of 
the microscopic states contributing to the entropy. In fact, the lore (mentioned above) that the physically 
fundamental representation j should be the mathematically fundamental one, of j = i, is actually not 
true, as states labeled by higher spins also contribute for the entropy calculation [36, 37 1. This has clear 



implications in the derivation of the BH black hole entropy within loop quantum gravity; of particular 
interest to our discussion above is the fact that it affects the derivation of the Barbero-Immirzi parameter 
for Schwarzschild black holes. Indeed, instead of having 

log 2 

for the Barbero-Immirzi parameter, one actually has that 7 obeys a complicated functional relationship 
and is a number which lies in the range p~ 



log 2 log 3 
< 71 < • 

However, it does not seem that this number can be written as a logarithm [^, |87|] which would even 
invalidate an attempt to match it against the Schwarzschild asymptotic quasinormal frequencies. There 
has been some attempts to try to keep the Schwarzschild quasinormal relation alive in |]BS| , |8^ , |90| , but at 
the moment the dust is far from settled. It thus seems that asymptotic quasinormal frequencies will not be 
of much help in loop quantum gravity. Their future role in quantum gravity is at this moment unknown. 



6. Conclusions and Future Directions 



In this paper we have learned how to analytically compute asymptotic quasinormal frequencies in a wide 
variety of black hole spacetimes. This made use of a monodromy technique which was first introduced in 
p4| for asymptotically flat spacetimes, and later extended to non-asymptotically flat spacetimes in fHf . 
Such a complete classification as the one we have performed in this paper certainly proves the analytical 
power of the monodromy technique and opens way to many new results in the calculation of quasinormal 
modes. The directions in which one can proceed are multiple. Perhaps the first question to address should 
be the Kerr black hole spacetime, as one would like to know whether the monodromy technique we have 
developed at length in this work may be extended to non-static spacetimes. Another question to address, 
and closer to the discussion in this paper, would be to compute asymptotic quasinormal frequencies of 
the many extremal situations which arise in our setting and that we have not considered. To do this, one 
should bear in mind our example of the extremal RN black hole, i.e., the extremal solution calculation may 
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need to be considered from scratch, independently of the non-extremal one. Also, the full significance of 
the continuous spectrum for scalar type gravitational perturbations in dimension d = 5, for asymptotically 
AdS spacetimes, which we have found in this paper, should be further explored. 

On what concerns quantum gravity our results are not very encouraging, at least on a first reading. 
Indeed, they seem to indicate that the results in |ll], [l^] are nothing but numerical coincidences. In spite of 
this, we believe that the calculation we have performed in this paper is opening the door for the calculation 
of asymptotic greybody factors, and these will certainly have a key role to play on the road to quantum 
gravity. Thus, one can view this paper as the first step towards a complete classification of greybody factors 
for black hole spacetimes. These greybody factors will describe scattering at high imaginary frequencies 
and may yet have a role to play on the identification of the dual CFT which microscopically describe the 
several black holes we have considered. The greybody factors may also provide for additional information 
on the statistical mechanical behavior of these black hole spacetimes. 

In spite of the fact that there is an intrinsic general relativistic interest in the computation and 
classification of asymptotic quasinormal frequencies, the goal of late on this line of problems has been 
to find new applications on what concerns quantum gravity. It thus seems like a natural step to take 
this calculation and apply it to string theoretic black holes. Here, one is immediately faced with a wide 
range of open problems, starting from higher derivative stringy corrections to the Einstein equations, 
proceeding towards many D-brane constructions of higher dimensional black holes, and passing through 
the understanding of the asymptotic quasinormal frequencies of the Schwarzschild AdS black hole in terms 
of the dual finite temperature gauge theory. 

If one is to make the step towards string theory, we should finish by mentioning that there is another 



line of research which deals with the so-called non-quasinormal frequencies [91] and the study of horizon 



CFT |92], |93|, |9J, [9f|]. Here, it seems that it is frequencies closely related to the quasinormal frequencies 
that are actually able to determine the correct quantization of the Virasoro operators associated to the 
horizon CFT, at least in a dilute gas regime, thus yielding the correct BH entropy via use of the Cardy 
formula 91, 95 1. Let us thus finish by reviewing one last application of quasinormal modes to quantum 
gravity, in this context of Cardy's formula for the asymptotic growth of states in CFT To put what 



we wish to do in perspective, let us start by reviewing a remark by Strominger [92], itself based on an 
earlier observation of Brown and Henneaux W^ - For the moment, we concentrate on three dimensional 
quantum gravity. What pTj showed is that any theory of quantum gravity in AdS3 must be described by 



a CFT with central charge (as before, we use units where all physical constants are set to unity) 



2\/|AT 



32[ takes this remarkable observation one step further by making use of Cardy's formula for the asymptotic 
growth of states in this particular CFT, thus being able to [microscopically] compute black hole entropy. 
Because the results in p7| guarantee that black hole states belong to the Hilbert space of the given CFT, 



one may concentrate in the BTZ case [98] and show that mass M and angular momentum J can be related 



to the Virasoro operators, so as to yield the eigenvalues of the Virasoro operators themselves |K| , 



M+J\\\J , - M-J\\\J 

Ln = t — — and Ln = } — — , 

2x/\X\ 2^r\X\ 



where one should take into account that these semi-classical considerations hold only when the central 
charge is c S> 1. Now Cardy's formula for the asymptotic growth of the number of states, p(Lq, Ln), of a 
given CFT states that (when Lq\0) = 0, |0) being the ground state) 
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logp(L 0) Io)~27r^+27rJ^, 
and the entropy can be computed from 5 = log p(Lq, Lq). Direct application to the BTZ black hole yields 

S-- 

in precise agreement with the BH result for this black hole. This line of reasoning is actually quite 
generic (see [92] for further details) but applies only to black holes whose near-horizon geometry is three 




dimensional. What one would wish for would be similar arguments valid in any dimension. 

Such an approach has been suggested by Carlip |9^, |94], |95|| , in an attempt to generalize the previous 
arguments to ci-dimensional black holes. In |94| one concentrates on the symmetries at the horizon of 
an arbitrary d-dimensional black hole. Using the ADM formalism one studies the algebra of constraints of 
general relativity in the presence of a boundary and finds that the presence of such boundary actually leads 
to a central extension of the algebra of spacetime diffeomorphisms. By then concentrating in a subalgebra 
isomorphic to Diff (S 1 ) , one obtains a Virasoro algebra where one can then apply Cardy's formula as before. 
The algebra in question, that of black hole horizon surface deformations, is the algebra of deformations of 
the (t, r) plane which leave the horizon fixed. In this "horizon CFT" holomorphic and anti-holomorphic 
functions correspond to functions of t ± x, where x is the tortoise coordinate J95|. Working within classical 



general relativity alone, |93|, |)4[ showed that the Virasoro algebra in question should have (with A the 
horizon area) 

3A A T A 

c = — and L = —-, 
2tt loir 

with the right moving modes having vanishing central charge and Virasoro operator. Applying Cardy's 
formula one then immediately obtains S = 4. It is important to realize one of the main messages in 



p3| , 94, 63] : that by making use of general relativity's classical symmetries alone, one may be able to 
determine the black hole entropy, independently of any microscopic analysis. |}l]] now proposes a semi- 
classical way to fix the eigenvalues of the Virasoro operator Lq with the help of quasinormal frequencies. 

Let us review where the authors centered upon the concept of non-quasinormal frequencies. 
These non-quasinormal frequencies are defined using monodromies at inner and outer horizons alone, 
rather than involving any condition at infinity. Recall that quasinormal mode boundary conditions are 
defined in terms of the behavior of perturbations at the horizon and at infinity. It so happens that for 
the three dimensional BTZ black hole these conditions can be recast in terms of monodromy conditions 
at the inner and outer event horizons alone [|9l] (further insight into the origin of these non-quasinormal 
frequencies has appeared recently in |99|]). One is thus naturally lead to a definition of non-quasinormal 
frequencies, for d-dimensional black holes, only in terms of these monodromies, via 9JT (R^) 931 {Rjj) = 1- 
This is accomplished by choosing ingoing boundary conditions at the outer horizon and a general solution 
at the inner horizon. This monodromy requirement for solutions leads to two non-quasinormal frequencies, 



uj±. The idea of [91] is then to relate these frequencies, oj±, to the quantization of the Virasoro operators, 
Lq and Lq, via a natural application of the correspondence principle to this situation. Surprisingly, in 
a dilute gas regime, these non-quasinormal frequencies lead to the correct quantization of the Virasoro 



operators as well as the correct quantizations of energy and angular momentum |91[]. It should be very 
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interesting to find an extension of those results to a broader situation, such as the one in this paper. It 
may be that quasinormal modes hold more information about quantum gravity than expected at first! 

To finish, let us stress that the study of asymptotic quasinormal modes is an active area of present 
research and has also led to many other applications besides the ones we have studied in this paper. Some 
recent references one should keep in mind include jlOC 



101, 102, 103, 104, 105, 106, 107, 108, 109]. 
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A. List of Conventions for Black Hole Spacetimes 



In this paper we consider spherically symmetric, static, black holes in <i-dimensional spacetime, with mass 
M, charge Q and background cosmological constant A. This set (M, Q, A) of parameters thus describes 
the black hole. In this appendix we review and set notation on basic facts for d-dimensional (M, Q, A) 
black holes. The Einstein-Maxwell equations describing our backgrounds are the well known 



7A 



V A F Xfl = 0, 

where Gd is the ci-dimensional Newton constant, and the energy-momentum tensor for the electromagnetic 
field strength F = dA is given by 

Tfiv = F^ F u \ ~^9fj,vF(xpF . 
The background spacetime metric g, solution to the previous Einstein-Maxwell system, is given by 

g = -f( r ) dt®dt + fir)" 1 dr®dr + r 2 dVL 2 d _ 2 , (A.l) 
where the function /(r) is the well known (we only consider d > 3) 



/W = l-^3+i"Ar 2 . (A.2) 



In terms of the background parameters fi, 9 and A in the metric, one can further compute the ADM mass, 
M, of the black hole as 

where A n is the area of an unit n-sphere, 

n+l 

Likewise one can obtain the charge, Q, of the black hole as 

and the value of the cosmological constant, A, as 

A=-(d-l)(d-2)A. (A.6) 

These parameters /x, and A, are thus a convenient way of keeping track of the black hole defining 
parameters, M, Q and A, in the d-dimensional case [||, [l(], |1 XCf ] . 

Let us denote — for the moment — the black hole background spacetime metric by g^, and let us 
perturb it with some disturbance h^ u (where <C 1 in some appropriate sense). Quantities with a hat 
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on top will here arise from the background metric, while un-hatted quantities will arise from the linear 
perturbation. The metric to consider is thus g^ u = g^ u + h^ u . While the background g^ u satisfies Einstein's 
equations (or the Einstein-Maxwell equations, depending on the conditions one is considering), one still 
needs to compute the linearized equations of motion describing h^ u . The perturbed components of the 
Levi-Civita connection are 

where the variation of the connection is actually a tensor, 



1 

2 

From this result one immediately has the variation of the Ricci tensor 



and, for vacuum solutions, Einstein's equations reduce to dR^y = (in the Einstein-Maxwell case there 
will be other terms). These are linear differential equations for the perturbation, huv, which can still be 
further simplified using the spherical symmetry of the background g^ u . For (f-dimensional spacetimes these 
equations were recently studied by Ishibashi and Kodama in g, ^, |l(]] and we shall return to their results 
in appendix B. 

The black hole backgrounds defined by the function f(r) can, generically, have both event horizons 
and cosmological horizons, which shall be denoted by Rh and Rc, respectively. They will be simply given 
by the zeros of f(r), i.e., the locations in which f(r) changes sign. Also, the parameters /x, 9 and A in the 
background metric cannot generically be taken as arbitrary, but are constrained if one is not to consider 
naked singularities. The general analysis can be quite intricate [10], and we shall in here mainly list the 
results of these authors, referring the reader to for the details. The several possible cases are: 

The Schwarzschild Solution: In here f(r) = 1 — -|=j> and the event horizon, determined by f(r) = 0, is 
located at 

R H = (2/x)^ , 

so that f(r) > for r > Rh- The parameter fi can take any positive value and there is no cosmological 
horizon. 



The RN Solution: In here f(r) = 1 — id-a , and the event horizons, determined by f{r) = 0, correspond 
to 



i 

d-3 



where one is usually interested in the outer horizon, R H . In this case f(r) > for r > R H . This result 
also leads to the constraint 6 2 < /j, 2 ; the extremal solution corresponding to 2 = fi 2 in which case there 
is only one solution for the horizon, 



i 



R H = R H = ^ d - i - 



Again, there is no cosmological horizon. 
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Solutions with a Cosmological Constant: For solutions with cosmological constant alone, f(r) = 1 — Ar 2 , 
there is only a cosmological horizon, corresponding to 

*° = js: 

It follows that f(r) > for r < Re- This clearly constraints A > 0, so that the AdS solution has no 
cosmological horizon — it will only exist for the dS solution. Because there is no black hole, there is no 
event horizon. 

The Schwarzschild AdS Solution: In here /(r) = 1 — — Ar 2 and A < 0. There is an event horizon 
implicitly defined by f(r) = 0, i.e., by the unique real and positive root of 

|A|r d_1 + r d - 3 -2/u = 0. 

There is no cosmological horizon. While fi can take any positive value, A can take any negative value. 
The RN AdS Solution: In here f(r) = 1 - + - 

.2d- 6 — Ar 2 and A < 0. There are event horizons, implicitly 

defined by f(r) = 0, i.e., by 

\ X \r 2d -' l + r 2d - e -2nr d - 3 + 2 = 0, 

and no cosmological horizon. This equation generically has two real positive zeros (corresponding to an 
inner and an outer horizon, as is familiar in the RN class of solutions) and /(r) > for r > (as usual, 
one is interested in the outer horizon). The background parameters are constrained to obey 6 2 < [i 2 and 
the cosmological constant is negative but bounded from below flC 



12 



> A > -2^=3 



d _ 3 J(d-l) 2 fJ , 2 -4(d-2)9 2 -(d-3) t M 



(d-l)fi- ^J(d-l) 2 n 2 -4 (d-2) 2 ^ 



d-l - 

d-3 



The condition 2 = [i 2 cannot be achieved in here as that would imply A = and the solution would no 
longer be AdS. Rather, the extremal solution for this case is when the cosmological constant reaches its 
minimal value 



2 d-3 J(d-l) 2 ^-4(d-2)6 2 -(d-3)fj, 

_2d=3 ¥ 

d-2 / i \ fci 

l(d-l)n-^/(d-l) 2 t j, 2 -4(d-2)e 2 ) 



From [1C], one finds that at extremality of the AdS cosmological constant, inner and outer horizons will 
coincide. This point of coincidence also corresponds to a minimal value for the [outer] horizon radius, and 
is located at 



R H = R H 



2 {d-2) e 2 



i 

d-3 



(d - 1) n + J(d - l) 2 /x 2 - 4 (d - 2) 2 



As the AdS cosmological constant grows from extremality to zero, the [outer] radius of the RN AdS black 
hole grows from this minimum to its maximum, which corresponds to the pure RN solution. 
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The Schwarzschild dS Solution: In here f(r) = 1 — — Ar 2 and A > 0. There is both an event horizon 
and a cosmological horizon, implicitly defined by f(r) = 0, i.e., by 

\ r d-i _ r d-3 + 2tl = 0. 

This equation generically has two real and positive zeros (the event horizon and the cosmological horizon) 
and f(r) > for Rh < r < Re- While /x can take any positive value, A is constrained to be positive but 
smaller than 

d-3 1 

< A < 



12 



d - 1 [(d-l)/x]3=3 

Because of this constraint it follows that there is an extremal solution in this case, precisely corresponding 
to the case where A attains its maximum value, 

d-3 1 



2 

d-3 



d - l [{d-l)^ 

At this extremal point, the event horizon and the cosmological horizon will actually coincide. The solution 
at this point can be found to be 



1 

-3 



R H = R c = [(d-l)n}*-- 
Observe that this is a quite particular situation, as there is now no range of r where f(r) > 0. 

The RN dS Solution: In here f(r) = 1 - -3^ + - Ar 2 and A > 0. There are both event horizons and 
a cosmological horizon, implicitly defined by f(r) = 0, i.e., by 

Ar 2d-4 _ r 2d-6 + 2fir d-3 _ e 2 = Q 

This equation generically has three positive real roots (inner and outer event horizons plus the cosmological 
horizon) and f(r) > for R^j < r < Re- The constraints on the background parameters are a bit more 
intricate now 110]: one can choose 2 < fi 2 in which case A is constrained to be positive but smaller than 



2 d _ 3 (d-3)/x + V(d-l)V-4(d-2)^ 



0<x ^a-2, , : , a 



(d- l)fi+ yj(d- 1)V 2 -4(d-2)6>^ 

On this range of parameters the choice of extremality corresponds to making A extremal (the standard 
cosmological constant extremality condition — observe that the standard RN extremality condition, of 2 = 
[i 2 , is in general not extremal for A / 0). One thus requires 



j_ d - 3 (d-3)n+ J(d-l)V"4(d-2)0 2 

3 d 

(d - 1) n + y/(d - 1) 2 n 2 - 4 (d - 2) 6> 2 

Instead of three positive real roots one now only finds two positive roots, corresponding to the dS extremal 
condition in which the [outer] event horizon and the cosmological horizon coincide. The solution at this 
point is [l0|| 
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(d - 1) n + J (d - if /x 2 - 4 (d - 2) 6 2 



i 

d-3 



As we have said previously, in this RN dS case the constraints on the background parameters are more 



intricate. Thus, besides the choice we have just studied, one can also study the case where [1C 



„2 / 02 < 1 (d-1) 2 2 
4 ~~d~^2~ ' 

in which case A is constrained to be a positive value bounded both from above and below as 

A < A < A+, 

where 



\± = 2^-3 



d 



3 (d-3)n±J(d-l) 2 fi 2 -4:(d-2)6 2 



d-2 



(d-l)n± yj(d- 1)V -4(d-2)0 2 ^ 



d-1 ' 

d-3 



In this case, the equation for the horizon radius will generically have three real zeros (again corresponding to 
both an inner and an outer event horizons plus the cosmological horizon) where f(r) > in R^j < r < Rc- 
There now seem to be several choices for extremality. The first choice we will make, exploring extremality 
conditions on 6, is to make the extremal choice of 



= i (d-iy 

4 d-2 ' 



This situation is actually interesting as in this case 



A_ = 2^ < d - 3 ' 2 



(d-2) (d- 1)3=5 fj,d-s ' 

so that this automatically also makes the cosmological constant extremal (at the above value). In this 
case, one finds a single real and positive root, corresponding to a situation where inner event horizon, 
outer event horizon and cosmological horizon, all coincide at the very same point. This is the dS type of 
extremal condition, and the point of coincidence of the three possible horizons is at 

r h = r h = R c= (~ 

The other possible extremal choices one can make leave 6 in its parameter range, and saturate A at either 
A_ or A + . These extremal conditions on A will make two of the three real positive solutions for the horizon 
radius coincide. If one saturates A = A_, one will only find two real positive roots, and this corresponds 
to a RN type of extremality, where inner and outer event horizons coincide. On the other hand, if one 
saturates A = A + , one will again only find two real positive roots, but this will now correspond to a dS 
type of extremality, where the outer event horizon and the cosmological horizon coincide. Observe that 
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this is coherent with the previous picture of full saturation at A + = A_ where all three horizons coincide 
at a point. At the A = A_ extremal point, inner and outer horizons coincide at the point 



R H = R H 



(d-l)n- \/(d-l) 2 fi 2 -4(d-2)e 2 " ' 



2 

Likewise, at the A = A+ extremal point, outer event horizon and cosmological horizon coincide at 



R + H = Rc 



(d - 1) n + J(d - l) 2 /x 2 - 4 (d - 2) e 1 



i 

d-3 



Recall that situations where the [outer] event horizon and the cosmological horizon coincide are very 
particular, as there will be no range of r where f{r) > 0. 

B. Formulae on Ishibashi— Kodama Master Equations 

A set of equations describing the perturbations to background black hole d-dimensional spacetimes were 
derived by Ishibashi and Kodama in ||, ||, [k|. The most general case including (M, Q, A) non-vanishing 
parameters is studied in |l0| and in this appendix we shall schematically review the results in that paper 
for the Schrddinger-like equations describing the quasinormal modes. In d dimensions, the IK master 
equations present a gauge invariant formalism for perturbations to the background spacetime, and are 
of the Schrodinger-like form fl3.2|) . These variables are grouped in three types, according to their tenso- 
rial behavior on the S d_2 sphere: tensor type perturbations, vector type perturbations and scalar type 
perturbations ||, [HJ. Perturbations of the gravitational metric field will include all these three types of 
perturbations, while perturbations of the electromagnetic vector field will only include vector type and 
scalar type perturbations As the Einstein-Maxwell system couples gravitational and electromagnetic 
fields, one expects the master equations to be coupled as well. What [^] shows is that one can always de- 
couple the master equations and these are the relevant gauge invariant equations to be used when studying 
quasinormal modes. Because quasinormal modes are studied as perturbations to the electrovacuum, the 
IK master equations will all be homogeneous. Let us list the IK equations and respective Schrodinger-like 
potentials: 

Tensor Type Perturbations: Only the gravitational metric field displays this type of perturbations. The IK 
master equation is of the Schrodinger-like form (as for the previous scalar wave equation case), where the 
potential to consider is now the following H [HJ 

*w = m ( 1 11 +d - 3) + id - 2) <j; 4) m + ld - l f ' (T) ) . (b.i) 

One immediately realizes that this is precisely the same potential as in the case of the massless, uncharged, 
scalar wave equation. It thus follows at once that quasinormal modes of tensor type perturbations of 
the gravitational metric field will coincide with the quasinormal modes one obtains from the scalar wave 
equation. In this situation there is also a special mode with t = 0, corresponding to a constant [tensorial] 
eigenfunction of the laplacian H 1C]. One thus disregards this mode and takes i = 1, 2, 3, . . .. 



Vector Type Perturbations: Both gravitational metric field and electromagnetic vector field display this type 



of perturbations, and one thus expects the master equations to be coupled. It was shown in llOf] that these 
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equations can actually be decoupled so that the IK master equations for these fields are again of the 
Schrodinger-like form. In the case where Q = there is of course only a gravitational metric field to worry 
about. Let us see both cases in turn. In the uncharged Q = case there is only one IK master equation 
||], with potential 

Vy (r) = fir) ( £ ( e + d -V + (ci-2) (d-4)/(r) _ (d - 1) (d - 2) g 

f £(£ + d-3) (d-2) (d-4)/(r) _ rf"(r) \ 
" Aj V r2 + 4r2 2(d-3)J- [ ' 

In the charged Q ^ case, there are two decoupled IK master equations fllOfl, with potentials 



l(l + rf-3) (d-2)(d-4)/(r) (d-l)/x , (d-2Ye\ A A 
Vv±(r) = /(r)( ^ + — 2 + r 2d-4 ^^ll' ( B - 3 ) 



where the $ + equation represents the electromagnetic mode and the $ equation represents the gravita- 
tional mode. Thus, when 6 = 0, one has that Vy-(r) = Vy{r). In the expression above 



A= J (d - l) 2 (d - 3)V 2 + 2(d - 2) (d-3) (£(£ + d-3) - (d-2) )0 : 



These potentials are clearly different from the previous ones characterizing tensor type perturbations of the 
gravitational metric field and characterizing the scalar wave equation. However, in some cases, they will 
have the same leading singularities at the origin as those other potentials. For vector type perturbations 



there are two special modes ||, 1C]. One is £ = 0, and corresponds to a constant [vectorial] eigenfunction 
of the laplacian. The other special mode has £ = 1, whose solution corresponds to addition of a small 
rotation to the black hole background, and which leads to no gravitational dynamical degrees of freedom. 
One therefore disregards both these modes and takes £ = 2, 3, 4, . . .. 

Scalar Type Perturbations: Again, both gravitational metric field and electromagnetic vector field will dis- 
play this type of perturbations, so that the master equations will be coupled. Also in this case it was 
shown that these equations can be decoupled and that the IK master equations for these fields are of 
Schrodinger-like form. In the Q = case one clearly only needs to be concerned about the gravitational 
metric field. Let us see both cases in turn. In the uncharged Q = case there is only one IK master 
equation [||, with potential 

V5{r) ~ 16r^(r)' (B - 4) 



where 



H {r) = e(e + 4-3)-(i-2) + {d - 1 lii- 2 '» i , 



r 

and 
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U(r) 



2 

4d (d - l) 2 (d - 2) 3 - 24 (d - 1) (d - 2) 2 (d - 4) <j^ + d - 3) - (d - 2) + 



+4 (d - 4) (d - 6) {,£ (t + d - 3) - (d - 2) }' 



Ar 2 + 8(d- l) 2 (d-2)' 



„3d-9 



+ 4(d-l)(d-2) 



4 (2d 2 - lid + 18) (£ + d - 3) - (d - 2) } + (d - 1) (d - 2) (d - 4) (d - 6) 



A* 



-2d-6 



24 (d - 2) 



(d- 6) + 3)- (d- 2) } + (d- 1) (d- 2) (d-4) + d - 3) - (d-2)}^ + 

+16^ (£ + d - 3) - (d - 2) } 3 + 4d (d - 2) (£ + d - 3) - (d - 2) } 2 . 



In the charged Q ^ case, there are two decoupled IK master equations [10], with potentials 



f(r)U±(r) 
64r 2 H±(r] 



(B.5) 



where 



H + (r) = 1 + 



(d — 1) (d — 2) (1 — «) 



2<j^(£ + d-3)-(d-2)} 



„d-3 ' 



flL(r) = I (£ + d- 3) - (d - 2) + \ (d- 1) (d - 2) (1 + 0) , 



with the definition, here and below, that 



4{e(e + d-3)-(d-2)\ 6 2 
1 -I L 1 



(d-ir 



and where 
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17+ (r) 



i d {d-lf(d-2f{l-nf (j^ | ig (d-l)(d-2) 2 (d-4)(l-0) im 



{^(£ + d-3)-(d-2)} 



2 r 2d-6 



£ (£ + d - 3) - (d - 2) r rf - 3 



-16 (d-4) (d-6) 



o , (d-l) 4 (d-2) 3 (3d-8)(l-ft) 3 (l + ft) /x 4 



{£ + d - 3) - (d - 2) } v 



Ad-12 



+ 



,(d- l) 2 (d-2) 2 (l-ft) 2 



|^(£ + d-3)- (d-2)}' 



- (d - 1) (3d - 8) (1 - fi) + 4d 2 - 15d + 12 



A* 



„3a!-9 



+ 



8 (d — 1) (1 — n) 
'£(£ + (1-3) - (d- 2) 



(d - 1) (d - 4) (d - 6) + d - 3) - (d - 2) + (d - 2) 2 ) (1 - n) 
2{^ + d-3)-(d-2)} 



+ 



+7d 3 - 49d 2 + 126d - 120 



„2d-6 



+ 



16 (d - 1) (^3 (d - 2) 2 (d-4) - 4^£ (£ + d - 3) - (d - 2) }) (1 - n) 



I {£ + d - 3) - (d - 2) 



+ 



+32 (d - 4) (3d - 8) 



-^3 + 64{£ {£ + d - 3) - (d - 2) } + 16d (d - 2) , 



Mr) = 



- 4d (d - l) 2 (d - 2) 3 (1 + n) 2 + 48 (d - 1) (d - 2) 2 (d - 4) + d - 3) - (d - 2) } • 



■ (1 + fi) ^3 - 16 (d - 4) (d - 6) (£ + d - 3) - (d - 2) }' 



,2 , (d-l) 4 (d-2) 3 (3d-8 
£(£ + d-3)-(d-2) 



(1 - n) (1 + o) 3 + s (d - 1) 2 (d - 2) 2 (1 + 17) 2 



1 



(d-l)(3d-8)(l-fi) + 



+ (d-2Y 



„3d— 9 



+ 8(d-l)(l + n) 



1 



2 (d - 1) (d - 4) (d - 6) y {£ + d - 3) - (d - 2) + 



+ (d - 2) 2 ^ (1 - 0) + 4 (d - 2) (2d 2 - lid + 18) \t {£ + d - 3) - (d - 2) } + 

^-32{^ + d-3)-(d-2)}- 



+ (d-l)(d-2) z (d-4) (d-6) 



i (d - 1) (^3 (d - 2) 2 (d-4) - 4{£ {£ + d - 3) - (d - 2) }) (1 - il) + 



+3 (d - 2) (d - 6) (£ + d - 3) - (d - 2) } + 3 (d - 1) (d - 2) 2 (d - 4) 

+64h {£ + d - 3) - (d - 2) } 3 + 16d (d - 2) + d - 3) - (d - 2) } 2 . 
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r d— 3 



+ 



In these expressions, the <3? + equation represents the electromagnetic mode and the $ equation represents 
the gravitational mode. Indeed, when 6 = 0, one observes that 



H-(r) 
U(r) 



0=0 



0=0 



H(r) 



4U(r) 



so that V$- (r) = Vs(r) at 9 = as expected. This is simple to check as when 6 = one has fi = 1. These 
scalar potentials are different from the previous ones characterizing tensor and vector type perturbations 
of the gravitational and electromagnetic fields. However, in some cases, they will have the same leading 
singularities at the origin as those other potentials. Also in here there are two special modes [||, 10]. 
The mode with I = corresponds to a constant [scalar] eigenfunction of the laplacian and represents a 
small change in the background parameters, fi, 6 and A. The other mode, with £ = 1, leads to an ill- 
defined master equation involving gauge dependent quantities, and to no gravitational dynamical degrees 
of freedom. One again naturally disregards both these modes and takes I = 2, 3, 4, . . .. 



C. The Tortoise and the Master Equation Potentials 

As explained elsewhere in this paper, the master equation for perturbations takes a Schrodinger-like 
form (|3.2| ) when written in terms of the tortoise coordinate. This coordinate explicitly depends upon 
the background spacetime in consideration and at first sight it would seem useful to obtain closed form 
expressions — relating the tortoise coordinate x with the standard r coordinate appearing in the spacetime 
metric — for each of the black holes we consider in this paper. This we shall do in the following. However, 
such analytic expressions are not needed for the calculation we carry out. As explained in the main text, the 
monodromy method we use requires only information around the region at infinity, the region at the origin, 
and the region at the event horizons (nevertheless we obtain such analytic expressions for completeness). 

Another point which needs to be taken into account concerns the master equation potentials described 
in the previous appendix. As we have just said, the monodromy method only makes use of their values in 
three specific regions: infinity, the origin and the horizons. Their full expressions, as presented earlier, are 
not required except for an exact analysis of the quasinormal modes (rather than an asymptotic analysis as 
the one we perform). Thus, one needs to list the values of the master potentials in the different regions of 
interest, both in terms of the metric coordinate r and the tortoise coordinate x. In this appendix we make 
an exhaustive study of the tortoise coordinate for the spacetimes considered in this paper, alongside with 
a list of the master potential values in the regions of interest for the monodromy calculation. 

Starting with the master equation potentials, it is simple to obtain their value at the black hole horizon: 
they all vanish. Indeed, as can be seen from the previous appendix, for all types of perturbations, tensor, 
vector and scalar, and both for charged and uncharged black holes, all master potentials are given by the 
product of f(r) by a function which is regular at the horizon Rh- Thus, since /(Rh) vanishes, so do all 
master potentials, V(Rh) = 0. Let us next analyze all the cases we address in the main text: 

The Schwarzschild Solution: The horizons are defined via r d ~ 3 — 2fi = with the roots 



Rn 



(2/x)- 



2iri 

exp(^— -n], n = 0, l,...,d-4, 
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so that one finds (d — 3) complex horizons (only one real positive root). The arguments of the horizons 
are those of the (d — 3) roots of unity. Given these roots one can proceed and factorize f(r). Then, all 
one has to do in order to find the tortoise coordinate is an integration of a rational function. Choosing 
x[r = 0] = it follows 

7M =r + S^' og V K) 



RnJ ' 



where k n = ^/'(-Rn) is the surface gravity at the given horizon and the sum goes through all possible 
horizons. Two special regions are of interest for our calculation: near r = one finds 

I r r d-2 
x\r] ~ -— / dr r d ~ 3 = -- 



2/i J 2{d-2)n' 
while near r = oo one finds 



[r]~ J 



dr 1 = r. 



Recalling that quasinormal modes live in the region r > Rjj, one can study the asymptotics of the tortoise 
and find 



lim x[r] = —oo, 

r\R H 

lim x[r] = +oo. 

r— >+oo 

Next we turn to the master equation potentials, as listed in the previous appendix. Using the general 
formulae presented earlier it is simple to obtain the following asymptotics. Near r = one finds 



V y (r) ~ 3(d-2) 
V s (r) ~ -(d-2) 



r 2d-4 4 X 2 4x 2 ' 

2 JJ?_ J_ _ 2 2 -l 

2 A* 2 1 2 -l 

r 2d-4 4 X 2 4 X 2 ' 



while near r = oo one finds 



V T (r) ~ 0, 
Vv(r) ~ 0, 
V s (r) ~ 0. 

The RN Solution: The horizons are defined via r M_6 — 2fir d ~ 3 + 6 2 = with the roots 
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so that one finds (2d — 6) complex horizons (only two real positive roots). On what concerns the angular 
direction in the complex r-plane, the horizons are located at the (d — 3) roots of unity. Given these roots 
one can proceed and factorize f(r), and the tortoise again follows from integration of a rational function. 
Choosing x[r = 0] =0 one obtains 

f dr d ~ 4 1 / r \ d ~ 4 1 / r \ 

*i -j W ) - r+ s w iog i 1 - + s « ios i 1 - • 

where = ^f'(R^) is the surface gravity at the given horizon and the sum goes through all possible 
horizons. Two special regions are of interest for our calculation: near r = one finds 

x[r] ~ 4 / dr ^ = T 1 - 

LJ e 2 J (2d-5)e 2 

while near r = oo one finds 




Recalling that quasinormal modes live in the region r > Rjj, one can study the asymptotics of the tortoise 
and find 



lim x[r] 
r\R+ 

lim x\r] 

r-»+oo 



-oo, 



+oo. 



Note that when addressing the extremal RN solution, the above formulae for the tortoise will change. First 
of all, the horizons will now be defined via r d ~ 3 — [i = (at extremality = 6), with the roots 

2m 



Rn 



1 






exp ( 



d- 3 



n 



n 



0,l,...,d-4, 



so that one finds (d — 3) complex horizons (only one real positive root). The arguments of the horizons 
are those of the (d — 3) roots of unity. Given these roots one factorizes f(r) and the tortoise again follows 
from integration of a rational function. Choosing x[r = 0] =0 one obtains 



x[r]= Iw 



d-A 2 



1 d _ 2 d ~ 4 ( r \ 

^ + (^3FS fl " log l 1 -^J' 



except in d = 4 where one still has to subtract fi to the above expression in order to satisfy x[r = 0] =0. 

Next we turn to the master equation potentials, as listed in the previous appendix. Using the general 
formulae presented earlier it is simple to obtain the following asymptotics. Near r = one finds 



W„ _' (( ,_ 2) ( M - g) * = _ ("-«)(«- 
4 V ;v y r 4d - 10 4(2d-5) 2 

y v± (r) ~ i (d - 2) (5d - 12) ^ = 
v W 4 l M ; r 4d-io 4 (2d-5) 2 x 5 

. , 1 , , w , , 4 (d- 2) (3d- 

F s ±(r T (d-2 (3d -8)^-^ = -- J --^ r 

5 W 4 V M V 4d - 10 4(2d-5) 2 

where jj = js± = 4i=5 anc ^ Jv± = Ijrf = 2 — jy- Near r = oo one finds 



Vy(r) ~ 0, 
F v ±(r) ~ 0, 
V s ±(r) ~ 0. 

The Schwarzschild dS Solution: The horizons are defined via — Ar d_1 + r d ~ 3 — 2/x = 0, and there is no 
general analytic solutions for the roots. In even dimension there is an odd number of roots yielding (d — 1) 
complex horizons (only two real positive roots) . If we denote by Rh and Rc the black hole and cosmological 
horizons, respectively, the roots are 

Rn = Rh, Rc, ■ ■ ■ ,i±±,l±±, R, 

2 2 

( — - \ 

where R = — I Rh + Rc + Yli=i (7« + 7i) ) ( as * ne ro °t s must add up to zero). Determination of the 7$ 

and 7j must be numerical for each given case. In odd dimension there is an even number of roots again 
yielding (d — 1) complex horizons (only two real positive roots). With the same conventions as before, the 
roots are 



R n = Rh, Rc, -Rh, -Rc, Ji, 7i> • • • > 7^-5 , 7j-5 , 

2 2 

where again the roots add up to zero (in fact cancel in pairs). Determination of the 7$ and 7j must be 
numerical for each given case. Even without explicit analytic solutions for the roots one can proceed 
and factorize f(r) so that again the tortoise follows from integration of a rational function. Choosing 
x[r = 0] =0 one obtains 

d-l 



where k n = ^/'(-Rn) is the surface gravity at the given horizon and the sum goes through all possible 
horizons. Two special regions are of interest for our calculation: near r = one finds 



I r r d-2 

x\r] r / dr r d ~ 3 = — — , 

L J 2/j.J 2(d-2)/z' 
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just like in the pure Schwarzschild case. Near r = oo one finds 



1 f dr 1 

Observe that the constant xq G C is required so that we keep the choice of x[r = 0] =0 fixed. As 
quasinormal modes live in the region Rh < r < Rc, one can study the asymptotics of the tortoise to find 



lim x[r] = — oo, 

r\Rn 

lim x\r] = +oo. 
r/R c 

Next we turn to the master equation potentials, as listed in the previous appendix. Using the general 
formulae presented earlier it is simple to obtain the following asymptotics. Near r = one finds 

Vl (r) ~ - (d - 2? Ji_ = __L = 2_I 

*W~ 3(^-2)^ = ^ = ^, 

/i 2 1 2 - 1 

Vs(r)~-{d-2f 7MZl = -— 2 =^ r , 

just like in the pure Schwarzschild case. Near r = oo one finds 

trM 1 „ , nu2 2 d(d- 2) (d- l) 2 - 1 
V T (r) ~ yKd - 2)AV = _ ^_ ; io)2 , 

tw n 1/, „w, ,u2 2 (d-2)(d-4) (d-3) 2 -l 

H (r) ~ ±«J - 4)(d - 6)AV = V""-" = ( ?r 

v 7 4 V A 7 4(x-x ) 2 4(x-x ) 2 

The RN dS Solution: The horizons are defined via — Xr 2d ~ A + r 2d ~ & — 2fxr d ~ 3 + 2 = 0, and there is no 
general analytic solutions for the roots. In even dimension there is an even number of roots yielding (2d— 4) 
complex horizons (only three real positive roots). If we denote the black hole [outer] radius by R H , the 
inner horizon by R H and the cosmological horizon by Rc, the roots are 

Rn = R^j,Rh,Rc,11:1i, ■ • ■ ,7d-4,7d-4)-R) 

where R = — (^R H + R H + Rc + Yli=i {li + 7i)) ( as the roots must add up to zero). Determination of 
the 7j and 7 y i must be numerical for each given case. In odd dimension there is an even number of roots 
again yielding (2d — 4) complex horizons (only three real positive roots). With the same conventions as 
before, the roots are 

Rn = R H , ~Rh,Rh, -R h > r C -RcYuli, • • • ,7d-5,7d-5> 
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where again the roots add up to zero (in fact cancel in pairs). Determination of the 7$ and 7^ must be 
numerical for each given case. Even without explicit analytic solutions for the roots one can proceed 
and factorize f(r) so that again the tortoise follows from integration of a rational function. Choosing 
x[r = 0] =0 one obtains 




where k n = \f'(R n ) is the surface gravity at the given horizon and the sum goes through all possible 
horizons. Two special regions are of interest for our calculation: near r = one finds 

1 r r 2d-5 
x[r] ~— dr r 2d - 6 = — 

LJ e 2 J (2d-5)e 2 

just like in the pure RN case. Near r = 00 one finds 




The constant xo £ C is required so that we keep the choice of x[r = 0] = fixed. As quasinormal modes 
live in the region R^j < r < Rc, one can study the asymptotics of the tortoise to find 



lim x[r] = —00, 

r\R+ 

lim x\r\ = +00. 

r/R c 

Next we turn to the master equation potentials, as listed in the previous appendix. Using the general 
formulae presented earlier it is simple to obtain the following asymptotics. Near r = one finds 



4 V ,K > r ±d-w ~ 4{2d-5) 2 x 2 4x 2 ' 

■st f \ 1 M 9 w,, 19 , ^ (d-2) (5d-12) j 2 ± - 1 
W ± (r) ~ - (d - 2) (5d - 12) = 4(2 ,_ 5)V = 

lim F S± (r) ~ -\ id- 2) (3d - 8) = - («* ~ 3 (** ~ *0 _ 



,._„ - 4 V Jy > r Ad-io 4(2d-5) 2 x 2 4x 2 

_ 3d-7 

2d- 5 JV^ — 2d- 5 



where jj = j$± = and j v ± = f^jrf = 2 — jj. This is just like in the pure RN case. Near r = 00 one 
finds as before 



. . 1„, . , 2 2 d(d- 2) (d- l) 2 - 1 

t, ,n i/j «wj ,n2 2 (d-2)(d-4) - 3) 2 - 1 
F V± (r) ~ -(d - 2)(d - 4)AV = 4(g = 1 4(x J Xq)2 , 

. , 1., JW , , l2 2 (d-4)(d-6) (d-5) 2 -l 
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The Schwarzschild AdS Solution: The horizons are defined via jAjr^" 1 + r d ~ 3 — 2/i = 0, and there is no 
general analytic solutions for the roots. In even dimension there is an odd number of roots yielding (d— 1) 
complex horizons (only one real positive root). If we denote the black hole radius by Rh, the roots are 

R n = Rh, 7i> 7i> ■ ■ ■ > li^l , , 

2 2 

and again must add up to zero. Determination of the 7$ and 7 y i must be numerically for each given case. 
In odd dimension there is an even number of roots again yielding {d — 1) complex horizons (only one real 
positive root). With the same conventions as before, the roots are 

Rn = Rh, —Rh,Ji,1i, ■ ■ ■ , 7^-3 , 7d-3 , 

2 2 

where again the roots add up to zero (in fact cancel in pairs). Determination of the 7, and 7^ must be 
numerical for each given case. One case of interest is that of large Schwarzschild black holes in AdS, where 
one basically drops the 1 in the expression for f(r). For these large black holes the horizons are defined 
via |A|r d_1 — 2fx = 0, with the roots 



Rn 



2/x\ 



exp ( — - n ) , n = 0, 1, ...,d-2, 



so that one finds (d — 1) complex horizons (only one real positive root). The arguments of the horizons 
are those of the (d — 1) roots of unity. Even without an analytic solution for the roots, in the general case, 
one can proceed and factorize f(r) so that the tortoise will follow from integration of a rational function, 
as usual. Choosing x[r = 0] = we get 



x r 



, d-l 
dr 



where k n = ^f'(R n ) is the surface gravity at the given horizon and the sum goes through all possible 
horizons. Two special regions are of interest for our calculation: near r = one finds 

I r r d—2 
x i r ] ~ ~7T- / dr r d ~ 3 = --— — — — , 

just like in the pure Schwarzschild case. Near r = 00 one finds 

1 f dr 1 

x[r] ^\x\J ^ = X0 ~W 

The constant xo £ C is required so that we keep the choice of x[r = 0] = fixed. As quasinormal modes 
live in the region r > Rh, one can study the asymptotics of the tortoise to find 7 



lim x[r] = —00, 

r\R H 

lim x[r] = xq. 

r— >+oo 



7 Here one has to use a different choice of branch cuts resulting in a real value for xq. 
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Next we turn to the master equation potentials, as listed in the previous appendix. Using the general 
formulae presented earlier it is simple to obtain the following asymptotics. Near r = one finds 



r 2d— 4 ' 

A* 2 3 _ 2 2 - 1 

j,2d— 4 4x 2 4x 2 ' 

2 A* 2 1 _0 2 -l 

r 2d-4 4 X 2 4 X 2 ' 

just like in the pure Schwarzschild case. Near r = oo one finds 



W(r) ~ 3(d-2) 
y s ( r ) „ -(d-2) 



^ / \ 1 ou2 2 d(d-2) (d-l) 2 -l 

tw n 1/, „w, ,u2 2 (d-2)(d-4) (d- 3) 2 -l 

v 7 4 V n 1 4(x-x ) 2 4(x-x ) 2 

The RN AdS Solution: The horizons are defined via |A|r 2d ~ 4 + r 2d ~ 6 — 2/xr d ~ 3 + 2 = 0, and there is no 
general analytic solutions for the roots. In even dimension there is an even number of roots yielding (2d— 4) 
complex horizons (only two real positive roots). If we denote the black hole (outer) radius by R~jj and 
denote by RJj the inner horizon, the roots are 

R n = #+,#",71,7!,..., 7<i-3,7d-3> 

where again the roots must add up to zero. Determination of the 7, and 7 y i must be numerical for each 
given case. In odd dimension there is an even number of roots again yielding (2d — 4) complex horizons 
(only two real positive roots). With the same conventions as before, the roots are 

Rn = R~H, —Kh,Rhi - -Rff>7l>7lJ • • • )7d-4,7ci-4i 

where again the roots add up to zero (and in fact cancel in pairs). Determination of the ji and 7 y i must 
be numerical for each given case. Even without explicit analytic solutions for the roots one can proceed 
and factorize f(r) so that the tortoise follows from integration of a rational function, as usual. Choosing 
x[r = 0] =0 one obtains 

where k n = ^f'(R n ) is the surface gravity at the given horizon and the sum goes through all possible 
horizons. Two special regions are of interest for our calculation: near r = one finds 
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just like in the pure RN case. Near r = oo one finds 




The constant xq £ C keeps the choice of x[r = 0] =0 fixed. As quasinormal modes live in the region 
r > Rjj, one can study the asymptotics of the tortoise to find 8 



lim x[r] = — oo, 

r\R+ 

lim x[r] = xo- 



Next we turn to the master equation potentials, as listed in the previous appendix. Using the general 
formulae presented earlier it is simple to obtain the following asymptotics. Near r = one finds 



, . 1,, w , , 4 (d-2)(3d-8) jl-l 

VWD ~ - (rf - 2) ( H - 12) - 2) (M ^ _ %1 

r^o b v 7 V V^ 10 4(2d-5) 2 x 2 4x 2 

where jy = is± = an d Jv± = f^zf = 2 — jj, just like in pure RN. Near r = oo 



. » 1 , . , 2 2 d(d-2) (d- l) 2 - 1 

, , 1,, .,, jxx2 2 (d-2)(d-4) (d-3) 2 -l 
W ± (r) ~ - (i - W - 4)AV = l 4( 444 = ^4^' 

. , 1., JW , , l2 2 (d-4)(d-6) (d-5) 2 -l 



Here one has to use a different choice of branch cuts resulting in a real value for xo- 
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